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Abstract—The continuous drying of wheat in a spouted bed is idealizable as the drying of uniform, 
nearly spherical, solid particles in a well-mixed, isothermal bed, with the drying rate controlled 
by moisture diffusion within the particles and the diffusion coeflicient independent of concentra- 
tion. A rigorous mathematical solution to this general problem is developed. The non-isochronal 


diffusion equation giving the drying rate is shown to be 
M <1 — 1-04 X, exp (— 0-44.X,) 


where / | (/7i — m,)/ my — m,) is the average, reduced free moisture content of the dryer product, 
x, (SUV) \ (Dé@,), SV is the particle surface-to-volume ratio, D is the diffusion coeflicient 
and @, is the weight-average residence time in the drying bed. Diffusion coeflicients calculated 
by this equation from data on drying wheat in spouted beds agree closely with the established 
value D 207 exp (— 21960 RT). The effective surface moisture content at relative air humi- 
dities below about 40 per cent and grain temperatures above about 100° F is a constant 
m, = 0-103 Ib Ib, dry basis. The critical grain temperature above which the baking qualities of 
wheat are thermally injured is given by the relation 


t, Is 115 (m, + 7), (CF) 
where m, and /7i are respectively the initial and final moisture contents in drying. The average 
heat of desorption of moisture from wheat is estimated to be 150 B.t.u Ib higher than the normal 
heat of vaporization from a free water surface. The correlations presented facilitate the prediction 
of drying rates, energy requirements and safe and economical limits of operation for spouted 


wheat dryers. 


Résumé—Dans un lit maintenu en suspension par introduction centrale de Pair, le séchage continu 
du blé peut-étre assimile a celui de grains solides, uniformes, pratiquement spheériques dans un 
lit isotherme parfaitement agité ; la vitesse de séchage est controlée par la diffusion de Thumidité 
a Vintéricur des grains et du coeflicient de diffusion indépendant de la concentration. Etablisse- 
ment dune solution mathématique rigoureuse pour le probleme général. 

Dans le cas of les grains n’ont pas un temps de résidence uniforme, Féequation de diffusion 
exprimant la vitesse de séchage prend la forme : 

M = 1 — 1,04 XN, exp (— 0,44 X,) 


avec I] = (fi — my my, — my): titre moyen en humidité libre de produit du séchoir, X, = (S/W) 
y (DO,), SV représente le rapport surface volume, D est le coeflicient de diffusion, @,, représente 
la durée de séjour moyenne dans le lit. Les coeflicients de diffusion calculés avec cette équation 
a partir des données expérimentales pour le séchage du blé dans un lit suspendu sont en étroit 
accord avec la valeur admise D = 237 exp (— 21960 RT). Pour des humidités relatives de 
lair inféricures 4 40 pour cent et pour des températures du grain au-dessus de 100 “F, le degré 
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Vhumidité superticiclle devient constant m, =~ 0,103 Ib/Ib, produit sec, La température critique 


ou débute Paltération thermique des propriétés panfiables est donnee par 


t 18? — 115(m, + 7), ( F) 


e 
m, ct#i étant respectivement les titres en humidité avant et aprés séchage. La chaleur moyenne 
d'élimination de l humidité est supérieure d°environ 150 B.t.u/ lb a la chaleur normale de vaporisa- 
tion A partir dune surface libre d'eau. Les corrélations établies dans le mémoire permettent 
la prediction des vitesses de séchage, consommation d’énergic, les limites de sécurité et d’économie 


pour la conduite de séchoirs 4 blé, a lit suspendu, 


Zusammenfassung—Dic kontinuierliche Trocknung von Weizen in cinem Wirbelbett kann 
man idealisiert auffassen als Trocknung gleichfOrmiger, fast kugeliger Feststoffteilchen in einem 
isothermen, gut durchmischten Bett, wobei die Trocknungsgeschwindigkeit bestimmt wird 
durch die Feuchtigkeitsdiffusion in den Teilchen und wobei der Diffusionskoeflizient unabhiangig 
von der Konzentration ist. Eine mathematische Lésung wird fiir dieses allgemeine Problem 
entwickelt, Die nicht zeitlich konstante Diffusionsgleichung fiir die Trocknungsyeschwindigkcit 
wird angegeben als 

M 1 
9 — mi) der mittlere reduzierte Teuchtigkeitsgehalt des Trockengutes ist. 
xX, (S byy (DO). S70 das Verhiltnis Oberfliiche zu Volumen der Teilchen, D der Diffusions- 
koeflizient und @, die auf das Gewicht bezogene mittlere Verweilzeit im Bett. Die mit Hilfe dieser 
Gleichung aus Werten der Wetentrocknung im Wirbelbett berechneten Diffusionskoeflizienten 
stimmen gut Gberein mit der aufgestellten Beziehung D 21960 RT). Der effektive 
Feuchtigkcitsgehalt an der Oberfliche bei relativen Luftfeuchtigkeiten unter ca. 40 prozent 
und bei Korntemperaturen fiber ca. 100 °F - 0-103 Ib/lb bezogen auf Trocken- 
Dic kritische Korntemperatur, iber welcher die Backqualitaten des Weizens thermisch 


144 XN, exp (— 0-44 X,) 


wobei VV (m mi)om 


207 exp ( 


ist konstant m, 
substanz 


vermindert werden, ist gegeben durch die Beziehung 


( F) 


wobei m, und i die Anfangs- und Endfeuchtigkeitsgehalte bei der Trocknung sind. 
Die mittlere Desorptionswiirme des Wassers im Weizen wurde mit 150 B.t.u./lb hoher 


t 189 — 115 (my, + Mi), 


‘ 


gegeniiber der normalen Verdampfungswiirme von einer freien Wasseroberfliiche abgeschitzt. 
Die angegebenen Korrelationen crieichtern die Voraussage der Trocknungsgeschwindigkeiten, 
des Energiebedarfs sowie der Sicherheits- und Wirtschaftlichkeitsgrenzen in Weizen-Wirbelschicht- 
trocknern. 


INTRODUCTION 
TH 


Marucr and Gisuier [1], is a modification of 


spourED bed technique, as described by 
conventional fluidization techniques which facili- 


tates agitation of beds of relatively coarse 


particles. Instead of distributing the gas flow 
at inlet over a major portion of the bed, the gas 
is made to enter through an orifice or pipe 
the the 
The resulting pattern of agitation is a fountain- 
like turnover of the bed in which the jet of gas 
from the inlet 


centrally located in bottom of bed. 


constricted causes a stream of 


particles to jet upward through a well-defined 


core within the bed, while in the surrounding 
annulus there is a balancing downward movement 
of the bed in bulk. The advantage is gained (a 
characteristic of well-agitated beds) that particle 
temperatures are uniform throughout the bed. 
In the drying of wheat, grain temperature is the 


controlling limitation on operation [2], and this 
advantage is of decided importance. Maruur 
and Gisuier [8] therefore made an exploratory 
study of the continuous drying of wheat in a 
The 
undertaken to obtain a fundamental analysis of 
this drying process to guide the design and 
The analysis 


spouted bed. present investigation was 


operation of practical dryers. 
demonstrates the rigorous mathematical treat- 
ment of the continuous drying of solids in a well- 
mixed, isothermal bed in the case where diffusion 
within the solid governs the drying rate and the 
diffusion coefficient is independent of concentra- 


tion. 


Tue Dryer 

Fig. 1 shows the salient features of the experimental 
spouted-bed wheat dryer. The heater, C, was an insulated 
column 9 in. in diameter by 6 ft. high with a 120° conical 


bottom and a 1}in. diameter standard pipe inlet. The 
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cooler, 1), was a geometrically similar unit with a column 
diameter of 6 in. 


feeder; B, 
heater; D, 


Fic. 1. 
variable speed drive and motor; C, 


Spouted bed wheat dryer. A, 
cooler; E, pressure reducing valve; F, heat ex- 
changer ; G, H, gate valves; and I, orifice meters. 


Wueat 


Three series of drying runs were made. 
seed wheat (series 1, a No. 1 Rescue wheat; series 2, a 
No. 1 Thatcher ; and series 3, a No. 2 Thatcher) was used 
throughout to ensure that the wheat was sound, clean, 
and of high and uniform baking quality. The wheat of 
series 3 was also used in fundamental studies of moisture 
movement within the wheat kernel [4], the moisture 
desorption isotherm of wheat [5] and the correlation of the 
loaf volume of thermally damaged wheat with time, tem- 
perature and moisture content [2]. 

The mean kernel diameter of the Thatcher wheats, 
expressed as the diameter of an equi-volumed sphere and 
measured at a moisture content of 0-10 Ib /Ib, dry basis, 
was D, = 0-118 ft. The surface sphericity (shape factor) 
of Thatcher wheat, based on the visibly exposed surface 
area, has been estimated [6] as ¥ = 0-91, and may be 
assumed to be the same for all varieties of wheat. The 
kernel diameter of the Rescue wheat was 0-0110 ft. The 
distribution of kernel size in wheat of good grade is narrow, 
and is not a factor in the practical analysis of a drying 
operation [6]. 

The density of moisture-free wheat is 90-5 1b/ft® [6] 
and its heat capacity is 0-30 B.t.u./lb °F [7, 8]. Water 
held by wheat as moisture has practically the same density 
[6] and heat capacity [7] as normal liquid water. 


Registered 


Meruops 

Tempering procedure 

The moisture content of the wheat was raised to approxi- 
mately 25 per cent, dry basis, by adding a measured 
quantity of liquid water. A conditioning period of 48 hr 
at 75 °F was found adequate for the absorption and uni- 
form distribution of the added water. 
thoroughly mixed at intervals during this period. 


The wheat was 


Starting procedure 

The dryer was not emptied between runs. Operation 
was begun by introducing to the heater and cooler beds 
a flow of air slightly exceeding the minimum required for 
spouting. Feeding of wheat was commenced when the 
heater bed approached the expected steady state tempera- 
ture. 
before measurements were begun. 


The unit was then operated for three more hours 


Grain temperature measurement 

A sampling cup mounted on a rod and provided with a 
cover removable by pulling a string was used to take 
A fine 
thermocouple was immediately inserted in the sample 
and the temperature quickly taken. (Grain temperature 
must be measured immediately on removing a sample 


wheat samples from any part of the heater bed. 


from the environment under study, before thermal effects 
accompanying the equilibration of moisture within and 
between kernels cause an appreciable rise in temperature 
[5]. 
Air temperature measurement 

Local air temperatures in the heater bed were measured 
with a thermocouple probe. The cylindrical probe tip 
housing the junction was made of wire gauze, and was 


loosely packed with glass wool to provide a radiation shield. 


Humidity measurement 

Humidities were determined gravimetrically. In deter- 
mining the average humidity at the heater outlet, a stack 
of baffles was placed inside the column to ensure that the 


air sampled represented the average outlet composition. 


Moisture determination 
Moisture contents were determined by the one hour, 
130 °C oven method [9], and are reported in Ib Ib, dry 


basis, unless otherwise noted. 


Baking test 

Samples of feed and product were retained for baking 
tests. The feed sample was first air dried to 14 per cent 
moisture content, wet basis, at room temperature; both 
samples were then stored at 10 C to minimize any change 
in baking quality between the end of a drying run and the 
time at which the baking tests were done. Baking tests 
were carried out using the standard test formula with 
the addition of two p.p.m. potassium bromate, since 
the bromate formula best reveals thermal damage to 
baking quality [2]. 
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Table 1 Operating data for series 3 


L = bed depth 3 ft 


é, O13 hr 





GA Ww 


Run (ib hr) (ib hr) 


87-8 
354 vO-0 
S25 oe2 
S21 61° 
B21 68: 
821 Tt 
318 77 
318 re 
S16 
327 
B82 


$25 


342 


0, Me 71 
(Ib Ib) (Ib Ib) 
0-2580 00-1960 
0-2530 00-1810 
0-252 01560 
02433 00-1743 
00-2546 1612 
2470 O75 
2517 O-1572 
2480 01450 
2583 Ol4i4 
00-2560 1520 
02500 1572 
O2555 1582 


(hr) 


0-639 
0-591 
571 
O-SS4 
784 
0-703 
O-705 
0-608 
707 
1-100 
0-923 
0-697 





It was not known at the outset that the uniformity of 
grain temperature in the spouted heater would permit 
inlet air vreatly exceeding 


operation at temperatures 


the limits assigned to conventional dryers. As a result 
the first two series of drying runs were largely wasted 
because the conditions chosen were not severe enough to 
damage the wheat. Attention will therefore 
chiefly on the third series, in which significant damage 
occurred in five out of twelve runs. 
this series are given in Table I, 


be focused 


Operating data for 


LocaL Conprrions IN THE Srourep Hearer 


Bep 


A knowledge of local conditions in a dryer is 
a requisite for a fundamental analysis of the 
drying operation. Exhaustive observations of 
local grain and air temperatures and air humudi- 
ties in the heater bed were therefore made under 
The then 


analysed to discover their general mode and 


all operating conditions. data were 


significance. 


Local air temperatures 


Fig. 2, which is roughly applicable to all runs, 
shows that the temperature of air flowing through 
the heater bed fell very rapidly, and in the annulus 
(i.e., the bed 
surrounding moving jet 


downward moving, 
the 


* spout - 


quiescent 
central, upward 
stream or particles) approached 
thermal equilibrium with the bed at a distance 


of about 15 in. above the gas inlet. In the spout, 


the drop in temperature was stretched over a 


slightly greater distance, owing to the higher air 
velocity and sparser particle population. It is 
evident that, except in impractically shallow beds, 
air leaving a spouted wheat bed can be expected 
to be in thermal equilibrium with the bed at exit, 

The transfer of heat in a spouted bed is not 
amenable to exact analysis, owing to the impossi- 
bility of accurately determining the density of 


Fic. 2. Air isothermals in the heater bed. 


100 
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particle population and the distribution of air 
However, 
the 
integral equation of heat transfer (with a constant 


flow in the zone of heat exchange. 
data can be roughly correlated as follows : 
stream of 


transfer coeflicient) for a hot 


flowing uniformly through a deep, isothermal bed 


gas 
is easily shown to be 


haz (1) 


exp | 

e Wp 
If as in the present case this equation is not 
strictly applicable, a correlation dependent on 
the gas and particle flow patterns will still exist 
the 


t.) and az Ge, Fig. 3 shows 


between appropriately chosen values of 
groups (¢ te) (t 
such a correlation of the air temperature in the 
annulus with the height above the gas inlet. The 
differences in this relation series between 1, 2 and 
8 must be ascribed to the effects of small differ- 
ences in particle size and shape on the grain and air 


flow patterns. The data of series 3 show a long 


oz/Gc,, hr ft °F/ Btu 


Correlation of annular air temperatures in 
the heater bed. 





a* 

(Ib /hr ft?) (ft? /ft®) 
730-800 319 
850 é 208 
70 208 





*Based on an average fraction void of 0-45. 


region of linearity in which a mean heat trans- 
B.t.u_ ft® hr °F 
calculated by equation (1). The expected value 
of the from 
the correlation of Gamson et al. [10] for packed 


fer coeflicient of 2-9 can be 


true local coefficient calculated 
beds for a mean Reynolds number of 150 is 53 
B.t.u ft? hr F. The 


values cannot be explained through the gas flow 


difference between these 
pattern ; e.g. if it is assumed that the gas flow 
rate in the annulus, G,, increased linearly with 
increasing height z the inlet, from 
-0to Gate 1 ft, the differential 
heat transfer equation for the annular bed with 


above gas 


zero at 2 


this gas flow pattern 


din[(t —#,) (t; — t,)] ha 


e - 
d\inz Cy (dG, dz) 
coeltlicients which 


heat transfer are still 


an order of magnitude smaller than predicted. 


gives 


The only other explanation that could adequately 
account for the observed low rate of heat transfer 
is that the surface temperature of the grain in the 
zone of high average grain-air temperature 
differences was considerably above the average 
bed temperature. It therefore appears that in 
the zone of heat exchange of the heater bed 
sharp temperature gradients existed in wheat 
kernels in the immediate neighbourhood of their 


exposed surfaces. 


Local grain temperatures 

Fig. 4 shows the longitudinal profile of measured 
grain temperatures in the annulus of the heater 
bed. As the wheat moved downwards through 


the zone of heat exchange (Fig. 2), it absorbed 


O OtF 


Fic. 4. Longitudinal profile of grain temperatures in 


the annulus of the heater bed. 
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heat and its temperature rose to a maximum at 
the base of the bed. At this point the wheat was 
entrained by the spout and swept to the top of 
the bed, its residence in the spout being so 
fleeting as to permit no change in its thermal 
condition. 
annulus back to the zone of heat exchange, the 


The grain then journeyed down the 


heat absorbed earlier in that zone being converted 
to the evaporation of moisture in the course of 
the journey. 
the 
minimum average grain temperatures in the bed 
heat for the removal of 
moisture, its maximum magnitude is given by the 
quotient of the total heat capacity of the bed into 


Since difference between maximum and 


represents absorbed 


the heat used for evaporation per bed turnover 
evele : 
h,) (0.) 


Wy Me m) (hy, 


Al 
’ Be, 


mic) 


ha) 0, 


c ~ Hic, 0, 


“a 


("lo 


mi) (h, 


(2) 


Observed values of At, may be expected to be 
smaller than this by an amount roughly propor- 
tional to the fraction of the bed volume occupied 
by the zone of heat exchange, since the evapora- 
tion of moisture continues unabated during that 
part of the bed turnover cycle in which heat is 
being absorbed. Thus, the present experimental 
values of At, were approximately 20 per cent 
The observed values in 
all runs were in the range 24 F. 

It may be enquired how the measured tempera- 


lower than predicted, 


ture of grain taken from the zone of heat exchange 
was affected by gradients in the local temperature 
within kernels. Relevant to this question, suppose 
that a wheat kernel! at an uniform initial tempera- 
ture ¢, is suddenly put into an environment where 
its surface temperature assumes a constant value 
t,, How long will it take until the difference 

the and =the kernel 
temperatures is only 1 per cent of the initial 
difference, t, — 1,’ The diffusion equation (3), 
which is also a heat conduction equation, predicts 
a reduced time for this event of XY? — 3-75. To 
translate this figure dimensional time, 
consider the following data: the thermal con- 


between surface average 


into 


ductivity of wheat may be conservatively assumed 
to be at least 0-1 B.t.u ft hr F, (the value for 
rubber, leather, celluloid and hardwood). Average 
values of the heat capacity and density of moist 
wheat are, respectively, 0-45 B.t.u/Ib °F and 
85 lb ft®, giving for the thermal diffusivity of 
moist wheat the value 0-0026 ft? hr. The volume- 
to-surface ratio of an average wheat kernel is 
0-00175 ft. Thus the time in question is 


N?(V/S) = 3-75 (0-00175)" 
dD 0-0026 


0-00-44 hr or 16 sec 


Since the grain temperature was experimentally 
recorded 20-30 sec after removing a sample from 
the heater bed (this being the time required for 
the manipulations involved in the measurement). 
there is no doubt that temperature gradients 
within kernels effectively disappeared before the 
grain temperature was taken. There was cer- 
tainly no tendency observed for the measured 
temperature of a sample to fall with time; 
instead, it tended to rise, very slowly but still 
perceptibly, owing to the slight evolution of heat 
accompanying the gradual dissipation of moisture 
gradients [5]. The quickly measured grain 
temperatures are therefore correctly interpreted 
as weight-average values assumed upon thermal 
equilibration of samples. 

The foregoing calculation also confirms that 
significant elevations of the surface temperature 
of wheat kernels in the heater bed were certainly 
likely in the presence of large grain-air tempera- 
ture differences, but could not have persisted for 
any appreciable time outside the zone of heat 
exchange. 


Local air humidities 


The humidity of air the 
annulus of the heater bed increased at an accelera- 


flowing through 
ting rate up to about 12 in. above the gas inlet, 
and from there on varied linearly with distance 
above the inlet. Fig. 5 shows the humidity 
profile at bed outlet and shows that air from the 
spout emerged with a lower humidity than air 
from the annulus. This pattern is due to the fact 
that air in the spout moved at a higher velocity 
and contacted fewer particles than air in the 
annulus. 
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Fig. 5. Radial profile of air humidities at outlet of 


the heater bed. 


CorReELATION OF THE Dryine Rare 


The rate of drying of the wheat kernel under 
conditions in the to external 
mass transfer is small is governed by moisture 
diffusion within the kernel [4, 6, 11]. The ex- 
ternal resistance in the normal atmosphere is 
easily reduced by circulating air past the kernel, 
Reynolds numbers 


which resistance 


and is negligible even at 
DG as low as twenty-five [8}. 
isothermal drying in air of constant humidity has 
been described [4, 12] by the integral equation 
diffusion out of spheres 


for unsteady state 


subject to the conditions 

m = Me O = 0 

m = m, (constant) at the bounding surface, 
@>0 


D = constant 


More recently [6], it has been shown that the 
non-spherical shape of the wheat kernel can be 
held account of if the equation is written in the 


general form 


_s 


W > = uM exp | = = 7 xy 


At reduced drying times X* smaller than unity, 
the sum on the r.h.s. closely approximates [6] 


9 
i —,X +033 X2 (4) 
\ 7 


The rate of 


while at values of X? larger than unity the series 
The diffusion 


coeflicient has been shown to be independent of 


reduces to one significant term. 


moisture level in the range 0-13-0-25 lb Ib [4, 6] 
in which the data of three investigations [4, 6, 11] 
are represented by the Arrhenius equation 


D — 297 exp (— 21960 RT) (5) 


The effective value found for the surface moisture 
content m, when drying wheat in a stream of dry 
air was 0-103 lb Ib [4], which agrees reasonably 
well with the values of 0-099 to 0-105 lb Ib later 
found for vacuum drying [6]. 

The diffusion equations (3) and (4) provide a 
rational basis for correlating drying rates in the 
the 


reasonably 


spouted heater bed. Conditions in heater 


approached the prescribed ideal 


ck wsely. 


(i) Isothermicity 

It has been seen that the maximum variation 
in measured grain temperatures in the heater 
4 °F. 
room temperature but was, it can be confidently 


bed was only 2 Feed entered the bed at 
assumed, very quickly heated to the temperature 
of the bed, (singly exposed wheat kernels placed 
in a hot air stream approach the air temperature 
within a few minutes [4] while in the spouted 
heater heat transfer to the feed was augmented 
by direct thermal contact with the surrounding 
hot grain). Significant elevations in the surface 
temperature of grain evidently occured in the 
zone of heat exchange, but this zone comprised 
only a minor fraction of the volume of the bed. 
Altogether, it is quite reasonable to suppose 
that the weight fraction of bed material with a 
differing the 
average was negligibly small insofar as diffusional 


temperature profoundly from 
analysis of the drying process is concerned. The 
temperature dependency of the diffusion coe- 
ficient should not be expected to be a factor in 
the analysis; the coeflicient should exhibit the 
constant the 


value corresponding to average 


bed temperature. 


(ii) Boundary conditions and constancy of the 
effective surface moisture content 
An experiment was carried out to ascertain 
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the effect of the relative humidity of circulated 
air on the effective 
Fully exposed wheat kernels were dried 20 min 
F and constant relative humidity in an 


surface moisture content. 
at 77 
apparatus previously described [4], the humidity 
being regulated by carrying out the experiments 
in a controlled humidity chamber. The effective 
surface moisture content was then determined 
from the regression of the moisture loss on the 
initial moisture content [4, 6]. Fig. 6 shows 
the values so the 
relative humidity. It is seen that at high relative 


humidities the dvnamic equilibrium data merge 


obtained as a function of 


Fic. 6. Correlation of the effective surface moisture 


content of isothermally drying wheat with the relative 
water vapour pressure of the drying air. 


with the static equilibrium desorption isotherm 
previously determined [5], while at low relative 
humidities they decline to a limiting minimum 
value of 0-103 lb Ib, which is the effective value 
of the surface moisture content when drying in 
moisture-free air [4]. The data are rather closely 
represented by Smirn’s isotherm equation [13] 


, P 
w, — w' In (1 ‘ 

Po 
which was derived for the sorption of normally 
condensed water on a monolayer of strongly 
This coincidence is 
theoretically plausible, in view of the fact that 


adsorbed water molecules. 


the dynamic surface moisture content at zero 
relative humidity is closely related to the number 


SALLANS 


of primary sites available for adsorption [14], 
and it is quite reasonable to postulate that in 
general the relation between the effective surface 
moisture content in air drying and the relative 
humidity is closely approximated by merging 
the constant value m, — 0-103 prevailing at low 
relative humidities with Surru’s isotherm at high 
relative humidities (the broken curve for 122 °F 
in Fig. 6 was constructed on this basis). Relative 
humidities at the outlet of the spouted heater 
were low, being in the range 7-25 per cent, and 
Fig. 6 
indicates that one is entirely justified in assuming 


the lowest bed temperature was 100 °F. 


that in these circumstances the surface moisture 
content takes the constant value effective at low 
0-103 lb Ib. Hence the 

constant at 0 > 0 may 


relative humidities, 
boundary condition m, 
be assumed to be generally satisfied, 
(iii) Tsochronicity 

The condition of isochronal drying of wheat 
kernels in the spouted heater, necessary for 
simple application of the diffusion equations (3) 
and (4), was clearly not fulfilled. However, no 
serious difliculty is raised thereby, the mathema- 
tical form of the distribution of drying times being 
easily ascertainable as follows if it is assumed 
that mixing in a spouted bed is instantaneous 
and perfect. Suppose that the flow of feed is 
quantized in infinitely divisible slugs of dry 
mass 59, 5, being vanishingly small compared 
with the hold-up in the bed, B. If 6 is the part 
of a slug still remaining in the bed at a time @ 
after injection, then the rate of erosion of the 
di dO, Simultaneously, 
the rate of efflux of the slug in the product stream 
is Wd B= 6 @,. where 6 B is the concentration 
of slug material in the product stream. Hence, 


slug at that instant is 


by material balance, 


dé 


dg (7) 


Integration of this equation subject to the 
condition 
5=— 6, at @ 


tion 


0 gives the residence time distribu- 


exp ( ==) (8) 


r 





Drying wheat in a spouted bed 


where 6 4, is the fraction of material having a 
reduced residence time in the bed of @/@, or 
longer. 

The assumption of perfect mixing in the spouted 
heater bed was tested experimentally by injecting 
slugs of dyed wheat into the bed at the feed point 
and following the decrease in concentration in 
Fig. 7 


range of weight-average residence times, @ 


indicates that in the 
of 


the drying runs the assumption is fully justified. 


the bed with time. 


r° 


The moderate scatter of the data about the theore- 
tical curve of equation (8) is attributable to the 
vagaries of the mixing process whereby a slug 
loses its identity in the bed, and may be assumed 
to have negligible significance for a diffusional 
analysis of the drying process. 


Fic. 7. 
practical validity of the residence time distribution 


Correlation of tracer data demonstrating the 


equation. 

0,,hr O O78 O47 0-20 
To analyse the cffect of the drying time distribu- 
tion on the diffusional behaviour, recall the fate 
of a slug of feed 6, on entering the bed and 
consider the distribution of moisture content 
among successive fractions of the slug, d (8, 4,) 
in size, leaving the bed at intervals of reduced 


time d(@,80,). 


extended to time infinity the reduced free moisture 


If the period of consideration is 


content of wheat leaving the heater is seen to have 


the average value 
1 


M [ Md | >) 


. 
0 


(10) 


Substituting the diffusion equation (3) for M 
and the residence time distribution equation (8) 
for 6 do. this becomes 


(11) 


1); dd 


The integral on the r.h.s. of equation (10) can be 
evaluated graphically, using published solutions 
[15, 16] of the diffusion equation (3) or solutions 
calculated from the limiting forms discussed in 
connexion with equation (4). It was thus found 
that in the practically important range of values 
X, of 0-3 to 1-5, equation (11) is closely represented 
by the empirical equation 


M 1-04 XV, exp (— 0-44 X,) (12) 
This it may be noted is a general integral diffusion 
equation for continuous drying in a well-mixed 
bed of spherical or near-spherical particles. Com- 
parison with the equations for isochronal drying 
(e.g. equation 4) shows as should be expected 
that the distribution in residence time leads to a 
smaller average moisture reduction in a given 
weight-average residence time than if the ex- 
posure were uniform. 

The validity of the non-isochronal diffusion 
equation (12) for the spouted heater may be 
the 
diffusion coetlicients calculated thereby from the 


verified by demonstrating normality of 
drying data. Fig. 8 shows coetlicients calculated 
from the present and from Maruur and GisHLer’s 
[3] data as a function of the reciprocal of the 
average absolute grain temperature, 


T 160 


and compares these with the very accurately de- 
termined relation for isochronal vacuum drying [6] 
and with the previous data on isochronal air dry- 
ing [4]. The particle volume-to-surface ratio, VS, 
was evaluated for these calculations at the surface 
moisture content (0-103 lb Ib), which has been 
recommended for the case in which the temporal 
variation in kernel volume is small [6]. It is 
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Fie. 8. Arrhenius correlation of the moisture diffusion 

coefficient, showing a general comparison between 

values calculated from data on spouted heaters and 

previous values for isochronal air and vacuum drying. 
© Series 2 @ Series 3 

J, Becxer and SALLANs Matucer and GisuLer 


seen that the data of series 3 are in excellent 
agreement with the data on isochronal drying. 
So also are the data of series 2 and of Maruur 
and Gisuier [3], with the exception of eight 
points v .ose irregularity is explained by abnormal 
the the 
positive deviation of the two highest 


conditions in runs concerned : large 
points 
(marked with slashed squares in Fig. 8) accom- 
panies the use of shallow beds (1 and 2ft. deep) 
together with low weight-average residence times 
(0-067 and 0-140 hr), and is undoubtedly due to 
failure of the assumptions of isothermicity and 
perfect mixing under these impractical conditions, 
It is noteworthy that for other runs of Marnur 
and GisuLer [3] with weight-average residence 
times as low as 0-200 hr, but with a normal bed 
depth of 4 ft, the diffusion coetlicients are entirely 
normal. The consistent negative deviation of the 
six lowest points (those at 10° 7 > 1-77) is 
explained by the fact that the outlet relative 
humidities in these runs were in the high range 
of 50-60 per cent. 
high average relative humidity and low bed 
temperature the effective surface moistures were 


Due to the combination of 


assuredly higher than the assumed value of 


0-103 Ib/Ib, as may be judged from Fig. 6. The 
standard deviation of the normal data from the 
relation for vacuum drying is 
Altogether then the diffusion coeflicients calculated 


8:5 per cent. 


from the data on drying in spouted heaters 
correlate very satisfactorily in a normal Arrhenius 
manner with the reciprocal of the average absolute 
grain temperature, the exceptions only serving to 
emphasize the expected conditions for deviation, 
Moreover, they agree closely with the Arrhenius 
relation D = 297 exp (— 21960 RT) which has 
been accurately established by vacuum drying 
experiments [6], the application of vacuum having 
no demonstrable effect on the diffusion coetlicient 
[6]. These results conclusively prove the validity 
of the non-isochronal diffusion equation (12) and 
its underlying assumptions in interpreting the 
present (and similar) data. 

It should be noted that Fig. 8 also constitutes 
The 


standard deviation of the experimental drying 


an ultimate correlation of drying rates. 


rates from the values calculated by the conjune- 
tion of equations (5) and (12) is approximately 

5 per cent, the drying rate being more nearly 
proportional to the square root of the diffusion 
coellicient than to its absolute value. 


CorRELATION OF THE THermMaL Damace Rate 

The rate at which an organic solid can be dried 
is usually limited by the fact that its temperature 
held certain levels if 
irreversible damage to product quality is to be 


must be below critical 


avoided, Excessive thermal treatment of wheat 
leads to deterioration of its bread-making quality, 
due to thermal denaturation of the proteins. An 
earlier study [2] demonstrated that the extent 
of such damage is most clearly reflected by the 
decrease in volume of test loaves baked by the 
bromate formula. A first-order rate expression 
was found to apply : 


(13) 


The rate constant followed the Arrhenius tempera- 
ture dependency 


(14) 


——? 
RT 


A exp | 
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The frequency factor A was found to be a function 
of moisture content, and is closely approximated 
at moistures above 0-11 lb Ib by the formula 


log A 46-59 + 8-45 log m (15) 


These equations for the rate of thermal damage 
can be applied to the drying process in the 
spouted heater and, as in the case of the diffusion 
equations the analysis should take into account 
the distribution in residence time. However in 
this case the complexity of the situation, which 
also involves moisture distributions and gradients, 
is too great to permit of a fundamental analysis, 
and so equation (13) will be used directly in 
conjunction with the weight-average residence 
time, @,. This should not prove seriously objec- 
tionable, since the effects of temperature and 


moisture content ordinarily outweigh those of 


time [2]. 

Values of the rate constant were calculated for 
those runs (all in series 3) in which significant 
damage occurred, and are shown in Fig. 9 as a 
function of the absolute grain temperature at bed 


ae eo le ey 
r 4 


- 4 


O-~ tpt KBr Os 
@-2pt KBrO,y 








rez 
' on” 5 
——» °R i 
= x10 


Fic. 9. Arrhenius correlation of the rate constant for 
thermal damage to the baking properties of wheat. 
Oo -I pt K Bro, e -2ptK Bro, 


outlet, T’ = 460 +. t,. This, the maximum average 
grain temperature, was chosen in preference to 
the average bed temperature in order to assure a 


maximum safety factor in the correlations which 
follow. The curve representing the data was 
drawn for an energy of activation of 116,300 B.t.u. 
mole, the value in equation (14), and gives for the 
frequency factor in this equation the value 
40-56. This 
inserted in equation (15) to define a mean moisture 
content 


log A value of log A can be 
one weighted for thermal damage to 
baking quality—for the runs concerned, and 
gives a value of 0-193 lb lb for this parameter. 
This corresponds closely with the arithmetic 
mean moisture content in the heater, (m, + 7) /2, 
which was in the range 0-196-0-205 Ib ‘Ib. 

Fig. 10 shows an easily visualized interpretation 
of the data for all three series of drying runs. 
The group (1 
the rate constant which is practically linear with 


V,)/@, is an approximation to 


grain temperature in the region of significant 
aetna 7p 
/ 4 
mm 

/ 


0 6} 
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@-2 ot K6r0, 
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100 20 140 16c 


GRAIN TEMPERATURE, °F 


Fic. 10. Simplified correlation of damage rates in 
the spouted heater, showing the critical damage 
temperature. 


o-I1pt KBO, e -ipt K Bro, 


damage. It is clearly seen that the limiting safe 
grain temperature in the heater was 146 °F. 
At lower temperatures significant improvements 
in baking quality occurred in many runs. An 
earlier study [2] quoted critical grain temperatures 
for drying, calculated from equations (13), (14) 
and (15) for the heat treatment at constant 
moisture required to produce a 3 per cent decrease 
in loaf volume in 1 hr, (3 per cent is considered 
to be the minimum significant deviation in loaf 
volume, while at temperatures close to the critical 
drying times will not ordinarily exceed 1 hr). 
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0-ll lb Ib, these 


calculated critical temperatures are represented 


At moisture contents above 
by the formula 


t 180 230 m (16) 


The this 
formula for the arithmetic mean moisture in the 
heater of 0-20 lb lb is 1483 °F, in good accord 
with the value given by Fig. 10. One 
therefore conclude that the critical bed tempera- 


critical temperature predicted by 


may 


ture for drying wheat in a spouted heater can in 
general be estimated by writing equation (16) 
in the form 


t 180 


; (17) 


115 (my, in) 


That this equation can be safely relied on is 
vasily seen by considering the fact that in practice 
moisture reductions, m, — %, do not ordinarily 
0-10 lb Ib; the effect of 


adopting the arithmetic mean moisture content 


exceed any error in 
in equation (17) could not therefore cause an 
error greater than at most 5 F in estimating ¢,. 
To be conservative, one could in practice subtract 
5° F from the critical temperature given by the 
equation. 

The practical significance of equation (17) is 
that it expresses the limit on the rate at which 
wheat can be dried in a spouted heater if the 
breadmaking quality of the wheat is to be 
preserved unimpaired. If however the property 
which must be preserved during the process of 
drying is the viability of the wheat as seed grain, 
then an expression differing from equation (17), 
though most likely similar in form, will apply 
(the damage mechanism will still be thermal 


denaturation, but a different set of molecules 
with different Arrhenius constants will probably 
be of dominant importance). Commercial drying 
practice seems to indicate that critical tempera- 
tures for the preservation of viability are perhaps 
20°-30 °F lower than those for the preservation 


of breadmaking quality. 


Tue Exercy Baancer 

The diffusion equation (12) in conjunction with 
the Arrhenius equation (5) gives the drying rate 
of wheat in a spouted heater at any selected bed 
temperature. The energy requirements at a given 
drying rate are then determined by the energy 


balance. The total energy balance for the heater 
in steady operation is (omitting insignificant 
terms) 


Q;, 
| Wie, 
E (hyg 
W (ec, 


Mics) 
ha, 
Mg Cy) (tg — t,) 


GA (ce, 
GA (¢, 


where, by material balance, 


E = Wim, — fi) = GA(#, — #;) (19) 


The energy balance was used in the present 
analysis to caleulate experimental values of the 
average net heat of desorption, hy. The caleu- 
lated net heats ranged in value from zero to 
300 B.t.u_ Ib, 
runs of 150 B.t.u per pound of moisture removed, 


with an average value over all 
Since the contribution of the net heat toward 
the total energy requirement in drying was small, 
and since all errors in the energy balance were 
accumulated in its computation, these results 
cannot be said to be very accurate. However 
the average value (150 B.t.u Ib) is of the right 


magnitude (see [5]), gives satisfactory closure of 


the energy balance and is therefore recommended 


for all calculations on drying wheat in spouted 
beds. 


Coounc or THe Hearer Propuce 
Wheat should be stored at a temperature within 10 °F 
of the prevailing ambient temperature [17], and it is 
the object of the cooling operation to reduce the tempera- 
of the heater product to this level. It is desirable that 
the cooling effect be derived as far as possible from the 
To encourage 
this, the cooler should provide a long residence time 
(for the slow diffusion of moisture at low temperature 
levels) and should be operated at low air velocities (to 
depress heat transfer to the air). The air should emerge 
with a relative humidity of 50-60 per cent, judging from 
Fig. 6. These conditions are probably best met by a cooler 
employing crossflow of a shallow bed of wheat and a gentle 
stream of air. Matraur and Gisuiers’ data [3] on a 
counterflow cooler indicate that 50-80 per cent of the 
reduction in sensible heat of wheat in an efficient cooler 
may be utilized for the evaporation of moisture. There 
is no justification for spouting the cooler bed: uniformity 
of bed temperature, a virtue in a heater, is a defect in a 
cooler, for, since the air may be expected to leave at the 
temperature of the bed, it is evident that an exorbitant 


evaporation of moisture from the grain. 
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flow of air will be required if the cooled grain is to be no 
more than 10 °F above the ambient temperature. More- 
over, the air would have to be supplied at a needlessly 
high pressure. The spouted cooler used in the present 
study was very ineflicient in all respects and will not, 
therefore, be reported on. 


CONCLUSION 


It is outside the present scope, which has been 
limited to a fundamental analysis of the drying 
process, to discuss the prediction of air flow rates, 
pressure drops and other factors connected with 
the physical operation of spouted beds. These 
have been partially investigated by Maruur 
and Gisuier [1], and will be treated in sullicient 
detail in a forthcoming paper to permit the 
design of smoothly functioning, optimumly sized 
spouted heaters. 

Owing to its own fundamentality of approach 
and to the fact that the drying rate is moisture 
diffusion controlled, the analysis which has been 
presented of the drying process is a generally 
valid one, and is wholly independent of such 


scale-up problems as may present themselves in 


the physical operation of spouted beds, provided 
only that assumptions of isothermicity, perfect 
mixing, ete. are satisfied. It may be noted that 
there are many other polymeric materials besides 
wheat, both of natural and synthetic origin, 
whose moisture diffusion coeflicients are indepen- 
dent of concentration, and to whose continuous 
drying in a_ well-mixed, isothermal bed the 
present analysis is therefore applicable. 


APPENDICES 
Appendix A illustrates the practical use of the present 
correlations by the solution of an example problem. 
Appendix B gives a brief comparison between the merits 
of spouted and standard commercial wheat dryers. 
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NOTATION 


area of heat transfer surface in unit volume of a 
bed ft? /ft3 
cross-sectional area of a column ft? 
frequency factor hr-} 
hold-up in a bed, dry basis Ib 
heat capacity of water (liquid) B.t.u, lb °F 
heat capacity of water (vapour) B.t.u/lb °F 
heat capacity of dry air B.t.u/lb °F 
heat capacity of dry wheat B.t.u/lb °F 
diffusion coeflicient ft? /hr 
particle diameter measured as the diameter of an 
equi-volumed sphere ft 
rate of moisture removal in the heater Ib hr 
mass dry air velocity Ib hr ft? 
average mass air velocity in the annulus of a 
spouting bed Ib /hr ft? 
B.t.u hr ft? °F 
latent heat of vaporization B.t.u. Ib 


heat transfer coeflicient 


average net heat of desorption B.t.u /Ib 
rate constant hr? 
moisture content, dry basis Ib /Ib 
average moisture content of adrying particle — Ib ‘Ib 
average moisture content of a statistical 
population of drying particles Ib /Ib 
uniform moisture content prior to drying Ib Ib 
effective moisture content at the surface of 
a drying particle Ib, lb 
relative vapour pressure 
heat loss from heater to surroundings B.t.u/hr 
radial co-ordinate originating at the column axis ft 
radius of column ft 
gas constant (1-087 B.t.u mole “R) 
particle surface area ft? 
local average temperature of air or particles F 
equilibrium temperature of air and particles 
leaving the heater 
temperature of air entering the heater 
initial temperature of particles ; temperature of 
particles entering the heater F 
reference temperature F 
surface temperature of particles F 
difference between the maximum and minimum 
average particle temperatures in the heater F 
absolute temperature R 
particle volume ft 
feed rate of particles to heater, dry basis — lb/hr 
dimensionless constants in Smrru’s isotherm 
equation 
height above the gas inlet ; the axial co-ordinate 
in the heater column ft 
mass of a slug of feed, dry basis Ib 
part of 5, still remaining in the heater bed 
@ hrs after injection Ib 
local humid moisture content of air, dry basis 
Ib/Ib 
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HH, ~ average humid moisture content of air leaving 
the heater Ib /Ib 
A; = humid moisture content of air entering the heater 
Ib /Ib 
0 = time hr 
0. duration of the heater bed turnover cycle hr 
0 B/W 


4 weight-average residence time hr 


Ib ft? 
Ib ft hr 


p = air density 
fe = air viscosity 
Reduced variables in the rate equations 
M (m reduced free moisture 
content of a drying particle 
MN = (im m,)/(my - reduced free moisture 
content of a statistical population of drying 


mt.) (Me m,) 
mm.) = 


particles 
Vy ratio between the volumes of test loaves baked 
from the heater product and the feed wheat 
X = (S/V)y(De) 
isochronal drying time 
X, = (S/V)y/(De,) 
weight-average residence time 
é = 0/0, 
Y = 6V/SD, 
(shape factor) 


square root of the reduced 
square root of the reduced 


reduced residence time 
particle surface sphericity 


Appenpix A 


A problem illustrating the practical use of the correlations 
for drying wheat in a spouted heater. 

A feasible heater for practical drying consists of a 
column 2 ft in diameter with a 90° conical bottom and a 
4in. diameter air inlet. A dryer has been built which 
consists of this heater and a cross-flow cooler. The heater 
Wheat at 70 °F 
with a moisture content of 0-230 Ib Ib is to be dried to 
the maximum 0-163 Ib /Ib 
(14 per cent wet basis) and cooled to the maximum safe 


is to be operated at a bed depth of 4 ft. 


safe storage moisture of 


storage temperature of 10 °F above the atmospheric 
temperature. The kernel diameter of the wheat 
moisture of O-101b Ib is 0-0115 ft. 
temperature, pressure and relative humidity are respec- 
tively 64°F, 20in. Hg and 10 per cent. 


cooler may be assumed to be 50 per cent efficient in 


at a 
The atmospheric 


The crossflow 


converting the sensible heat of the heater product to the 
evaporation of moisture. 

Problem. ¥Find the critical safe grain temperature in 
the heater, the maximum safe drying rate, the inlet 
air temperature and the thermal energy requirements 
of the heater. Check the validity of the assumptions 
concerning boundary conditions and local conditions in 
the heater bed. 

Solution. (i) Assume that the moisture content of the 
wheat is reduced 0-010 Ib/Ib in the cooler. 


. mH = 0-163 + 0-010 = 0-173 Ib Ib 
.*. by equation (17) 
t, = 189 — 115 (my, + Mi) = 


189 — 115 (0-230 + 0-173) = 148 °F 


Check: for the maximum safe storage temperature of 
74 °F, the sensible heat removed in the cooler is approxi- 

mately 
(cy + cpm) (t, - T4) 
(0-30 


+ O17) (143 — 74) = 42 B.t.u. 


per pound of wheat, dry basis. 


.". the quantity of moisture removed in the cooler is, at 
50 per cent efficiency, 

32 16 

d (22) 0-013 Tb /Ib 

Iiyg + ha 1030 
.*. the corrected heater product moisture content is 


0-176 Ib Ib 


150 


ni 0-163 0-013 


.*. the corrected maximum safe grain temperature in the 
heater is, by equation (17), 
— 115 (0-230 


t, = 189 + 176) = 142 °F 


(ii) Assume that the effective surface moisture content 
is 0-103 Ib tb (the value obtaining at low relative humidi- 
ties). The reduced free moisture content of the heater 


product is then 
q."- 0-176 — 0-108 
oom, «= 230 — 0-108 


by equation (12), the square root of the reduced 
0-506. 
maximum 


weight-average residence time is X, 

Assume that the difference and 
minimum average grain temperatures in the heater is 
8 °F. The arithematic average absolute grain temperature 
in the heater is then 


between 


160 142 460 = 598 °R 


.. by equation (5), the diffusion coefficient is 


21960 
D = 207 exp 
RT 
21960 


28 x 10°° ft? hr 
~« 598 


207 exp | 


1-87 
The volume to surface ratio of the wheat is (the surface 
sphericity being ¥ = 0-91) 

yr DY OOS « ODT 
Ss 6 6 


0-O0O174 ft 


.. the weight-average residence time in the heater is 


- ves . 2 y (0-00174)" 
P (X, V/S) (0-506) ( 174) 0-278 hr 
, D 28 « 10 


(iii) The volume of the heater bed is 

A(L — $%) = 7(4 — 9) = lose te? 
The fraction void of a spouting bed is approximately 
0-45. 
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. the volume of wheat in the heater is 


10-44 (1 — 0-45) = 5-75 ft® 


The absolute density of the heater product is 
90-5 (1 + mH 90-5 (1 + 0-176 
ha ee. = O5 thn? 

(1 + 4537 (1 + (1-453 x 0-176) 


.*. the hold-up in the heater is 5-75 x 85 189 Ib, wet 
basis, or 
489 


(1 + 0-176) 


480 
(1 + 7) 


B 416 lb (dry basis) 


.. the feed rate is 
416 


We Be, 
(278 


1500 lb/hr (dry basis) 


(iv) The thermal energy requirements of the heater 


are 


(a) To heat the feed to heater outlet temperature 


W (ec, + Mg ey) (t, — ty) 
1500 (0-30 0-23) (142 70) 
54000 B.t.u. hr 


(b) To evaporate moisture in the heater 


W (my MH) (Neg + hg) 


1500 (0-230 0-176) (1018 150) 


94000 B.t.u/ hr 


Total 148000 B.t.u hr 


(v) A correlation on the physical operation of spouted 
beds (to be published) gives as the duration of the bed 
turnover cycle in the heater @, 0-020 hr. 

. by equation (2) the maximum difference in grain 
temperatures in the heater bed is 


i) (hyg + hig) (O-) 


W (m, 
At 
. Bie, 


} Hic,) 
84000 o-oo 


416 « 0476 


This agrees satisfactorily with the assumed value of 8 'F ; 
hence the preceding calculations need not be repeated. 


(vi) A correlation on the physical operation of spouted 
beds (to be published) gives for the minimum mass air 
580 Ib he ft®. 
Allowing a mass air flow 10 per cent in excess of this 


velocity required to spout the heater bed G 


minimum (to ensure stability of the spout) the total air 
flow to the heater is 


GA =11 K S802 2000 Ib hr 


.*. the heater inlet air temperature is 
‘ t, + (total heat from item iv) 
GAc, 
142 + 148000 


—— = 447 °F 
2000 x 0-24 


(vii) The relative humidity at heater outlet is 


Ww (my — i) 
GA FF, 


1500 (0-230 — 0-176) 2 

2000 « 0-164 
.. by Fig. 6, the assumption of an effective surface 
moisture content of 0-103 Ib/lb is well justified. 


(viii) At a fraction void of 0-45, the specified heat 
transfer surface in the heater bed is 


6(1 — €) 6(1 — 0-45) 


dD, ¥ O-O115 


200 ft? ft® 
a < 0-91 

By Fig. 2, the mean heat transfer coefficient is of the 
order of 3 B.t.u/hr ft? °F. Equation (1) then gives a 
reduced air temperature of (f — t,)/t; — t,) = 0-002 at a 


height above the gas inlet of 


< 0-24 x 2:3 


~ 


Ge (in - ‘| (1-1 & 580 
‘—#, 


ha (3 290) 


1-1 ft 


At this point t—¢, = 06°F. Hence it is certain that 


air and grain left the heater at the same temperature (,. 


Appenpix B 


A comparison of the performance of the spouted heater 
described in Appendix A with that of a standard, commercial- 
type pilot model 


The performance of the spouted dryer described in 
Appendix A may be compared with that of an experi- 
mental dryer described by SrTansriecp and CooKk*. This 
dryer was designed to simulate the performance of 
commercial wheat dryers of the kind still in general use 
in North America, (the Ellis, Hess, Morris and Randolph 
dryers). The purpose of their study was to define safe 
operating limits for commercial drying, and as a result 
of their and other investigations the maximum inlet 
air temperature in commercial drying is officially limited 
to 180 °F. 

We choose from Table 5, 
subseries A-19, a run in which the dryer was operated to 


STANSFIELD and CooKs’ 
maximum capacity at the prescribed inlet air tempera- 
ture of 180 °F without damage to the wheat, the initial 
and final moisture contents being the same as in the 
problem in Appendix A. The feed grain was at approxi- 
mately 70°F. However, the cooling air was at the low 
winter level of 9° F. It may be assumed that with air 
at this low temperature the amount of drying done in the 
cooler was negligible. Hence the moisture content was 
reduced from 0-230 to 0-163 1b/Ib in the heater alone. 
Recalculation of the quantities involved in the operation 
of the spouted heater in Appendix A for this moisture 
reduction gives t, = 147 °F, @, = 0-416 hr, W = 1000 lb/hr 


*Stansrietp E. and Cook W. B. The Drying of Wheat. 
Report No. 25. National Research Council of Canada 1932. 
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and ¢; = 395 F. 
performance of the heaters : 


The following is a comparison of the 





Spouted Standard 
Air flow to heater (Ib hr per Ib) 
wheat charged, dry basis 20 ou 
Nominal length of path followed 
by air through heater bed (ft) 
Air pressure drop through heater 
(in HO) 
Product of air pressure drop 
times air flow (proportional to 
blower power) 34 
Air inlet temperature ((F) 180 
Air outlet temperature ('F) S4 
Air outlet humidity (Ib ‘Ib) O10 
Air outlet relative humidity (°,,) 
Maximum temperature 
(CF) 
Minimum 
(°F) 
Heater outlet grain temperature 
(PF) 

Weight average residence time 
(hr) 

Enthalpy of inlet air above 70 F 
(B.t.u, 
evaporated) 


grain 


grain =temperature 


per pound moisture 
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This table permits a fair comparison of the merits of 
standard and spouted heaters. 

Although the air flow required to operate a spouted 
heater (item 1) is less than one third that for a standard 
dryer, it is not low enough to offset the difference in 
delivery pressures (item 3) and so the spouted heater 
requires more than twice as much blower power as the 
standard model (item 4). 

On the other points the spouted heater compares more 
favourably. Due to the uniformity of grain temperature, 
(items 9, 10 and 11) it can be operated at a very much 
higher average bed temperature. This increases the drying 
rate to such an extent that the weight-average residence 
time is little more than half that in the standard model. 
Due to the high inlet temperature of the air (item 5), 
the heat requirements are slightly lower (item 13) despite 
the higher outlet temperature (item 6). 

It may be noted, concerning STANSFIELD and CooKs’ 
work, that the official prescription of a maximum inlet 
air temperature rather than a maximum grain tempera- 
ture is due to the fact that, because of the absence of 
efficient mixing, the maximum grain temperature = in 
commercial heaters is rather closely associated with the 
maximum air temperature. This association does not hold 
in spouted heaters, and hence the official limit on air 
temperature should not be applied to restrict their opera- 
tion, 

The foregoing comparisons underline the salient features 
of spouted wheat dryers, and give direction for further 
evaluation of their potentialities. 
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Simultaneous heat and mass transfer in free convection 
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Abstract. The problem of simultaneous heat and mass transfer in laminar free convection from 
a vertical flat plate is investigated theoretically. An integral method is used to find solutions 
for zero wall velocity and for a mass-transfer velocity at the wall, with the wall either insulated 
or held at a constant temperature, The zero-wall-velocity solution is compared with existing 
solutions and data for simultaneous heat and mass transfer, heat transfer alone and mass transfer 


alone. 


Résumé—L auteur fait une étude théorique du probleme du transfert simultané de chaleur et 
de masse dans une convection libre laminaire a partir dune surface plane verticale. Une méthode 
intégrale est utilisée pour trouver des solutions dans le cas d'une vitesse nulle, et dans le cas 
dune vitesse finie de transfert de masse sur la paroi, celle-ci étant isolée ou maintenue a tempeér- 
ature constante. La solution du premier cas est comparée aux solutions existantes et aux données 
relatives au transfert de masse et de chaleur, simultanées, de chaleur seule, et de masse seule. 


Zusammenfassung—-Das Problem des simultanen Wiirme- und Stoffiibergangs von einer 
senkrechten flachen Platte bei laminarer freier Konvektion wird theoretisch untersucht. Ks 
wurde eine Integralmethode benutzt, um Loésungen zu finden fiir den Fall einer Wandgeschwindig- 
keit gleich Null und fiir eine Stofftransportgeschwindigkeit and der Wand, welche isoliert ist oder 
konstante Temperatur besitzt. Die Lésung fiir den Fall der Wandgeschwindigkeit gleich Null 
wird verglichen mit den bestehenden Ldsungen und Daten fiir simultanen Wiirme- und Stoffiiber- 
gang sowie fiir alleinigen Warmetransport und fiir alleinigen Stofftransport. 


INTRODUCTION Maruens et al. used an analogue computer to 
IN MANY processes mass transfer and heat solve the problem for zero wall velocity [1). 


transfer occur simultaneously. In free convection However in the process they dropped the inertia 


these may either hinder or aid one another. For terms in the momentum equation. Thus their 


example mass-transfer rates that are 100 per cent solution strictly holds only for high Prandtl and 


greater than those predicted by considering only Schmidt numbers. They also obtained data which 
the effects of mass transfer were obtained by 


Maruens et al. [1] for sublimation of spheres with 


shows that the general form of the solution is 
valid for constant Prandtl and Schmidt numbers. 


simultaneous heat transfer. Seemingly anomalous It is the purpose of this paper to derive simple 
mass-transfer results have also been observed in  Téstlts which will be expected to apply for all 
zone melting work performed by this author [2] Prandtl and Schmidt numbers and for various 


Two theoretical solutions have been carried out Situations. In part it is intended to bridge the 


gap between the work of Somers and that of 
Maruers et al, 


for this situation. Somers used an integral 
method for the problem of evaporation and 
condensation, taking into account the non-zero ; 
wall velocity arising from the mass transfer [3]. INTEGRAL EQuations 

His solution, although perfectly correct, is complex The basic equations for the approximate 
and difficult to apply to specific problems. integral solution can be shown to be [3, 4, 5]: 


*Present address: Pacific Semiconductors, Inc., Los Angeles, California. 
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“- 
Mes |S | b dy =v : (1) 
dx J ‘ 


0 
ra 


( ud dy ' (2) 


Energy 


. 
0 
4 


1 
Momentum — [ ut dy 
dr . 


Bg (T, - r.) | d 
0 
, 


ag(Cy— Cx) [ vay (3) 
0 
where the various symbols are detined in the 
tabulation of nomenclature at the end of this 
paper. 

Equation (3) includes the assumption that 
density is a linear function of both concentration 
and temperature. This will be examined in detail 
later for perfect gases. Note also that in the 
momentum equation the limits of the last two 
integrals are 5’ and 6”, not 6. This is justified by 
the following argument. The concentration and 
temperature profiles provide the driving force for 
the velocity profile. Thus we have 6 > 3, 
5 > 8”. Since é = O for y > 3’, 0 for 
y > 6” the limits on the integrals are correct as 


and uw 


shown. The only way this would not be true is 
if either « or 8 were zero. If this were true, 
however, the integral concerned would drop out 
of the momentum equation and so the limit there 
is not important. 

So that these equations may be solved it is 
necessary to assume a form for the velocity, 
concentration and temperature profiles. It is 
tacitly assumed that the shape of the velocity 
profile is independent of the ratio of the mass 
and heat transfer, while actually, of course, it is 
not. However if the Schmidt Prandtl 
numbers are not too different the result should 
be quite good anyway because the equations will 
still be satisfied on the average. 


and 


SOLUTION FOR ZERO VELOCITY AT THE WALL 


Strictly speaking the condition of zero wall 


diffusion. 


velocity is true only for equal-molal counter- 
However for low rates of mass transfer 
the wall velocity is small enough to be safely 
neglected. 

Accordingly we assume the following profiles : 


“= (s)-2() +( 
+1 -3(6)+3(8) 


(3 . 


WW 
2 \0 


ys 





If these profiles are substituted into equations 


(1), (2) and (3), the following equations result : 


Momentum 


1 ad 
105 da 


7 Ba(T, 


(6 u,*) 
T..) : s’ + ag(C, — CC.) ! &” =(5) 


Knergy 


4 ,{2)1 
~ 
105 


where 


To solve these equations we must assume a 
relationship between 6, 8° and 6”. For Se > Pr 
the temperature profile must extend farther into 
the fluid than the concentration profile, i.e. 
5’ > 6”. Furthermore the velocity must extend 
out to the end of the temperature profile, i.e. 
5 = 8’. This situation is here called “ heat- 
transfer controlling,’ and is now solved in detail 
as an example. 
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If u,, 6 and 6” are each assumed to be propor- 
tional to various powers of x and it is further 
required that the terms in # drop out when these 
are substituted into equations (5), (6) and (7) 
then we have 


1/2 é' “ 1a 
“ue=-a,ze", 6=8 ,a''*, a, a? 


where 2,, %, %3 are constants to be determined 
by substitution into equations (5), (6) and (7). 
This substitution vields 


Momentum 


Energy 


Mass 


63 4 
112 


hey 


2 as 


where = Mg Bs 


From equations (9) and (10) we obtain 


D Pr 3 | 21 105 63 4 
| = —- , 4 a 
a se 8 48 112 


which defines rvs. Pr Se. This relationship will 


(11) 


be calcuated in detail later. 

From equations (8) and (9) and the definition 
of r the value of a, is obtained. However what is 
usually desired is the Nusselt numbers, hr k 
and Kr, D. It can be shown that these definitions 
vield [3] 

Nu, 


* Nu,’ 


From the assumed velocity profiles and the form 
assumed for 5’ and 8’ we obtain 


*Primed dimensionless parameters are for mass transfer, 
unprimed for heat transfer. 


Nu 


a 


Nu,’ 
Thus fer heat-transfer controlling we have 
0-625, 34 
Nu, — 0-435 (Pr)!4 (1 —P 
Pr 


(Gr, + rGr,)''* (14) 


where 


Gr local Grashof number for heat transfer 


(= g2* B(T, — T,)/v*), 


Gr,’ = local Grashof number for mass transfer 
[= gx%a(C, —C,) v*]. 


Somers [3] and Maruers et al. [1] have found that 
r =~ (Pr Se)'?, which will also be demonstrated 
later. With this fact and by summation over the 
height of the plate we obtain the over-all Nusselt 
AL k KL D, for heat-transfer 


numbers, and 


controlling 
Nu 0-58 (Pr)! 4 [Gr — (Pr Se)! ? Gr'}!/* (15) 
pt ‘J . 
[1 + (0-625/Pr)]'/4 


— Se\ 12 
Nu Nu 
43 
where 


Gr over-all Grashof number for heat transfer 


(—gl* 8(T, — T,,)'v*), 


Gr’ over-all Grashof number for mass transfer 
[= gL? «(C, — C,,) v*]. 


In a like manner for mass-transfer controlling 


(Se Pr)! ? Gr}''* 
Se)}/4 


+ 1/8 (Ce! 
0-58 (s« ) [Gr 
[1 + (0-625 


Nu’ (16) 

These results can also be reduced from Somers’ 
solution by making the same simplifying assump- 
tions and the same definitions as made here [6]. 


A. Parabolic temperature and concentration profiles 
The constants in the equations just obtained 
are slightly different than those usually obtained 
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by the integral method [4]. Therefore a solution 
will now be obtained using the usual parabolic 
temperature and concentration profiles; i.e. let 


= (1-5). o=(1- f) 


Proceeding as before, we obtain a new defining 


(17) 


relation for r: 
Pr 


Sse 


2r+4r*] (18) 


r>[§ 


Likewise there is obtained for heat-transfer 


controlling : 
. (Pr)!/* [Gr + (Pr/Se)!? Gr}! 


Nu 0-667 
fl + (0-952/Pr))' 


Nu’ = (Se Pr)! ? Nu (19) 


and for mass-transfer controlling : 


(Se)! *[Gr' + (Se, Pr)? Gr}!/4 


fu’ <= 0-667 »/Se) ™ 
Nu . (1 + (0-952/Se)}! 4 = 


B. Comparison of results 


In form, these results approach those of 
Maruers et al. [1] as the Schmidt and Prandtl 
numbers grow large and as they approach one 
another in value. For these conditions equations 


(19) and (20) reduce to 


Nu = 0-677 (Pr)! *[Gr — (Pr Se)! *Gr'})'4 
» 

Nu’ = 0-677 (Se)! *[Gr + (Se/Pr)'* Gr}! 4 (3?) 
At this point it is interesting to compare these 
various results reduced to the case of heat transfer 
alone (Gr’' = 0) with Ostrach’s exact numerical 
solution for free convection heat transfer [7]. 
This is done in Table 1. 

Inspection of Table 1 shows that equation (15) 
is best for low Prandtl numbers, while equation 
(19) is best for high Prandtl numbers. Equation 
(21) and therefore Matures’ analogue solution is 
suitable only for relatively high Prandtl! numbers, 
It is expected that the same statements will hold 
true for simultaneous heat and mass transfer with 
* Prandtl numbers ” replaced in the foregoing by 


* Prandtl and Schmidt numbers.”’ 


Table 1. 


Comparison of theoretical results for heat 
transfer by free convection 





Calculated values of Nu/Gr' 4 
Prandtl Ostracn 


number. [7] (nu- Equation (15) Equation (19) Equation (21) 


merical) (Gr = 6) (Gr = 0) (Gr’ = 0) 


0-0766 
0-479 
0-535 
0-676 
1/103 
207 


374 


09-0650 0-0648 
0-460 0-509 
0-514 0-573 
O-O44 0-730 
1-016 1-177 
1-835 2-09 
3-26 3-81 


0-212 
0-627 
0-677 
0-806 
1-204 
214 

3-81 





Relationship between r and Pr Se 
Previously it was stated that r ~ (Pr/Se)!®, 
Table 2 shows the exact relation between r and 
PrSe as given by equations (11) and (18). 
Table 2. Relation between r and Pr/Se as given 

by equations (27 )and (35) 


(Pr Se) 


(Pr, Sey 2 


Equation Equation Equation Equation 
(11) (18) (11) (18) 
1”) 1) 1-0 10) 
0-625 0-632 0-790 O-795 
0-320 O-327 0-566 O-572 
O-114 0-118 0-337 O34 
OO170 00-0196 0-130 0-140 
oo oOo oo oo 





This table shows that r—(Pr Se)! is a 
reasonable approximation to the relation between 
r and Pr Se as given by the integral method. 
This approximation improves as Pr /Sc approaches 
1. 


D. Expansion coefficients for a perfect gas 

As stated previously, it has been assumed that 
density is a linear function of concentration and 
temperature, Over a limited range this is often 
true for liquids. For a perfect gas the exact 
relationship may be derived as follows. 
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Let C, M, 


molecular 
weight of the concerned com- 


concentration and 
ponent 
concentration and average mol- 
ecular weight of the rest of the 
mixture 

total 
C 


molecular concentration, 


Cy 


Then we have p — (CM, + C,M,) 


P 
RT 


(22) 


but (C C,) Cy (23) 





(C,, —C)(M, — M,) 


cC(M, —_ 


n 


M,) 


where 
- number of moles 


= volume 
total pressure 
ideal gas constant 


absolute temperature 


Therefore we have 
PM, 


C(M, Ta 


-M,) + 
and so 


- {1 RT.) —(/RT)| PM, 


-(P/RT) M, 





To get this into a form linear in C and T we 
must assume that the density varies only slightly 
so that in the denominator an average temperature 
T and an average concentration C can be used to 
the first approximation. With this assumption 


we obtain 
(C, —C) 
(PRT) {1 —(M, M,)) 
(T — T..) 
1— (CTR, P) {I 


‘erwi be - 


T . } 


= (26) 
(uM, M,))} 


From this we find by definition of « and § that 


C — (PRT) {1 — (M, ™,)) 


1 


e (CTR P){l 


(28) 


T {1 (M, M,))} 


In general therefore for gases or liquids x and 8 
must be evaluated at constant temperature and 


Both must be the 
value at some average value of both concentration 


and at constant concentration. 


and temperature. 


FINITE WALL-VELOCITY SOLUTIONS 
In many practical cases there is a bulk flow of 
material to or from the interface. 
by crystallization, 


This flow is 


caused such processes as 


dissolution, evaporation and condensation. The 
Cc. . « 


D Cr-C, 


T 


problem for an arbitrary flow through the bound- 
ary is extremely difficult however and has not 
been solved. As mentioned before however this 
flow may often be ignored to the first approxima- 
tion. In fact this author has used the form of 
equation (20) for correlating mass transfer results 
in zone melting [2]. Velocities through the 
boundary layer of several cm/hr were used. 

Simpler is the situation when only the com- 
ponent of interest partakes in these processes, 
i.e. it is dissolved in an inert solvent. In this case 
the bulk flow is caused by diffusion. A material 


balance at the wall vields [3] : 


| 


where C, is the total molar concentration at the 


D 
Cy 


IC 


29 
dy a 


Ve , 
0 


Cy 


wall. In such processes the wall can either be 
held at some constant temperature or it can be 


insulated. 


This is the same problem as previously solved 


Specified wall temperature 


except that v, in equations (1) and (2) is given by 
equation (29), and is not zero. It is the situation 
examined by Somers, although in a slightly 
different form [3]. Assuming the parabolic profile, 
equation (17), and proceeding as before we obtain 
for heat transfer controlling 


CaWCr — Cy). 
2r + 4r°] 
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Nu 
{0-952 {1 + (Pr 80 r Se) [(C, 


and Nu’ Nu /r 


Actually, the parabolic profile is not entirely 
the 
cannot be corrected for the effect of velocity at 


satisfactory in this case, because protiles 
the wall. This effect could conceivably be quite 
large. It would be better to use a fourth order 
polynomial for the velocity profile and third-order 
polynomials for concentration and temperature, 


B. Adiabatic wall 

Occasionally convective mass-transfer processes 
may take place at an insulated wall [8]. In these 
cases the heat necessary for the phase transition 
is supplied or removed by the bulk fluid. A heat 
balance at the wall, expressing the fact that the 
heat conducted to the wall is used up as latent 
heat, is given by 


k | “3 ) Vo Po A 


oY 


(32) 


where A is the latent heat (taken as positive for 
melting, evaporation, etc.). Substitution for v, 
from equation (29) and assuming the parabolic 


concentration and temperature profiles we obtain 


. () (Ce Ce) A 
al \Cy; — Cy) e,(T. Ty 
This together with equations (30) and (31) and 
the equilibrium relation between C, and T, is 
the solution to the adiabatic wall problem. 


(33) 


C. Example calculation 


The following example calculation illustrates 
the adiabatic wall solution. 


Assume: Water evaporating into air. 
No water vapour in bulk air. 
Bulk temperature of 120 F. 


One atmospheric pressure. 


Problem: What is the wall temperature ? 


Solution: Note that in this case we have: 


Co) = Py (P 


PP 


(Cy 
(CP 


Cy) (Cr 


C.) (Cr — Cg) 


0-677 (Pr)! 2(Gr + rGry 4 
CL) (Cr — Cy}? + Pr 


(31) 


(Pr? 80 r Se) (Cy — C,,./Cy — Cy) * 


where p, — vapour pressure of water at Ty. 


When the following physical properties are used : 


¢, = 0°25 B.t.u_ Ib °F (Perry [9)) 


A = 1060 B.t.u/lb (Krenan and Keyes [10]) 
a D — 081 (Surrwoop and Picrorp [8}) 


Equations (32) and (35) become : 


Ps PAP — Ps) | _6 (36) 


“p P, OST r pas j \r?) 


The 
temperature 7’, is assumed, and p, is obtained 
With these values r is caleu- 


The method of solution is trial and error. 


from reference 10. 
lated from equation (29). The correct result is 
recognized when / as defined by equation (36) is 
zero. Table 3 summarises the calculations for 
this problem. The result is that the calculated 
wall temperature is about 54 F, which is about 
10 °F lower than the temperature of adiabatic 


saturation given by Perry [9]. 


Table 3. Sample calculation for adiabatic wall 


f 
(p.s.La.) (p.s.ia.) (°F) 
| 


ig | 
0-915 
1-186 


1-122 


1.009 
—0°358 
0-107 
0-038 


0-295 
O-178 
214 
0-206 


44 | 56 
145 70 
145 65 
145 66 


| 





CONCLUSION 


The problem of simultaneous heat and mass 


transfer in laminar free convection has been 
solved approximately by the integral method, 
The cases of equal-molal counterdiffusion, iso- 
thermal wall and adiabatic wall have been in- 
vestigated in detail. The expected validity of the 
results improves as Pr Se approached 1. The 
form of the final results are already proving useful 


in the correlation of data for situations in which 
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heat and mass transfer are both important in Velocity of the fluid parallel to the plate 
free convection [1, 2]. in a positive a direction 
= A reference velocity 
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= Volume 

- Velocity out from the wall, at a point and 
at the wall, respectively 
Distance along the plate, from the bottom 
edge up for flow up, and from the top 

Masmere seven edge down for flow down 


= Di ‘ .. f k +4 1; 
Thermal diffusivity of liquid (k, Cpp)- SeaaE Se See aes 


Concentration of the concerned component 
at a point, in the bulk fluid, at the wall and 
and the average in the boundary layer, = Temperature densification coefficient, 
respectively. (1/p) Qp/dIT)e. 
= Concentration of the remainder of the , &, &” Velocity, temperature and concentration 
fluid at a point. boundary-layer thicknesses, respectively 


- Concentration densification coefficient 
(1 /p) (dp )--. 


= Total molar concentration of solution at a’ Ratios of thermal and of concentration 
a point or at the wall boundary-layer thicknesses to the velocity 
Heat capacity of the fluid boundary-layer thickness, respectively 
Molecular diffusivity of the concerned A = Latent heat (taken as positive for melting, 
component evaporation, etc.) 

Acceleration due to gravity (taken as 
positive for flow up the plate and negative 
for flow down) 

Local and over-all heat-transfer coefficient 
respectively 

Local and over-all mass-transfer co- % =(C — Cy )ACy — Ca) 

eflicients, respectively Gr,, Gr = Local and over-all Grashof numbers for 
Thermal conductivity of the fluid heat transfer, respectively. Defined by 
Height of the plate g2z°8(T, — T.)/*, gL* B(T, —T.)/? 
Molecular weight of the concerned com- sr’. Gr’ = Local and over-all Grashof numbers for 
ponent and average molecular weight of mass transfer, respectively. Defined by 
the rest of the fluid, respectively. gaa(Cy —Cxo), gL « (C,—C ov. 
Number of moles in a given volume Nu., Nu = Local and over-all Nusselt numbers for 
Total pressure and partial pressure at the - heat transfer, respectively. Defined by 
wall, respectively hye k and kL /k 
= Ideal gas content 

Ratio of concentration to temperature 


Kinematic viscosity 
Density of the fluid at a point and in the 
bulk fluid, respectively 

=(T — T,)((T, — Tx) 


Nu,’, Nu’ = Local and over-all Nusselt numbers for 
mass transfer, respectively. Defined by 


boundary-layer thicknesses 

2 Bit , :; K.4/D and KL/D 
Temperature at a point, at the wall, in at : 
the bulk fluid and the average tempera- 
ture in the boundary-layer, respectively Se = Schmidt number (v/D) 


Pr = Prandtl number (v/a) 
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lon exchange kinetics : the removal of oxalic acid from glycol solutions 


Cur Trexn* and Grorce Tnopos 


The Technological Institute, Northwestern University, Evanston, Ilinois, U.S.A. 
(First received 17 July 1959; in revised form 18 May 1960) 


Abstract—The kinetics for the removal of oxalic acid from glycol solutions through the use of 
the ion exchange resin, Permutit SKB, has been studied. The mechanism was found to be con- 
trolled by the combined resistances of both the liquid and resin phases. The equilibrium relation- 
ship of the system was determined experimentally and the results were found to fit a Freundlich 
type of equation. Liquid film mass transfer coeflicients were established and correlated with the 
variables of this system. Simple experimental procedures together with the diffusion equation 
enabled the determination of the diffusion coeflicient for the resin phase. This information in 
conjunction with the proper mathematical developments has made possible the prediction of 
the break-through curve which has been found to be in agreement with experimental results. 


Résumeé— L étude cinétique de lélimination de ‘lacide oxalique des solutions de glycol par la 
résine échangeur d ions, la premutite SKB, montre que le phénomeéne est contrélé par les resistances 
combinées des deux phases liquide et résine. La relation d'équilibre du systéme est déterminée 
expérimentalement et les résultats s accordent avec Péquation du type Freundlich. Les coefficients 
déchange massique de film liquide ont été établis et exprimés en fonction des paramétres du 
Des procédes éxperimentaux simples ainsi que léquation de diffusion permettent la 
Ces résultats permettent de prévoir 


svsteme., 
détermination du coeflicient de diffusion pour la phase résine 
au moyen de déeveloppements mathematiques la cassure de la courbe en accord avec les résultats 


experimentaux. 


Zusammenfassung -Die Kinetik der Abtrennung von Ovxalsiure aus Glycollésungen bei 
Benutzung cines lonenaustauscherharzes, Permutit SKB, ist untersucht worden. 
gefunden, dass der Mechanismus durch die kombinierten Widerstande von Fliissigkeits- und 


Die Gleichgewichtsbezichungen des Systems wurden experimentell 


Ks wurde 


Harzphase bestimmt wird. 
bestimmt, und die Ergebnisse erfallten cine Gleichung vom Typ der Freundlich’ schen Bezichung. 
Es wurden Stofftransportkoeflizienten fiir den Flissigkeitsfilm aufgestellt und mit den 
Variablen des Systems korreliert. Einfache Experimente erlaubten zusammen mit der Diffusions- 
gleichung die Bestimmung des Diffusionskoeflizienten fiir die Harzphase. In Verbindung mit ent- 
sprechenden mathematischen Ableitungen war es mdOglich, aus diesen Werten die Durchbruch- 
skurven vorauszuberechnen, die gut mit den experimentellen Ergebnissen tibereinstimmten. 


Tur removal of the acidic constituents formed in 
antifreeze solutions has been a subject of con- 
siderable recent interest. GLaser and THopos (1) 
considered the process involving the recovery of 
ethylene glycol by means of vacuum distillation 
after pretreating the feed stock with excess 
sodium hydroxide. Although this procedure has 
proved successful the current market price of 
ethylene glycol does not make such a process 
profitable. A direct approach without the strin- 
gent requirements associated with vacuum dis- 


tillation suggests the use of ion exchange resins 


for the removal of these acids formed in glycol 
solutions. 

The object of this investigation concerns the 
kinetic study of the removal of acid components 
from glycol solutions with the use of ion exchange 
resins. In this connexion the physical properties 
of the system pertinent to this study have been 
established the 
equilibrium relationships, liquid film mass transfer 
coeflicients and the diffusion coeflicient for the 


and involve determination of 


resin phase. The various phases of this work are 
presented in the sections that follow. 


*Present Address : Essex College, Assumption University of Windsor, Windsor, Ontario, Canada. 
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MATERIAL FOR STUDY 


For this study Permutit SKB was used. This 
ion exchange resin, furnished by the Permutit 
Company, is a strong basic resin and has been 
used primarily for the concentration of uranium. 
The resin is supplied in the chloride form and 
ranges, on an air dry basis from 14 to 25 mesh in 
size. Before use it is converted to the hydroxy] 
form according to the instructions furnished by 
the manufacturer. The resin particles were found 
to swell when saturated with the glycol solution. 

The ethylene glycol used was of chemical grade 
and was supplied by the Dow Chemical Company. 


Aqueous solutions containing 30 per cent glycol 


were prepared, to which various quantities of 


oxalic acid were added to approximate a com- 
position similar to that of antifreeze solutions. 


Eauminricum RELATIONSHIP 


Measured volumes of glycol solutions of known 
acid concentrations were added to bottles con- 
taining a weighed amount of Permutit SKB resin. 
The volumes of solutions used varied from 500 
to 1000 em’. 
sionally agitated. Small quantities of the solution 


The bottles were sealed and occa- 


were removed at convenient time intervals and 
titrated hydroxide 
solution. This procedure was repeated until the 


with a standard sodium 
acid concentration no longer changed. To reach 
equilibrium one to several days were required. 
The equilibrium relationship obtained at room 
temperature (75 F) expressed as acid concentra- 
tion in the solution, C*, vs. the acid concentration 
in the resin phase, g*, is presented in Fig. 1. 
The results presented in Fig. 1 indicate that at 
the beginning g* increases rapidly with C* and 
asymptotically approaches a constant value. 
Different empirical equations can be used to fit 
the experimental results for different ranges of 
concentration. For up to 


C* — 0-10, the following expressions give the 


concentrations 


equilibrium relationship : 


q* 3-05 0W-6C* 0.01 < CC® <— 0.05 (1) 


g* — 1088 C8 0.05 < C* <— 0-10 (2) 


It may be argued that a better representation 
for these equations is possible. Nevertheless, a 


initial Concentration 

milhequivalents /cm> 
002632 
005149 
00734 
00736 
00896 
0108 
0158 


of oxalic ocd /grom of dry resin 


q®, Acid Concentration in Resin Phase, mul 





Or Oe OS 04 OS OS SOF 08 010 
. 


C Acid Concentration in Solution, milhequivolents of oxalic acid / cm3 


Fig. 1. 
ethylene glycol solutions and Permutit SKB resin. 


Equilibrium relationship of oxalic acid in 


straight line has been selected in order to simplify 
the subsequent developments. This assumption 
seems to be justified by the close agreement 
found between predicted and experimental values 
as presented later. 

COERFFICIENT 


MAssS TRANSFER FOR 


LIQUID FILM 


A schematic diagram of the equipment used 
for the determination of mass transfer coeflicients 
is shown in Fig. 2. The main part of this unit 
consists of the differential reactor presented in 
Fig. 3. 
the lower-half being permanently fixed to a base 


The reactor is made up of two sections ; 


while the upper-half is capable of sliding into 
position and thus permitting the flow of liquid 
through the reactor. The two parts are forced 


Flow Meter 
mn 





Nitrogen — 


Integral Reactor 








o 
al 














Sample 


Fic. 2. Schematic diagram of experimental 


equipment. 
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Fic. 3. Essential features of differential 


reactor, 


together with light springs to prevent leakage 
the 
movement of the upper half of the unit is guided 
either the 
reactor and running parallel to the motion of the 
The differential reactor 
diameter of Iyyin. The bed 
between the two screens can be varied by changing 


between contacting surfaces. The sliding 


by two yrooves located on side of 


slide. has an inside 


height contined 
the amount of glass beads placed immediately 
above the fritted glass plate. 

The experimental procedure consisted of placing 
between the two screens approximately lg of 
resin in the differential reactor. This quantity of 
resin corresponds to a bed height of less than 
lem and, therefore, could be treated as a differ- 
ential reactor. With the two parts of the reactor 
disengaged, the resin was first saturated with 
acid-free glycol solution, while the upper half of 
the reactor was filled with the acid-containing 
glycol solution. By sliding the upper half of the 
reactor in line with the lower half, the reactor 
became functional and permitted the flow of the 
acid solution through the differential bed. Periodic 
samples of the effluent liquid were analysed by 
titrating with standard sodium hydroxide solu- 
tion. The total effluent was collected and used 
to establish the average flow rate. 

To utilize the information obtained from these 


and Grorce THopos 


experiments the following mathematical approach 
has been developed. The rate of transfer con- 
sidered from the liquid side can be expressed as 
Ve : ‘ 
“fT — k, ( (3) 


. 


The continuity equation is given by the expression 


ure 
d ds M 


p 


- 4 
d » *) 


Combining equations (3) and (4) and neglecting 
the second term of equation (4) produces the 
following 

u de 


— k, (C — C,) (5) 
Pp e- 


The total cumulative liquid volume, V = ust, 


and the amount of resin, m — stp, when substi- 
tuted in equation (5), produce, 
dC . 
us k, (¢ 


om 


C,) (6) 


Since the bed is considered to be differential 


IC 


om 


(7) 


where M is the total amount of resin in the dif- 
ferential bed. Combining equations (6) and (7) 
results in 

us Cy, 


8 
MiC (8) 


where (C Cy)as the arithmetic 


average of the top and bottom conditions of the 


represents 


differential bed. Of the quantities appearing in 
equation (8) C, is unknown for the general case. 
However, for the case where the diffusion through 
the resin phase is exteemely rapid as compared 
with that taking place through the liquid film, 
the concentration gradient across the particle 
becomes insignificant and C, can be obtained 
from the equilibrium relationship, or, 


C* = f(q*) 


where g* can be established from the material 
balance. Letting k,’ be defined as 


us Cy, — C 
M(C 


k,’ 
; C* day 
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An essentially constant value of k,’ over the 
entire period of a run is an indication of the 
A de- 


with time gives evidence that the 


absence of internal diffusion resistance. 
creasing k,' 
resistance due to internal diffusion cannot be 
overlooked. However, at the beginning of a run, 
since the particle is free of solute, k,’ becomes the 
same as k, and, therefore, provides a means for 
computing the true mass transfer coeflicient from 
the experimental data. 

A total of 47 runs were made for the determina- 
tion of mass transfer coeflicients for the liquid 
film. The 


varied from 0-02 to 0-31 while the concentration 


low rate as expressed in g/sec em?*® 
of the glycol solution expressed as milliequivalents 
of oxalic acid/em® ranged from 0-03 to 0-09. 
Four different particle sizes varying from 0-0338 
to 0-0578 in. have been used. 

The values of hk,’ as defined by equation (9) 


were computed and found to be decreasing with 


time. The results of Run T-1 are presented in 
Fig. 4. This is an indication that internal diffusion 
is significant. The values of k,’ thus obtained 
were plotted vs. qg, and the limiting value at 








@, Average Acid Concentration in Resin Phase 


Fic. 4. Determination of liquid film 


coellicient. 


q = 0 corresponding to ¢ = 0 was taken as the 
value of the true mass transfer coellicient for the 
liquid film. The results of the experimental runs 
are summarized in Table 1. The numerical 
values of k, ranged from 0-012 to 0-249. When 
the quantity d)G is plotted vs. k,/G@ on log—log 
co-ordinates a linear relationship results as shown 
The equation for this correlation is 


in Fig. 5. 


ky 0-0708 


G (d,G)os% 


where d, 


G mass velocity, g/sec em? 


particle size, em 


coeflicient, milli- 


dry 


k, mass transfer 


equivalents /sec g of resin 


milliequivalents /em® 


Fic. 5. Correlation of liquid film mass transfer 


coetlicients. 


In order to show that k, is a true liquid film 
coellicient values of 1/k, were plotted against 
1/G°, Since the resistance of the liquid film 
decreases with flow rate it is expected that such 
a plot should go through the origin, if the resist- 
ance due to surface reaction is insignificant. 
This fact is substantiated from the results appear- 
ing in Fig. 6. In the light of this evidence together 


with the results presented in Fig. 4 a strong 








Fic. 6. Correlation of liquid film coefficients for 
establishment of absence of resistance due to surface 
reaction. 
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Table 1. Experimental results for determination of 


: THopos 


mass transfer 


coefficients of liquid film 





Run no. 


Zz 


* 
~ 


~~ oe 


> = 


d, 
O-0859 
O-O859 
O-O859 
00-0859 


O-0859 
O-0859 
O-O859 
0-0859 
0-0859 
O-0859 
0-1468 
0-1468 
0-1468 
0-1468 
00-1468 


1264 
“1264 
“1264 
“1264 
“1264 
1264 
1264 


1049 
1049 
“1049 
1049 
00-1049 
1049 
“1049 
“1049 
1049 
“1O49 
0-1049 
01049 
00-1049 
O-1049 


O-O859 
O-O859 
O-O859 
O-O859 
O-0859 
O-O85S9 
O-O859 
O-O859 
-O859 


239 
0-601 
0-39 


0-416 
0-238 
0-593 
‘TOS 
BAG 
O52 
‘72 
“5ST 
a) 
Osis 
O-446 


0-653 
O-322 
0-248 
0-555 
O-363 
723 


213 


OOS 
523 
0-690 
215 
420 
O-348 
251 
0-101 
O-4+42 
O-322 
0-166 
0-666 
O- 208 
O97 


O-416 
0-53 
0-675 
0-604 
O722 
1-09 
202 
O54 


25 


G 


00-0304 
O-0704 
O-O496 
2205 


0-0528 
00-0303 
O-0754 
(2279 
1711 
00-1337 
O-ools 
O-0587 
O-O435 
O-OS18 
O04 


0-083 

O42 
O-OS15 
00-0706 
00-0462 
o-O919 
0-027 


00-0766 
0-0665 
O-O877 
O-0284 
O-0545 
0-442 
00-0319 
O-O128 
0-0562 
0-409 
oriz 
0-08.46 
O-O3T9 
ool 


O-0529 
O-O67T4 
O-0858 
0768 
O-O918 
O-1585 
2163 
O-O451 
OS17S8 


Co 
0-0698 
0-0698 
0-07 
0-07 


0-0302 
00-0302 
0-0300 
0-0300 
0-0300 
0-03500 
0-282 
0-0587 
0-0278 
00-0278 
00-0278 


O-0397 
0-0397 
00-0397 
0-0397 
0-0397 
0-0397 
0-0397 


00-0963 
00-0963 
00-0963 
00-0963 
00-0963 
0-963 
O-0963 
00-0963 
0-0963 
00-0963 
0-0963 
0-0963 
0-0963 
0-0963 


0-0573 
OOSTE 
0-0573 
O-O5T! 
0-057: 
0-057: 
0-0573 
O-0573 
00-0573 


ky 


0-0453 
O-113 
00-0967 
O15 


O-O915 
0-0671 
119 
O-177 
O-147 
O14 
0-0933 
0-0900 
O-O12 
0-06 
0-0671 


0-104 
O-0733 
O-O401 
O-116 
0O-0885 
0-128 
00-0491 


0-121 
0-090 
O12: 
OOSTT 
0-0933 
0-087 
O-0658 
O-O537 
0-123 
0-166 
00-0607 
O-139 
0-100 
0-0503 


O-0851 
0-008 
0-126 
0-133 
O-130 
O171 
0-202 
0-023 
O24 


d,G 
2-61 
6-05 
1-26 

18-04 


4:54 
260 
6-48 
19-58 
14-70 
11-50 
13-43 
8-62 
6-39 
4-67 


6-55 


10-49 
5-33 
3-98 
8-02 
5-84 

112 
341 


S04 
6-08 
18 


3 


hy G 


1-491 
1-609 
1-950 


1-712 
1-968 
2-063 
4195 
2-104 
4-061 
2-863 
1-643 
2-641 
4-208 


1-009 
1-454 
1-472 
1-732 
1-422 
1-235 
O32 
2447 
O-786 
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inference exists that the rate-controlling step for 
this ion exchange operation is represented by a 
combination of resistances due to the liquid and 
resin phases. 


DIFFUSION COEFFICIENTS OF THE RESIN PHASE 

To obtain the necessary experimental data that 
would permit the evaluation of the diffusion 
coeflicient for the resin phase the following 
procedure was adopted. To a designated amount 
of glycol solution containing oxalic acid, a 
specified amount of resin was introduced and 
rigorous agitation was provided with a mechanical 
stirrer. This mechanical agitation was introduced 
in order to eliminate the liquid film resistance 
surrounding the particles. For this particular 
arrangement it was found that a stirring corres- 
ponding to 1300 r.p.m. completely eliminated 
such a resistance and therefore this work was 
conducted with agitation speeds above this value. 

Periodic samples were withdrawn and analysed 
for acid content. Since the liquid film resistance 
for these runs has been eliminated the concentra- 
tion of the solute on the surface of the resin is 
in equilibrium with that in the liquid phase. 
Therefore, once the concentration of the solute 
in the liquid phase in known it becomes possible 
to account for the concentration of the solute on 
the surface of the particles with the use of the 
equilibrium relationship. This information, to- 
gether with the subsequent mathematical de- 
velopments enables the calculation of the diffusion 
coetlicient. 

The diffusion equation expressed in spherical 
co-ordinates can be used to describe the physical 
problem 


(11) 


MM » 2 M 
Hy [HF OH 
ot Pl i rw 


with the initial condition 


q = fort <0 


and the boundary condition 


. 


qj qs at b t>O0 


and 


Equation (11) can be solved provided q, is a 
known function of time. From the experimental 


data values of g, are determined at different times 
and it is possible to approximate this relationship 
with a polynomial equation of the form 


q, = 4 t (12) 


where A, B and C are constants determined from 
the experimental results. 

The solution of equation (11) satisfying the 
boundary condition expressed by equation (12) 


is found to be 


2 Ab y (—1)" 
: | mw n a n 

hrs 2 
= D(=) ] 


qi * 


exp 


in — 
6D 


2 Bb y (— 1)" 
7 Dr <a 
n 1 


t Br? a 
3 180 D- 


exp ~D (“") | sin 7 y (13) 
b b 
Letting w represent the total amount of acid 
adsorbed per particle up to time ¢ and recognizing 
the fact that the concentration in the resin phase 
is expressed on a mass basis, we have 


t 


r > 
w= | Dp (*) (4orb*) dt (14) 
rn 


Differentiating equation (13) with respect to r 
and letting r — b, substituting in equation (14) 
and carrying out the integration, the following 
expression results 
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where 


$a A 


b 


and w Dr 


It can be seen that equation (15) ts a function 


of ¢,, d, and ¢, for a specified value of ¥. The ¢’s 


are functions of time and, therefore, it is possible 
to plot tvs. q/q, with & as a parameter. On the 
the other hand values of ¢/q, obtained from the 
experimental results when plotted against time 














0 1000 500 
t, Tine in seconds 





Analysis of experimental data for determina- 
tion of diffusional coefficients. 


produce a relationship which intersects the cal- 
culated curves corresponding to different values 
of ¥. From the values of time and the parameter 
¥ corresponding to each intersection point values 
of the diffusion coetlicient, D, are computed with 
the definition of wo b?/ Dt. 
presented in Figure 7, 


These results are 

To illustrate the procedure used to calculate 
the diffusion coeflicient, the experimental data of 
Run M-13 are analysed. 


To 500 ce of glycol solution with Cy 


00-0512 
milliequivalents of oxalic acid/em*, 2-7 ¢ (dry 
basis) of Permutit SKB were added and vigorously 
agitated to eliminate the resistance of the liquid 
film. The resin when fully wet was 0.049 in. in 
The history as 
obtained from the titration of periodic samples 


diameter. time—concentration 


was as follows: 





Time 
(sec) 4/45 

0-568 
OS1T 
0-700 
1800) 0-906 
O00 1-000 
4500 “ 


500 
600 
1000 





The values of g, are obtained from the equilibrium 
relationship corresponding to the concentration 
of the liquid solution, while g, the amount of acid 


10 


——— enperimenta! 





t, Tune 9 seconds 


Fic. 8. Comparison of experimental and predicted 


break-through curve. 
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Table 2. 


Diffusion coefficients for the adsorption of ovalic acid on Permutit 


SKB resin 





Run no. dy, 


(cm) 


0-1264 
00-1049 
0-0859 
00-1264 
0-1264 


00-0376 
0-361 

0O-O178 
0-0329 
0-03382 


0-1264 
0-1264 
0-1264 
01264 
0-1264 


O-O328 
0-0825 
O-O324 
O-O515 
0-108 


adsorbed per gram of dry resin, is calculated from 
a material balance. 

With this information g, can be approximated 
with equation (12) to be 
$83 — 206 « 100% ~ 2-12 « 100° 


Ys 
Values of ¢/q, are calculated from equation (15) 
The 
results of these calculations are presented in 
Fig. 8. 
from the experimental data are also plotted to 


for different times, ¢, and parameters, yw. 
In this Fig. the values of q/q, obtained 


produce the single curve which intersects the 
calculated values of q/q, at the following points : 





t (sec) b Dt D 
10°52 
10 Spe 
10° 5p? 


482 
802 
1340 


w 


x 10°52 





Since 0-0245 in, 


D = 2-256 « 10°7 em®/sec 

In this particular case D was found to be in- 
creasing with ¢; however, in other cases, the 
opposite trend was observed. Since no definite 
trend exists a simple arithmetic average was taken 
as the true value for the diffusion coeflicient. A 
total of ten runs was made and the results are 
summarized in Table 2. For these runs, an average 
value for the diffusion coeflicient was found to 


deviation 


D % 


(em? sec) 


(10 


16°48 
2-75 
14-29 
40 


2-75 


15-93 
21-98 

O34 
24-18 
22-53 
13-46 


‘8321 ww? 


This 


order of magnitude with other systems reported 


be 1-82 « 10‘ cm? sec. value agrees in 


by other investigators [2, 4}. 


PREDICTION OF BREAK-THROUGH CURVE 


The experimental information obtained for 
this system is used in conjunction with the math- 
ematical developments presented elsewhere [3] 
to predict the break-through curve for the removal 
of oxalic acid from glycol solutions with Permutit 
SKB resin. Experimental run P-1 has been carried 
out in the integral reactor shown in Fig. 2 and 
the break-through curve resulting for it has been 
compared with the predicted curve. The pertinent 
experimental data are as follows : 


Column diameter 1-625 in. 


Bed height 
Amount of resin 
Particle size 
Flow rate 


Concentration of influent solution ; Cy 


5-0 in. 

68-9 g (dry basis) 
0.0413 in. (fully wet) 
0-976 cm?* /see 


= 0-0461 


milliequivalents /em*®. The concentration history 


of the effluent stream was found to be: 





Time (sec) 


1000 
2000 
3000 
4000 
5M) 
6600 


00-0289 
0-245 
0-45 
0-65 
0-762 
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the initial concentration C, = 0-0461 


equation (1) applies under which conditions the 


Since 


equilibrium relationship is 


q* = 3-05 + 40-625 C* (1) 


It has been pointed out elsewhere [3] that the 
explicit expression for the concentration ratio 
depends on the form of the equilibrium relation- 
ship. For this particular system, the results of 
the general linear case [3] can be applied. This 
expression has been found to be 

A, 

KC 


1 | P(E, r) 


| ae 


(16) 


Numerical values for ®(€, +) are available in the 
literature [3]. Using equation (16) concentration 
The 


have been calculated. procedure 


involved is given as follows: 


ratios 


3-05 
Therefore, the parameter 


From the equilibrium relationship A, 
and A, = 40-625. 
K, _ 3-05 . 1 632 

KC, 40-625 © 0.0461 


For this run k, is calculated to be 0-0982 milli- 


equivalents of oxalic acid transferred sec gram of 


dry resin milliequivalent /em*, 
It has already been pointed out [4] that ®(£, + 
is related to the parameters 


BDK,p | 


N— 3D Ky 
bk, 


and 


Substituting the appropriate values : 


_ BCLS? & 107) (40-625 
N 0-082 
(0°05245)* (0-0082) 
10°*) (40.625) 
(0-05245)* 


3 (1-82 . 
8-003 10° 


Spr 68-9 ae 
= _ 70-99 
wa 
su wore 


(8-063 » 10°) (70.59) — 0.5602 
2D 10") 
bs (005245)? 


Adz 


2 (1-82 


1.3231 10 * 


The quantity A é=z represents the void volume 
of the bed which is found to be 44.em*. Therefore 


44 


= 45-1 
0.976 


Therefore +r 1-3231 wt 


oT 


1-3231 (17) 


ort 
Equation (17) relates the parameter or with the 
time, ¢. Using the solutions of ®(£, r), together 
with equation (16), the following concentration 


ratios have been calculated for different times : 





{ (sec) 


O-0000 
1578 
Oa052 
4262 
OSs21 
6651 
(7658 
(-SZS0 


2510 
2648 
wee 
bee 
S822 
4578 
A 
GOS 
6845 
Tuo 


OSTIO 
0-9000 





The predicted values of C/Cy are compared 
The 


made its 


with the experimental results in Fig. 8. 


experimental break-through curve 
appearance earlier than predicted, but its slope 
increased at a slower rate. This agreement is quite 
favourable and particularly since the correlations 
used to obtain the mass transfer coetlicient and 
the equilibrium relationship as well as the average 
value of the diffusion coetlicient can vary as much 


as 10 per cent. 


NOTATION 


constants for equation (12) 
radius of particle em 
concentration of solution milliequivalentscm® 
concentration of influent solution 
milliequivalents /em4 
concentration of solution at surface 
milliequivalents /em® 
equilibrium concentration in liquid phase 
milliequivalents of oxalic 
acid ‘cm® of solution 
particle diameter cm 
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diffusion coefficient in resin phase cm? /sec 
mass velocity g/sec cm* 
mass transfer coeflicient for liquid film milli- 
equivalents sec gram of dry resin 
milliequivalent /em® 
pseudo-mass transfer coeflicient for liquid film, 
milliequivalents, sec gram of dry resin 
milliequivalent /em? 
3DK, bh 
amount of resin in bed 
amount of resin present in differential bed 
average concentration in resin 
milliequivalents yg dry resin 
equilibrium concentration in resin phase 
milliequivalents/g dry resin 
point concentration in resin 
millequivalents, g dry resin 


concentration on resin surface 
milliequivalents /g dry resin 

cross-sectional area of reactor cm? 

time sec 

superticial velocity em //sec 

amount of acid adsorbed in resin 

milliequivalents / particle 

bed height em 

3DK, bv? 

pz/u 

bulk density or resin g dry resin/cm* 

2D /b 

t — oz/u 

porosity 

A/(A + Bt 

Bt/(A + Bt 

CPR/(A + Bt 

b?/ Dt 
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Heat transfer to N,O, gas flowing normal to a heated cylinder* 
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Abstract. The NO, NO, gas system is convenient for experimental studies at relatively low 
temperatures of heat transfer processes in equilibrium-type reacting systems. The thermal 
conductivity and specific heat of these systems are unusually large and are non-uniform functions 
of temperature making difficult the assignment in heat transfer correlations of effective film 
values of these thermal properties. Rates of forced-convective heat transfer from a cylinder to 
N,O,- NO, gas at 1 atm were experimentally determined at free stream velocities up to about 
40 ft/sec. If specific enthalpy instead of temperature is used as the driving force for heat transport 
a consistent correlation of the over-all heat-transfer rate data can be obtained by dimensional 
analysis, employing equilibrium properties in the conventional relationship, 

Nyy = (Nee. Nae) 
The advantages of using specific enthalpy as a driving force are justified from theoretical consid- 
erations. Due primarily to the high thermal conductivity of this system, values of conventional 
heat transfer coeflicients based on a temperature driving force may be as high as eight times 
those for a similar non-reacting system. 
Résumé —Le mélange gazeux NO,-N,O,, est utilisable pour des études d’échanges de chaleur 
des réactions d’équilibre 4 des températures relativement basses. Comme la conductivité therm- 
ique et la chaleur spécifique de ces systémes sont anormalement élevées, et ont une variation non 
uniforme en fonction de la température, il est difficile de relier la valeur effective du film a ses 
propriétés thermiques. Les vitesses d’échange de chaleur en convection forcée d'un cylindre de 
NO, gazeux sons une atmosphére ont été déterminées expérimentalement pour des vitesses de 
la veine gazeuse atteignant prés de 40 ft/sec. 

Si le potentiel pour ’échange de chaleur est lenthalpie spécifique au lieu de la température, 
on obtient une relation pour la vitesse de l'échange de chaleur global en se servant de l'analyse 
dimensionnelle, les propriétés d equilibre étant sous la forme conventionnelle : 

Nnu = ¢(Npr, Nre) 

L’emploi de lenthalpie spécifique comme potentiel se justifie d'un point de vue théorique. 
\ cause de la valeur élevée de la conductivité thermique de ce systéme, les valeurs habituelles 
des coefficients d’échange basés sur la température sont jusqu’é huit fois plus grands que ceux 
des systémes stables analogues. 

Zusammenfassung—Das gasfOrmige System N,O, — NO, ist fiir experimentelle Studien bei 
relativ tiefen Temperaturen von Wirmeiibergangsvorgingen bei Gleichgewichsreaktionen gut 
geeignet. 

Die Wiarmeleitfahigkeit und die spezifische Wirme dieses Systems sind ungewOhnlich gross 
und sind uneinheitliche Funktionen der Temperatur, wodurch die Zuordnung der effektiven 
Filmwerte dieser thermischen Eigenschaften bei der Korrelation des Wirmeiibergangs schwierig 
wird. 

Die Geschwindigkeit des Warmetibergangs bei erzwungener Konvektion von einem Zylinder 
in ein N,Og — NO, — Gas bei 1 atm wurden experimentell bestimmt bei freien Damfgesch- 
windigkeiten bis zu ca. 40 ft/sec. Wenn man an Stelle der Temperatur die spezifische Enthalpic 
als treibende Kraft fiir den Wiirmetransport benutzt, erhilt man mit Hilfe der Dimensionsanlayse 
eine Korrelation der Gesamt-Warmeiibergangsgeschwindigkeitsdaten, wobei in der konventionellen 
Beziehung 

Nnu = ¢(Npp Nre) 


*Presented at the American Institute of Chemical Engineers Annual Meeting, San Francisco, December 1959. 
+Present address: Convair Astronautics, San Diego, California, U.S.A. 
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Gleichgewichtseigenschaften benutzt wurden. 


Die Vorteile bei der Benutzung der spezifischen Enthaplie als treibende Kraft ergeben sich 


aus theoretischen Betrachtungen. 


Besonders wegen der hohen Wiirmeleitfihigkeit des Systems kinnen die Werte fiir die 
konventionellen Wirmetibergangskoeflizienten bezogen auf die Temperatur als treibende Kraft 
bis acht mal so hoch sein als bei ahnlichen nichtreagierenden Systemen. 


INTRODUCTION 


Tue detailed and theoretical experimental study 
of heat transfer in reacting systems has been of 
great interest to investigators in many different 
fields in recent years. A large part of the impetus 
for this study has arisen from the modern develop- 
ment of high-speed missiles, where heat transfer 
processes in rocket engines and re-entry aero- 
dynamic heating involve chemical dissociation. 
The problem is also of interest in the field of com- 
bustion and industrial chemical processes. In 
1904, Nernst [1] derived a theoretical expression 
for the thermal conductivity of a dissociating gas. 
Recently more exact expressions have been ob- 
tained by Hirscurevper [2], Franck [3], Frank 
and Srpavruorr [4], Haase [5] and Burier and 
Brokaw [6]. Studies of the internal problem of 
heat transfer to N,O, gas flowing at subsonic 
velocities in a tube have been made by Bra and 
Lyer.y [7], Scnorre [8] and Krieve and Mason 
[9]. Brokaw [10] has analysed the data of Bra 
and Lyerty and recommended cor- 
relations. Rose and Srark [11] have measured 
stagnation point heat transfer in dissociated air 
at velocities above sonic using a shock tube, and 
their results appear to support the theoretical 
analysis of Fay and Ripper. [12]. Rosner [13] 
has summarized much of the available theoretical 
and experimental work applied to the problem 
of dissociation accompanying aerodynamic heat- 
ing. 


several 


SENSITIVITY OF EQUILIBRIUM THERMAL 
PROPERTIES TO TEMPERATURE 

Equilibrium thermal properties are not only 
unusually large but are non-uniform functions of 
temperature making difficult the assignment of 
effective values to properties of a fluid surrounding 
a heated body. Experimental values of isobaric 
heat capacity and thermal conductivity of the 
N,0, system are shown in Fig. 1 from which this 
temperature dependence is evident. In subsequent 


discussion properties are assumed to be equilib- 
rium properties unless otherwise specified, and 
the term N,O, will refer to an equilibrium mixture 
of N,O, and NO, in accordance with 


N,O, <> 2NO, 


AH, = 2-47 x 10° B. Th. U. at 77°F (1a) 


All properties used in this paper are those compiled 
by Brokaw [10] for a total pressure of 1 atm. 
A system is designated as an equilibrium-type 
where chemical equilibrium is established so 
rapidly that reaction kinetics do affect 
processes occuring in the system. The reaction 
of equation (la) is in this category for subsonic 
flow processes as supported by the data of Baver 
and Gustarson [14] and Carrincron and 
Davipson [15] which indicate that equilibrium 
near room temperature is established in the order 
of a «sec. If heat is added isobarically and with 
no heat loss to a closed endothermic system, the 
temperature rise is smaller than would be pre- 
dicted from a weighted average of the heat 
capacities of the individual species of the mixture 
since a portion of the heat is absorbed as enthalpy 
of reaction. Thus the measured equilibrium iso- 
baric heat capacity C,, is larger than the sum of 
the weighted heat capacities of the individual 
species, Cp, the so-called frozen heat capacity. 
Similarly, the equilibrium thermal conductivity 
k, of an endothermically reacting system is 
greater than that of the same system frozen, kg, 
due to the added contribution of the diffusion of 
enthalpy of reaction. Thus a reacting mixture, 
flowing past a heat source, absorbs more energy 
per unit time than does a frozen system under 


not 


otherwise similar conditions. 
For a general endothermic equilibrium dissocia- 
tion of the type: 

X, 25 2X (1b) 
the isobaric heat capacity C, may be expressed 
as the sum of a frozen and reaction component 
as shown by Pricocine [16]: 
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+ Cr, (2a) 


AH, d€ (2b) 


"a Mx, d T Pp 


For systems in chemical equilibrium, the term 
(9€/dT), may be utilizing the 
van't Hoff relationship [17] giving: 


evaluated by 


€(1 — €) ; AH,\? 


‘ 2RMy, T (2c) 


Cp = Cp 
Depending on the magnitude of AH, the reaction 
component, Cp, may be much greater than the 
frozen component, Cp,, resulting in large values of 
equilibrium specific heats. 

Similarly the thermal conductivity, &, on the 
assumption of additivity of a frozen and reaction 
contributien and neglecting thermal and pressure 
diffusion may be expressed as described in Refs. 
[4] and [5]: 

k=k, +k (3a) 
where for the equilibrium-type system of equation 
(1b) kg is shown in Ref. [5] to be equivalent to 
Dy, p Crg, and thus 


k =k, (3b) 


. 


Dy, p rg 


Within the precision of the existing experimental 
thermal conductivity data of Franck and 
Spavrnorr [4] and Corrix and O'Neat [18] this 
assumption of addivity appears to be valid. 


Insensivity oF k/C, TO TEMPERATURE 

It will be shown that the ratio k/C, is relatively 
insensitive to temperature compared to either of 
these individual properties. This fact will prove 
to be of convenience in correlating cylinder heat 
transfer data where by using specific enthalpy 
instead of temperature as a driving force the ratio 
k/C, instead of k alone appears in the Nusselt 
number, making the corresponding heat transfer 
film coefficient much less sensitive to temperature 
than the coefficient based on a_ temperature 
driving force. By dividing equation (8b) by 
equation (2a), an expression for k/C,, results : 


— he + Dine Cog _ 
Cy Ong + Cog 
(k Cy) ¢ [(Crg (pg) + (Djs p Cog/k,)] 
((Cpg/Crg) + 1) 
Dy, p Cog 
ky 


(4a) 


Nie, ( th) 


Hixscuretper et al. [19] and Burier and 
Brokaw [6] show that for gaseous systems of the 
type of equation (Ib) Nte, is approximately 
unity and thus from equation (4a) and (4b) 


(4c) 


k k 
c= (eG), Nr 


Since calculations based on the kinetic theory 
[19] show that for gases Npr, is not a marked 
function of temperature and , varies approx- 
imately as 7? ?, equation (4c) indicates that k/C,, 
might be relatively insensitive to temperature. 
This conclusion is substantiated by data for the 
N,0, system shown in Fig. 2. 
CONSTANT SPECIFIC ENTHALPY GRADIENT 
Another result of the relative insensitivity of 
the ratio k/C,, to temperature is that the specific 
enthalpy gradient in a thin laminar boundary 
layer surrounding a cylinder will be shown to be 
approximately constant at a given position 
along the leading portion of the cylinder. This 
fact further enables one to assign appropriate 
effective film properties of the fluid by basing 
them on a state corresponding to a linear average 
of the free-stream and surface specific enthalpy. 
One may write for steady-state two-dimensional 
heat flux in the laminar boundary layer on a 
cylinder : 
«Sage 


A dr (5a) 


or for an isobaric, equilibrium process where 
dH = C,dt 

Q | k\ dH 

A Cp) dr 
Substituting & from equation (3b) into equation 


(5b) and eliminating Cp, by equation (2a) 


Q@ 1 — 
~ = Ge [he +e DualCr— Co] S 


(5b) 


(5c) 
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Substituting equation (4b) in (5c) and rearranging 


@ hk ) ’ _ ke . dH 
: —-}| Nie, + — (1 Nie _—- 6a 
A (c é ‘ Cp ( ? dr ( 
Since for gaseous systems of the type of equation 
(1b), Nueg ~ 1, equation (6a) becomes 


Q h ) dH 

on “06 Real aime 

Af (c ¢ dr 

It has been noted in the previous Section 
(k/Cp), for gases is not markedly a function of 
temperature, and thus for the case of steady-state 
heat transfer through a cylindrical laminar film 
so thin that A in equation (6b) is effectively 
constant, the specific enthalpy profile would be 
expected to be approximately linear throughout 
the film. That the laminar momentum boundary 
layer on a cylinder is relatively thin is verified 
from ScuiictinG [20] where for Nre, > 1500 the 
boundary layer thickness is less than 9 per cent 
of the radius of the cylinder. Shown in Fig. 3, 
assuming that the thicknesses of the thermal 
and momentum boundary layers are the same, 
are three temperature and corresponding specific 
enthalpy profiles for the N,O, system calculated 
from equations (5a) and (5b) respectively for 


(6b) 


average thermal boundary layer thicknesses of 


0-0094 in. (NRe, = 3220), 0-012 in. (NRer = 1980) 
and 0-014 in. (Nre, = 1470). It is evident from 
Fig. 3, as expected from equation (6b), that the 
specific enthalpy profile is quite linear, whereas, 
as might qualitatively be predicted from equation 
(5a) used in conjunction with the conductivity 
data in Fig. 1, the temperature profile is curved. 
Thus both thermal and non-thermal properties 
based on linearly averaged specific enthalpy 


should approximate the effective film state and 


lead to consistent heat transfer correlations, 
whereas properties based on linearly averaged 


temperature : 


(Tb) 


would be less accurate. 


HEAT FLUX THROUGH LAMINAR BOUNDARY 
LAYER 


For separated flow such as occurs around a 
cylinder for Reynolds numbers greater than the 
order of unity and less than 10° where the 
boundary layer becomes turbulent, it is necessary 
to consider the relative magnitudes of heat 
transfer through the laminar boundary layer and 
through the separated region. DouGLas and 
CnurcuiLt [21] have proposed a method of 
correlating heat transfer to cylinders for non- 
reacting systems by considering the heat transfer 
coefficient to consist of two thermal resistances 
in parallel: a laminar boundary layer in the front 
of the cylinder and a flow region beyond the point 
of separation. Even at the maximum Reynolds 
number of 18,000 


encountered in the present set of measurements, 


approximately which was 
the laminar film portion contributes 73 per cent 
to the film coetlicient calculated from equation 
(4) of Ref. [21]. Thus in the present investigation 
the heat flux through the laminar boundary 
layer predominates. At higher Reynolds number 
the separated region contributes more to the 
relative heat flux from the cylinder, with turbulent 
transport becoming more predominant; thus the 
sensitivity of correlations to the film properties 
might be expected to diminish. 


Heat 
SPECIFIC 


TRANSFER CORRELATION USING 


ENTHALPY AS DRIVING FORCE 
Cylinder heat transfer coeflicients may be de- 
fined either in terms of temperature or specific 


enthalpy as a driving force in accordance with 


h (t, h At (8) 


h’ AH 


t.) 
oh’ (H, — H,) 


@ 
{ 


Using temperature or specific enthalpy as driving 
forces together with average film properties based 
on an average between the surface and free stream 
temperature, or an average between the surface 
and free stream specific enthalpy respectively, 
and applying the techniques of dimensional 
analysis as described by McApams [22], one 
obtains respectively, the following relationships : 
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Pe [(PaL) (4) 


N Pry ) 


(Ce, 


Npry) 


ss = 


NNy = — > 


(9a) 


“e | 


(9b) 


= ¢ fh Rey, 


at, _* 
= d’ (N Rey, 


For given values of Ny, and Np», the value of 
Nx, is established. With temperature as a driving 
force the heat transfer coefficient, h, is directly 
proportional to ky (equation 9a), which has been 
seen to be temperature sensitive making the 
assignment of the effective state of the film quite 
critical, and _ particularly the inherent 
curvature of the temperature profiles (Fig. 3) 
precludes a simple linear average for effective 
film temperature (equation 7b). On the other 


since 


hand with specific enthalpy as a driving force the 
heat transfer coeflicient, h’, 
(k/C,)y equation (9b) which has been seen to be 


is proportional to 
not very temperature sensitive, making the 
assignment of the effective state of the film much 
less critical. Also, of the 
stancy in the film of specifie enthalpy gradients 
as compared to temperature gradients, properties 


because greater con- 


such as » and p and also Ny, and Ny», can be 
assigned precise values based on an average film 
specific enthalpy (equation 7a) since these pro- 
perties are more nearly linear with specific 
enthalpy than with temperature. 


EQUIPMENT AND METHODS 

These experiments were performed in a closed- 
cycle flow system shown schematically in Fig. 4 
The working section consisted of a two-dimen- 
sional flow channel, formed by two parallel 
aluminum plates 24 in. wide, 7 ft. long and sep- 
arated by a distance of 2 in. The maximum velo- 
city encountered in this set of measurements in 
the channel was 37-8 ft/sec. The intensity of 
turbulence measured with a hot wire anemometer 
set described by Fuqua [23] was less than 0-52 
per cent in air flowing at 33 ft/sec. just upstream 
of the cylinder. Bulk velocity was measured by 
means of a calibrated Venturi meter, the velocity 
being adjusted to the desired value by control of 
the blower speed. 


Boynton, 
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The cylinder used in these studies was a modi- 
fied General Electric Calrod Tubular Heater with 
an internal heater insulated with ceramic material 
from an Inconel sheath, 0-315in. in outside 
diameter, 0-255 in. in inner diameter and 18 in. 
in length. A pair of 36 B. and S. gauge, Teflon- 
insulated iron-constantan thermocouple wires 
were placed in each of two longitudinal grooves, 
approximately ,/; in. deep, machined on diamet- 
rically opposite sides of the sheath. Each of the 
two grooves extended from the opposite ends of 
the sheath to about lin. from the midpoint 
where the thermocouple junctions were spot- 
welded directly to the surface of the Inconel 
sheath. The wires were cemented in the grooves 
with Sauereisen low expansion cement which 
was polished flush with the cylinder surface in 
order to eliminate protuberances. The heated 
cylinder was inserted near the top of the channel 
normal to the flow of gas and midway between 
the channel walls as depicted in Fig. 4. At this 
position the maximum boundary layer thickness 
on the channel walls encountered in these experi- 
ments was less than 10 per cent of the distance 
between the walls. Thus the effect of the channel 
walls boundary layers on the cylinder was negli- 
gible. The cylinder could be rotated through 
180° permitting complete circumferential tem- 
peratures to be obtained from the two thermo- 
couples positioned 180° from one another. For 
each free stream gas velocity and heat input, the 
average cylinder surface temperature was deter- 
mined by averaging, arithmetically, the tempera- 
tures measured at 30° intervals around the circum- 
ference of the cylinder. An accurate measurement 
of the d.c. power delivered to the cylinder was 
obtained by means of a simple voltage divider 
network. 

Radiation losses from the polished Inconel 
cylinder surface to the polished aluminum walls 
of the channel section were calculated and found 
to be less than one per cent of the total heat 
flux and were therefore neglected. As has been 
mentioned previously in the literature [8], 
natural convection effects can be appreciable in 
reacting systems due to abnormally large co- 
of thermal The natural 
convection loss in the N,O, system will be 


efficients expansion. 
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relatively large at low Reynolds numbers due to 
its large Grashof number. Crude calculations for 
test no. 3, using N,O, equilibrium properties 
(based on T;) in the usual natural convection 
correlations for non-reacting systems, indicated 
that natural convection may account for as much 
as 10 per cent of the cylinder heat flux at this 
low Reynolds number, but because of the un- 
certainty in these calculations for reacting sys- 
tems, no corrections were made. The effect on 
U of the cylinder blockage of flow in the two- 
dimensional channel was in all cases less than 
2-6 per cent and was therefore neglected. 

Nitrogen tetroxide supplied by the Matheson 
Company was used in these studies. The purity 
of the gas was checked by noting that a frozen 
sample was straw-white indicating a lack of 
impurities such as N,O,, which imparts a bluish 
color to the solid. Before introduction to the 
tunnel the nitrogen tetroxide was passed through 
a tank containing solid drying agents. Gas density 
measurements performed between runs indicated 
that as much as 15 vol. % of air may have been 
present in the system. Weighted on a molal 
basis, properties of the system were insignificantly 
affected by this quantity of air. 


RESULTS AND DISCUSSION OF DATA 

As shown in Tables 1-3, sixteen tests were 
made, covering a range of free stream gas veloc- 
8-86 to 37-8 ft/sec. Corresponding 
average cylinder film Reynolds numbers ranged 
from about 1,300 to 19,000 and over all average 
Nusselt numbers from about 30 to 80, exact 
values depending upon whether temperature 
driving force and properties based on the average 
temperature (equation 7b) or specific enthalpy 
driving force and properties based on the average 
specific enthalpy were used. As shown in Table 1, 
heat flux densities (equation 7a) ranged from 
2,203 to 23,050 B.Th.u. /hr ft®?, producing tempera- 
ture differences from 12-6 
the cylinder surface and the free stream. It should 


ities from 


F to 553 °F between 


be noted that these large temperature differences 
resulted in large, non-uniform variation of thermal 
conductivity in the film (Fig. 1). Free stream 
temperatures were maintained constant during 
each test at about 75°F. That the cylinder 
surface temperature approached isothermal con- 
ditions is indicated in Table 1 where the maximum 
difference in circumferential temperature, Atcna, 
is presented. The largest value of Ate,,,,,/ts, 0-25, 
occurred in Test no. 8. 


Table 1. Summary of experimental heat transfer data 





| 


' t ts 
B.Th.U. 
hr ft? 


/ 


(°F) 





14-950 
2-208 
18-580 
2-208 
15-080 
18-810 
8-769 
22-910 
2-203 
14-900 
18-830 
23-050 
8-593 
14-780 
22-950 
19-040 


H, Hs Hy 


(ee) (aE) (=) 


47°5 
52-5 
50°3 
55°3 
48:1 
51-2 
52-0 
51°83 
56-1 
48:1 
51-4 
51°8 
58-1 
48-1 
51-6 
50-0 


t 
Cmax 


(°F) 











199 
69-3 
225 
65-9 
138 
158 
100 


120 


i 
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Table 2. Value of properties at 1 atm employed 





ky x 108 Cp, (u/p)y x 104 ' ky » 108 Coy 
aoe So ‘) (") — ) (am. Th.t 
see ft “F Ib F see . (— ft oF 


6-50 MH 143 i 21-5 
19-5 . 0-448 2 19-6 

4:80 “244 1-86 76 18-5 
19-0 4 0-443 . 19-2 
245 6 0-630 { 24-5 
24-0 . OTM 24-2 
22:3 ‘ 0-523 22-0 
18-5 . 0-965 5 23-1 
18-9 . 0-440 ! 19-0 
23-2 , 0-557 7 23-8 
24:5 ‘02 0-610 27 245 
24:5 . O-710 37 243 
21-9 a0 0-506 5 22:1 
22-3 OS 0-520 22-9 
245 { O-645 35 24-5 
25-6 “SO OS7T5 2 240 


(n/p * 104 
ft2 
(x) 
0-867 
0-453 
0-972 
O-445 
0-659 
0-720 
0-541 
0-793 
0-442 
O-S75 
0-635 
O-708 
O-516 
O-544 
0-686 
0-595 


a e- e- e- e e 


_— — 


~ 





Some heat transfer measurements were made heat transfer coetlicients for the frozen N,O, 
with air in the apparatus for comparison with system are included. These values were obtained 
N,O, The resulting experimental heat transfer using experimental values of U and t, from Tables 
coellicients, h, for air and N,O, are presented as 1 and 2 and calculating frozen properties there- 
functions of gas velocity in Fig. 5 where it may from. Values of Nwu were calculated from data 
be seen that values of h for N,O, are much higher in the literature, Ref. [21], using a value of 
and more scattered than for air. For comparison, 0-7 for Npr, since, for the frozen N,O, system, the 


Table 3. Values of dimensionless groups 





NNuy 


NRey NPry NNuy NPr' 3 NRey NPry 


1-630 o-vo2 2 47-8 2-680 0-006 
5-340 0-651 35-6 41-5 5-280 0-632 
1-310 O-T737 5 64 2510 O51 
10-500 0-625 be 66-4 10-400 0-627 
7-810 0-765 3-6 47-7 7 400 778 
7070 O-S34 4-2 w-9 7-770 O-821 
10-000 0-631 5 5s) 9-700 0-695 
5-610 O-9r0 j 47-2 6-820 O-865 
15-100 O-ozt iT! 7TO0 15-000 0-625 
13-400 0-704 2: 70-6 13-000 0-723 
12-700 0-746 i 69-2 12-200 0-766 
11-000 Os14 7-4 62-0 11-000 O-s10 
16-500 0-670 81-0 16-200 0-676 
18-800 0-682 7 87-5 1S-000 0-697 
15-300 O772 j 73-7 14-40 0-800 
17-300 0-720 3° S40 16-700 0-736 








Heat transfer in reacting systems 


Prandtl number varied between 0-71 and 0-73 
over the range of values of ft, in Table 2. 

Values of dimensionless groups are shown in 
Table 3. A log—log plot of Ny, /(Np,)'® vs. Na, 
was made employing both temperature and 
specific enthalpy as the driving force in accordance 
with equation (9a) and (9b) respectively. In 
Fig. 6 is shown the correlation with temperature 
as a driving force and with film properties based 
on a linear average of the free-stream and surface 
temperature. For comparison data from Ref. [21] 
for air and nitrogen with Nery taken as 0-70 are 
shown as a solid line. It is noteworthy that there 
is considerable scatter in this correlation parti- 
cularly in the tests (1 and 3) with the highest 
values of tg. An average deviation of 38 per cent 
from the line of Ref. [21] and a definite non-linear 
trend in the data are observed. Using specific 
enthalpy as a driving force and film properties 
based on a linear average specific enthalpy a 
This correlation 


similar plot is shown in Fig. 7. 
is seen to be much more consistent, with an aver- 
age deviation from the line for air and nitrogen 
of 5-8 per cent. A least-squares analysis of the 
data in the form of this latter correlation gives 


NNuy’ 


: (10) 
(N Pr)y 13 


0-559 (Np) 


It is noteworthy from Table 3 that where ft, 
is high, as in tests 1 and 3, there is a great variation 
in properties through the film due to the large 
difference between ty: and ty shown in Table 2, 
and thus there is a large difference between 
Nitey and Neer. Corresponding values of Nrr are 
less sensitive to a choice of film properties as 
suggested by Hansen [24]. 

The data were also plotted in a form suggested 
by Dovucias and Cuurcutis. [21] which is based 
on the concept of two parallel thermal resistances 
consisting of a laminar film resistance and a 
separated zone resistance. A better correlation 
results when specific enthalpy is employed as a 
driving force (with film state based on equation 
7a) than when the temperature is employed as 
driving force (with film state based on equation 
Tb). In the former case the data deviated on the 
average 5-3 per cent from the curve recommended 
in [21], equation (4) (Nery = 0-7). 


CONCLUSIONS 


The experimental results are in accord with 
the theory for the mechanism of heat transfer 
in a dissociating gas. Heat transfer rates are 
encountered which are several times greater than 
for non-reacting or frozen gases under similar 
conditions. The results indicate that heat transfer 
data from a heated cylinder in an equilibrium 
reacting system can be correlated in a simple 
manner using forms of correlation similar to those 
for 
correlation results when specific enthalpy is used 
as the driving force and the film properties are 


used non-reacting systems. A successful 


evaluated at the average specific enthalpy. Since 
the present study was restricted to low Reynolds 
numbers where the heat transfer through the 
laminar boundary layer accounts for the major 
contribution to heat transfer, further studies at 
higher Reynolds numbers up through the region 
of transition to a turbulent boundary layer should 
be made. Also a study at Reynolds numbers 
less than unity, where laminar flow around the 
cylinder without separation cecurs, would be 
of interest, where the assumption of a_ thin 
film used in the present investigation would be 
invalid, 

It would be of interest to study local behaviour 
by measuring local heat transfer coefficients 
around the cylinder circumference. Also, temp- 
eratures and velocities near the cylinder and in 
the wake should be measured and checked with 
theoretical predictions. Comparison of the width 
of temperature and velocity wakes should give 
relative magnitudes of the eddy conductivity and 
eddy viscosity and provide a test of the Reynolds 
Analogy (Ref. [22]), for reacting systems. 

A theoretical and experimental investigation 
of heat transfer in non-equilibrium reacting sys- 
tems should be made, and thus the role of reaction 
kinetics on the heai transfer processes determined. 
At the present time, there are no detailed data 
with regard to this problem. 
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at average film specific enthalpy. 


Heat transfer in reacting systems 


NoTATION 


area normal to heat flux ft? 
isobaric specific heat B.Th.U Ib F 
binary diffusion coefficient (ft? sec) 
outside diameter of cylinder (ft) 
specific enthalpy of equilibrium mixture in- 
cluding enthalpy of dissociation and relative to 
Hy,o, ~ © at 77°F (B.Th.U./Tb.) 
enthalpy of reaction per Ib. mol. of reactant 
equation (la or Ib) (B.Th.U./ Ib 
mole) 
heat transfer coefficient (temperature driving 
force) B.Th.U, ‘see ft® °F or B.Th.U. hr ft® “F 
heat transfer coefficient (enthalpy driving 
Ib, ft? sec 
B.Th.U, /ft 
length along cylinder ft 


force) 
thermal conductivity F sec 
molecular weight of associated species N,¢ Ys 
Ib Ib mole 
total pressure Dorce/ft® 
B.Th.U./see or B.Th.U. her 


radial distance ft 


heat flux 


outside radius of cylinder 
universal gas constant 
1-987 B.Th.U. Ib mole 
absolute temperature 
temperature 


free-stream velocity 


X or X, designation for |b mole of dissociated or associ- 
ated species respectively as defined by equation 
(1b) 
Greek letters 

pw = absolute viscosity Ib /ft sec 
p mass density Ib /ft® £)RT 
4 or @ function of 


€ degree of advancement, cf. Ref. [17], of reac- 


My,P/(1 


tion of type X, ——~ 2X for a system of mass, 
My.: in which case 1 > & >O 
Subscripts 
ce refers to circumferential point on cylinder 
d component of property attributable to dissocia- 
tion 
refers to film property based on arithmetic 
average temperature 
refers to film property based on arithmetic 
average specific enthalpy 
property at state of free-stream 
property at average state around cylinder 
surface 
thermal property frozen at constant degree of 
advancement of reaction 
Dimensionless parameters 
Nusselt 
Frozen 


NNuyg of f number 
N Leg 
N Pry 
NRey 


hDy ky or W’ Dy /(k, Cp)y 
Dy» pCpg ke 
(Cpp By orf’ 
(DyUp BP orf’ 


Lewis number 


Done . 
or g Prandtl number 


or Reynolds number 
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Mixing and contacting in gas-—solid fluidized beds 
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Shell Oil Company, Houston Research Laboratory, P.O. Box 2527, Houston 1, Texas 
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Abstract—The two-phase model of a gas-—solid fluidized bed is applied to steady-state and 
transient gas mixing tests. The two phases are : a dense phase consisting of particle aggregates and 
a dispersed phase of high vertical gas velocity. Gas and solids are exchanged continuously between 
the two phases by disintegration of particle aggregates and coalescence of dispersed phase solids. 

Expressions for gas back-mixing and for the standard deviation of the gas residence time 
distribution are derived and it is shown that the gas eddy diffusivity of the dense phase and the 
rate of mass-transfer between the two phases may be determined with the aid of suitable mixing 
tests. 


Zusammenfassung—Das Zweiphasenmodell eines Gas-Feststoff-Fliessbetts wird auf Versuche 
iiber stationiire und nicht-stationire Gasmischung angewandt. Die zwei Phasen sind: eine 
dichte Phase von Feststoffaggregaten und eine dispergierte Phas emit grosser Gasgeschwindigkeit. 
Zwischen beiden Phasen findet Austausch von Gas und Feststoff statt, zufolge Desintegration der 
Aggregaten und Bildung neuer Aggregaten von dispergierten Feststoffteilchen. 

Ausdriicke fiir Gasriickvermischung und fiir die Streuung in der Gasverweilzeit werden 
abgeleitet. Ferner wird gezeigt, dass sich die effektive Gaslingsdiffusion in der dichten Phase 
und der Gasstoffaustausch zwischen den beiden Phasen mit Hilfe geeigneter Mischversuche 
bestimmen lassen. 


Résumé—Le modéle en deux phases pour des lits fluidisés 4 l'aide de gaz est appliqué sur des 
essais pour déterminer le mélange de gaz dans ces lits. Les deux phases sont: une phase dense 
se composant d’agrégats de particules solides et une phase dispersée de grande vitesse de gaz. 
Le gaz et les solides sont échangés continuellement par suite de la désintégration des agrégats 
et de la coalescence des solides dispersés, 

L’auteur déduit des équations pour le mélange vertical de gaz et pour la dispersion du temps 
de séjour de gaz. Tl démontre que la diffusion verticale de gaz dans la phase dense et le transfert 
de matiére dans le gaz entre les deux phases peuvent étre déterminés a l'aide d’essais de mélange 


appropriés. 


INTRODUCTION 


Tux behaviour of gas—solid fluidized reactors has 
been the subject of many studies and discussions 
in the past decade. Although fluidization research 
has made considerable progress, the fluid mech- 
anical aspects are still not completely understood. 
This situation makes it difficult to explain in 
some detail what is going on in a fluidized bed 
when heterogeneous reactions take place. 

Earlier studies [1] on gas back-mixing in fluid- 
ized catalyst beds have shown that the down- 
mixing increases as the bed diameter is increased. 


*Present address: Koninklijke /Shell-Laboratorium, 
Maatschappij). 


Extrapolation to commercial scale leads to the 
conclusion that a large-diameter fluidized reactor 
may, as a first approximation, be considered as 
one well-mixed stage. 


Later information from tests on commercial 
units, as well as from experiments on a smaller 
scale, has indicated that the actual situation is 
much more complicated than was suggested by 
the earlier work. During sampling in the re- 
generator of a commercial catalytic cracking unit 
[1], it was observed that oxygen concentrations 
in the bed were appreciably lower than the 


Amsterdam, Netherlands. (Shell Internationale Research 
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oxygen concentration in the stack gas. In order 
to explain this phenomenon a two-phase fluidiza- 
tion mechanism was postulated, according to 
which the major portion of the gas passes through 
the bed in the form of * bubbles ” or * cavities.” 
In these discontinuities the oxygen concentration 
is considerably higher than the dense continuous 
phase, where most of the oxygen is rapidly con- 
sumed by burning. The low oxygen concentra- 
tions which were observed are then the result of 
preferential sampling from the dense phase. It 
was concluded that a fluidized reactor might be 
ineflicient in establishing gas-—solid contact. 

The 


developed 


theory of two-phase fluidization was 


more Toomey and 
(2) (3). 
Pansine [4] and Lewis et al. [5]. In the math- 
ematical model a transfer coetlicient is introduced 


explicitly by 


JOHNSTONE Suen and JoOuNsTone 


which describes the interchange of gas between 
bubble phase and dense phase. Gas back-mixing 
in the bubble phase is neglected. For the dense- 
phase gas only the two extremes of zero back- 
mixing or complete vertical mixing were consid- 
ered. The latter is evidently a serious restriction 
since it excludes the use of an eddy diffusivity as 
a separate parameter. May [6, 7] was the first 
to introduce a vertical diffusion coeflicient for 
the dense-phase gas in the two-phase model. In 
the evaluation of this model with the aid of gas 
and solids mixing tests it was assumed that the 
vertical solids diffusivity was equal to the dense- 
Although this 
assumption would seem to be reasonable it need 


phase vertical gas diffusivity. 


not necessarily be true. 

In the present paper the model of May is 
applied to gas mixing experiments to determine 
both dense-phase gas vertical diffusivity and gas 
interchange between the two phases without 
using the assumption of equal diffusivity of solids 
and dense-phase gas. 


THE PHYSICAL PICTURE 
The 


fluidization is the formation and disintegration of 


main feature of gas-solid aggregative 
aggregates of particles in a sort of dynamic 
equilibrium. Large aggregates fall down through 
the bed, carrying with them entrained gas until 
they are broken up into smaller fragments ; small 


VAN Deemwren 


aggregates or single particles are carried upward 
by the gas stream until they coalesce with other 
aggregates. The continuous disintegration and 


coalescence of aggregates imposes a random 
motion on the vertical gas flow. 

As the gas velocity is increased, the aggregates 
They 


The result 


are more easily disintegrated. have 
shorter lifetime and become smaller. 
is that the average bed density decreases. With 
further increase in velocity, finally all the aggre- 
gates and single particles are carried upward 
and dispersed-phase pneumatic transport becomes 
possible. 

When the gas velocity decreases, the aggregates 
become larger until they eventually form a con- 
tinuous phase through which the gas, together 
with a small fraction of the particles, flows up 
in “cavities.” With still smaller gas velocities 
the bed becomes quiescent and the “ cavities” 
disappear. 

In the 


boundaries, 


the wall and other solid 


formation of 


region of 
the 
with longer lifetime is promoted. 


larger aggregates 


ar 


Downtlow of 


gas and particles will, therefore, often take place 


the wall, although for large-diameter 


vessels it is just as possible that local downtlow 


along 


occurs in the centre, 

In general it will not be appropriate to speak 
of the dense phase as the continuous phase, 
because there is no sharply defined phase reversal. 
In the region of the wall, for instance, the dense 
phase might be continuous in character whereas, 
more centrally, dense aggregates are present in a 
continuous dilute phase. 

In order to define more clearly two phases in a 
fluidized bed, it is convenient to consider an 
instantaneous situation. Part of the gas will have 
high vertical velocities and will carry only very 
small aggregates and single particles. This will 
be called the dispersed-phase gas. The dense- 
phase gas consists then of the gas in the interstices 
of the large aggregates. 

In the dynamic the 
disintegration and formation of aggregates effects 


situation continuous 
an exchange of gas (and particles) between the 
dense and the dispersed phase. This mass-transfer 
mechanism is of importance for heterogeneous 
reactions because it brings unreacted gas into 
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contact with the bulk of the solids and it removes 
products from the dense phase back into the dis- 
persed phase. 

The density of the aggregates is usually assumed 
to be equal to the density of the fluidized bed at 
incipient fluidization. The volumes of dense and 
dispersed phase may then be easily calculated 
from the overall bed density at any higher gas 
velocity. It should be realized that the bed 
expansion, and therefore the volumes of the phases 
is dependent not only on gas velocity, but also 
on vessel diameter, gas and solid properties, etc. 

Now what will happen when a gas-mixing test 
is carried out or a heterogeneous reaction is taking 
place ? 

When a pulse of tracer gas is injected at the 
bottom of the bed, initially most of the tracer 
molecules will travel up in the dispersed phase 
at a rather high speed. Gradually part of them 
will be transferred to the dense phase, where 
down-mixing takes place, until they are removed 
(stripped) back into the dispersed phase. De- 
pending on the height of the bed, or more accur- 
ately, on the number of transfers between dense 
and dispersed phase and on the degree of random 
motion, the output signal will exhibit more or less 
spread in residence times. A considerable spread 
in residence times, for instance, may be the result 
of either considerable down-mixing in the dense 
phase combined with rapid transfer between the 
phases or little down-mixing combined with a slow 
rate of transfer or any other situation in between. 
In other words, a pulse test alone is insuflicient 
to evaluate the parameters of the model. 

In a steady-state gas-mixing test a tracer gas 
is injected at a steady rate higher up in the bed 
and the tracer 
determined by sampling at various locations in 


concentration distribution — is 
the bed. There will be no appreciable concentra- 
tion of tracer gas below the injection point when 
the dewn-mixing is slow and, in addition to that, 
the rate of transfer between dense and dispersed 
Otherwise the 
gradients in the bed may become small. 


phase is rapid. concentration 

For the case of a heterogeneous reaction, a 
fluidized bed may be considered as two parallel 
reactors with exchange of both gas and solids. 
Without any exchange the reactor corresponding 


14 


to the dense phase would operate at extremely 
small throughput and almost zero space velocity 
and, therefore, yield practically 100 per cent 
The reactor corresponding to the 
dispersed phase would be a high throughput, 
very-low-conversion reactor. The total conversion 
would, therefore, also be very low. As the rate of 
transfer between the two reactors increases, the 
total conversion will increase until it levels off 


conversion. 


when complete participation of all the solids is 
approached. It is evident that, for incomplete 
participation of the solids, possible secondary 
reactions are enhanced. 

True bypassing will only occur when the rate of 
transfer between the dense and the dispersed phase 
is very small and the dispersed phase contains 


very little solids. 


MATHEMATICAL FORMULATION 

Let the gas in the dispersed phase occupy a 
volume fraction, F, and let the average vertical 
velocity be VV. When f and v represent the same 
quantities for the gas in the dense phase, then the 

superticial gas velocity u, is equal to 
u=FV+fv (1) 
The rate of mass transfer between the phases 
is assumed to be proportional to the difference in 
the average concentrations C and ¢: rate of 

k(C —e), 
The random motions in the gas will be described 
with the aid of an eddy diffusivity, E. Although 
the vertical gas mixing has only partly to do 


transfer per unit height 


with the behaviour of the dense phase it will 
nevertheless be assumed in the following that the 
eddy motions take place in the dense phase only. 

The eddy diffusivity concept imposes certain 
restrictions on the applicability of the model. If, 
for instance, there should be a continuous fall of 
aggregates from top to bottom along the wall, 
the eddy currents would lose their random 
character and the eddy diffusivity would afford 
at the most a rather qualititive picture*. On the 


*Another possible model circumventing this difficulty 
is perhaps: countercurrent flow of gas in the phases com- 
bined with mass transfer. The rate of gas down-flow in the 
dense phase (caused by the falling large aggregates) would 
then be the important factor for the degree of down- 
mixing in the bed. 
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other hand it should be realized that the model 
should describe in the first the ideal 
fluidized bed. Deviations such as extreme wall 


instance 


effects or uneven gas distribution will then be 
characterized by unusual behaviour of the para- 
meters of the model. 

The material balances for a differential height. 
dx, in the fluidized bed lead to the following 
differential equations 


(2) 
» « eC 
iE 0 (3) 
A ioe 


Y and 


dC Mf are equal to zero. For the case of a chemical 


For a steady-state gas-mixing test d« 
reaction d¢ M and MM may be expressed in 
terms of rate constants and concentrations of the 
two phases. 

The 
(F = 0) and infinitely high rate of mass transfer 


limiting cases of incipient fluidization 


(k *) lead to the same simplification of (2) 
and (3). 


f E Ale 
‘ ° — 


UE oe 


STEADY-STATI rests 


GAS BACK-MINING 


In a gas back-mixing test a tracer gas is injected 
at a steady rate higher up in the bed and the 
concentration gradient below the pont of injection 


is observed. Fig. 1 shows the result of a typical 
= 
\ 


UNC ENTRATION 


IRAC ER « 


ATive 


e-5 - 


a 
a = ee 
a— 


BELOW INJECTION ABOVE +10 INCHES 


Fic. 1. Back-mixing in a 3 in. diameter fluidized bed 
according to Mason [8 


155 « glass beads 
© Centre 


Tracer injection int he centre. 
1-2 ft see superficial gas velocity 
e Wall 


small-scale back-mixing test. What happens 
apparently is that tracer is transferred from the 
centre to the wall from where it is carried down 
below the injection point until it is stripped from 
the dense phase by the uptlowing gas in the centre. 
This behaviour bears a close analogy to the theor- 
The 
profiles indicate that it takes about 1 ft. for the 


etical model. downstream concentration 
tracer to become transferred to the dense region 
at the wall: in other words. the height of a 
transfer unit is roughly equal to 1ft. It may, 
therefore, be expected that when the bed level 
above the injection point becomes smaller than 
1 ft.. the back-mixing decreases. 

It should also be expected that tracer injection 
at the wall will give more back-mixing than tracer 
injection in the centre since in the first case no 
additional transfer step is required, This has been 
confirmed by Mason [8, 9] who compared single- 
point injection at the centre with multiple-point 
injection over the cross-section. 

In large commercial vessels the analysis of a 
steady-state mixing test is complicated by radial 
mixing effects. The mathematical treatment will, 
therefore, be contined to small-scale experiments 
which, except for the transfer between dense and 
dispersed phase, are essentially one-dimensional. 

Under steady-state conditions, equations (2) 
and (3) reduce to: 


Subtracting (5) from (6) gives 


d?¢ de 
KE vf 
J dr? * & 


or, integrated, 


de 
SE vie | 


Ft 
dr 


constant (7) 
Applying equation (7) to a level in the bed far 
enough above the tracer injection point, 2 0, 
where de,dx may be considered as being substan- 
tially equal to zero and where C and ¢ will have 
the same value C,, it is seen that 
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» de 


E mC 
da ’ 


Si vfe V 


(vf VP), 


Below the injection point 
Ofora (0) 
the 


regions the concentration ¢ is continuous, but by 
(8) and (9) 


Across the boundary, ! 0, between two 


V F A c “ (, } kK A (< ACTOSS 7 0 (10) 
adi 


The fraction of tracer injected into the dense 

phase is equal to 

hy de 
4(z) 


= (11) 
ud dr) --0 


% 
0 
When all the tracer is injected into the dense 


phase : 


When all the tracer ts injected into the dispersed 


phase : 


a=0,V FAC =uC,, and A (““) 0 


di 


The boundary conditions at the bottom of the 


bed are 


de 
Cc Oand BE . 


ad 


vetora L, (12) 
The solution of the differential equations (8) 
and (9) subject to the boundary conditions (10), 
(11) (12) 
methods. 


and may be obtained by standard 
The concentration im the dense phase is found 


to be 
A exp | r) Bexp (‘ for. 0 (13) 


2£\,. 
Dexp (5 for. 0 (14) 
where A is the positive and y« is the negative root 
of the equation 


2, (AL Li 


kL? fv 
iin . —- 1+ - » (15 
+ (Ep EY ze (+ pr) -° 


The constants 4, B, and D are 
2(L/k)(u f) fa 
A fe 
, A (Li &) 
ft (Lv, BE) 
A (Lv/ EB) 
fe (Lv Rk) 


exp (ye A) 


exp (ys A) 1 (18) 
The concentration in the dispersed phase may 
be derived from (13) and (14) with the aid of 
equations (8) and (9), 
Under many conditions a considerable simpli- 
feasible. 


fication is When the superticial gas 


velocity is net too small (say 0-5 ft. sec.) most 
of the gass passes through the bed in the dispersed 
phase; Le. fv may be neglected compared with 
BV. 


the bottom of the bed is negligibly small. Then 


Finally, often the tracer concentration at 


L. is intinite for practical purposes, Under these 
the the 


conditions solutions become for dense 


phase 
l a) 


u“ 


and for the dispersed phase 


Cc B | l 3 
{« ) exp ( 
’ ». a . 
«/) - 


'e) 


(oo 


where 8? (23) 


 TEk 


Denoting 


Hi, height of a transfer unit (24) 


and 
SE 


Hy, 
, - 


height of an eddy diffusion unit (25) 


so that 
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a 


2 26 
ii, (26) 


p? 
it is seen that the concentration profiles are com- 
pletely determined by the dimensionless para- 
meters «, H,/H, and «/H,. 

Figs. 2 and 8 show as examples the 
Hi, He 6, 
persed phase) and H,/H, 
injection in the dense phase). 


cases 
0, (tracer injection in the dis- 
6, « 1 (tracer 


‘moms 
r 


SS 


DISPERSED 
PHASE 


DENSE PHASE 


DISPERSED PHASE 


RELATIVE TRACER CONCENTRATION 


Hy INJECTION 


Tracer injection in the dispersed phase, a = 0, 


H, = 6Hy. 


7 


ENTRATION 


DENSE , PHASE 


CONE 


DISPERSED 
PHASE 


RELATIVE TRACER 


H INJECTION 


8. Tracer injection in the dense phase, a = I, 
H, = 6H,p. 


The value of « appears to have a pronounced 
influence on the nature the concentration 
profiles. The close analogy between Figs. 1 and 
2 suggests that tracer injection in the centre is 
approximately equivalent to tracer injection in 
the dispersed phase. The high concentrations 
observed in the centre above the injection point 
might be the result of sampling from the plume 
emerging from the injection pipe. Below the 


of 


DrEemTER 


injection point the samples drawn in the centre 
might not be representative of the dispersed 
phase. These samples will then show too high a 
concentration. 
however, may be considered to be dense-phase 


Samples drawn along the wall, 


samples. 

When Figs. 1 and 2 
assumptions it would seem that H, ~ 10 in. and 
Hy, = 2 in. for the conditions of Fig. 1. It should 
be realized, however, that this estimate is rather 
inaccurate due to the scatter in the experimental 
results and particularly because of the uncertainty 
in the value of a, which was assumed to be equal 
to 0. For instance, the combination « — 0.02, 
H} = 20 in. and Hy = 1 in. would fit the experi- 
mental data just as well. It appears, therefore, 
0 would lead to a lower 


are compared under these 


that the assumption « 
limit for H, and an upper limit for H,. These 
limits may most conveniently be found by plotting 
the logarithm of the concentration along the wall 
as a function of the vertical position. Upstream, 
the experimental points should, according to 
equation (19), fall on a straight line with a slope 


p—1 


§ =F — 
2, 


The slope, 8, is usually termed the mixing slope. 
The intercept of this straight line with 2 = 0 is, 
for « = 0, equal to (8 — 1)/28. From slope and 
intercept the limiting values of H, and H, are 
easily calculated, 

In this way some of the older back-mixing 
experiments [1] where the tracer was injected in 
the centre have been evaluated. Since the results 
are rather qualitative details of this evaluation 
will be omitted. No systematic effect of the 
superticial gas velocity (0-1-1-1 ft/sec.) or particle 
diameter (cracking catalyst 55-110 «) was appar- 
ent within the limits of accuracy. The effect 
of bed diameter may be seen from Table 1 below : 


Table 1 





| 
| Hy (Maximum) (ft.) 1H), (Minimum) (ft.) 


Diameter (in.) 


| 


Os-12 
O-4-1°8 
20-255 


O35 
O-4+-1°2 
12-19 
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Although only (rather inaccurate) limiting 
values are obtained it appears from these results 
that in small-diameter fluidized reactors mass- 
transfer might be even more important than the 
effect of eddy diffusion. Bed heights of several 
feet are required to ensure that mass transfer is 
not a limiting factor. 

In many of Mason’s back-mixing experiments 
[8] the tracer was injected in both the dense and 
the dispersed phase. For such a case the value of 
a is too large to allow a similar calculation. It 
follows, however, from equations (19) and (21) 
that the mixing slope (27) is the same for the dense 
and the dispersed phase and, moreover, that it is 
independent of the value of « It will be shown 
in the following how H, and H, may be calculated 
time 


from the mixing slope and a_ residence 


distribution measurement. Fortunately Mason 
has carried out both types of mixing test for the 


same systems. 


RESIDENCE TIME DISTRIBUTION MEASUREMENTS 


For the evaluation of residence time distribution 
measurements an expression has been derived 
for the standard deviation, s, in terms of the 
parameters of the two-phase model (Appendix). 
To simplify the calculations the assumption has 
been made that there is no net vertical flow of 
gas through the dense phase. The expression 
reads 

att «<9 +—+— 
Nk ns | 
1)] 1} {exp 
1)] — exp | 


Ww 


{exp bnnlq 


exp [J m (q- 


hnq+1)—1} (28) 
jn, (q + 1)) 


where 


- relative standard deviation (29) 
(30) 
(31) 


number of transfer units* 


*n, is equal to the “cross flow ratio” introduced by 
May [6, 7). 


L ws ul 
H, fE 


= number of eddy diffusion units 


(32) 


fraction of gas present as dispersed phase (33) 


In Figs. 4, 5, 6 and 7 lines of constant standard 
deviation are drawn in ™% vs. my plots for four 
values of the parameter ¢. The positions of the 
n, values on the vertical axes were calculated by 


”» 


“ta — ¢)° zr" 


>)? HM, (34) 


Lu 
@M so 2 


Fic. 4. Standard deviation of the residence time distribu- 
tion and number of back-mixing units for two-phase 
fluidization ; ¢ = 0-2. 


In other words, the vertical co-ordinate is propor- 
tional to the height of a transfer unit. 
The ng positions on the horizontal axis are given 


Np) |] = - (35) 
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Fie. 5. Standard deviation of the residence time distribu- 


tion and number of back-mixing units for two-phase 
fluidization ; 6 = 0-3. 


The expression (35) results from (28) when ny 


becomes infinite. ” is equal to the number of 
completely mixed stages giving the same standard 


deviation (see Ref. 10). The horizontal co-ordinate 


Table 2. 





Fic. 6. Standard deviation of the residence time distribu- 


tion and number of back-mixing units for two-phase 
fluidization ; dé = 0-4. 


is, therefore, proportional to the height 7, of a 
mixing stage in the dense phase. 

Under the assumption of no net gas flow in the 
dense phase the expression for the mixing slope 


Experiments by Masowx in a 3 in. diameter, 37 in. high bed 





Particle size Superficial velocity 
Particles (nu) (ft sec) 
Cracking catalyst | 70 O-4 
70 0-6 
70 oo 


Cracking catalyst O4 
O-6 
0-9 





Glass beads O-4 
oO-o 
ov 


Glass beads 0-6 


Relative standard Reciprocal mixing, Hyp 


deviation slope (ft) (ft) 


O87 
OST 
OST 


Ow 
| O36 
0-30 


O-S4 
O-s84 
O-S4 


O40 
0-33 
0-23 


0-62 
0-62 
0-62 


0-06 
0-16 
0-27 


0-52 0-06 
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Fic. 7. Standard deviation of the residence time distribu- 
tion and number of back-mixing units for two-phase 
fluidization ; 6 = 0-5. 


obtained in a steady-state back-mixing experiment 


was found to be (ef. equation 27 


4 p/n) 
2H, 


s_VvU 


(36) 


It was pointed out earlier that the theoretical 
mixing slope (36) is not dependent on the relative 
amounts of tracer gas injected in each of the two 
phases. 

When both o and S are measured, n, and np, 


or H, and Hy, may be determined if the value of 


¢ is known. This may be most easily achieved 
graphically. To this end, lines are drawn in Figs. 
4, 5, 6 and 7 on which the product of bed height 
and mixing slope 


np = LS = bn, [a $ np/n,) 1] (37) 


This product may be called the 
number of back-mixing units. A value 


indicates little down-mixing ; 


-_? 


is constant. 
much 
greater than 1 a 
value smaller than 1 means considerable down- 


mixing. The value of ¢ may be estimated from 
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the measured bed expansion or from bed expan- 
With the aid of the grid, 
constant, the values of np 


sion correlations, 
ao = constant, ny 
and n, can then be read directly from the corres- 
ponding graph. In this way a few of Mason’s 
experiments could be analysed. The results are 
given in Table 2. 

Comparison of Tables 1 and 2 shows that the 
limiting values of Hy, and H, are reasonably 
close to the values calculated from the combined 
information of steady-state and transient mixing 
tests. Although fluidized beds of glass beads 
appear to be quieter than cracking catalyst beds 
the ratio H,/H, is roughly the same. 

May [6, 7] has carried out solid mixing tests 
with fluid cracking catalyst in vessels of 3 in., 
15in. and 5ft. diameter. The solid diffusivities 
were at 0-8 ft/sec, superficial gas velocity 0-6, 1-6 
and 5 ft®/see respectively. The present study 
indicates that at 0-8 ft/sec superficial gas velocity 
u H,/f) is of the same 
order of magnitude for comparable vessel diameter. 


the gas diffusivity (FE 


From residence distribution measurements in 
a 3in. diameter unit for 12 and 35 ft bed height 
May finds a height of transfer unit of about 3 ft. 
This value is higher than the values derived from 
Mason’s data, but May points out that the 
experimental conditions are not directly com- 
parable. 

Coming now to the question of the physical 
significance of the present results, some specula- 
tion is perhaps justified. To this end, let us go 
back to the definitions. H, is a measure of the 
distance over which the bigger particle aggregates 
fall before they are broken up into smaller pieces. 
the 
dispersed-phase gas between a place where it 


H, is the average distance travelled by 


escapes from disintegrating aggregates to a place 
where it is trapped by aggregates which are 
being formed. The ratio H,/H, will, therefore, 
be roughly equal to the ratio of the linear velocity 
of the dispersed-phase gas and the downward 
velocity of the bigger aggregates. This conclusion 
is in agreement with the picture we have of a 
small-diameter fluidized bed, viz. an_ internal 
circulation pattern consisting of fast-moving gas 
pockets in the centre and a slow downward 
movement of dense catalyst masses along the wall. 
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Table 1 suggests that for increasing diameter 


the ratio H,/H, decreases although the absolute 


values of H, and H, increase. It would seem, 
therefore, that in small-diameter vessels the size 
of the aggregates is determined by the diameter. 
This results in lower values for 1,» and H,, but 
also in a higher value for H,/H, because of a 
the the 


aggregates. It may be speculated that for vessels 


reduction in downward velocity of 
of very large diameter the size of the aggregates 
is determined by the tluid-dynamical parameters 
of the fluidized system rather than by wall effects. 
This would mean that H, and H, approach a 
maximum value if the diameter is increased under 
otherwise identical conditions. The ratio H,/H,. 
however, will approach a miminum value deter- 
mined by ftluid-mechanical conditions. 


HrTEROGENEOUS REACTIONS 
The two-phase model can directly be applied 
to the Huidized reactors. 
May [6, 7] has carried out this calculation for 


calculation of bed 
a first-order heterogeneous reaction. On the other 
hand it is possible to evaluate the parameters of 
the two-phase model from the performance of 
the fluid bed as a reactor if at the same time a 
suitable gas-mixing test is carried out. In Fig. 8 
the relationship between conversion, nm, and ng 
is plotted for a typical case of a first-order 
heterogeneous reaction.* 

In the same Fig. lines of constant relative 
standard taken from 
Fig. 5. The two families of lines form a grid 
which determines np and n, when conversion and 


deviation are drawn, 


standard deviation of residence times are known. 
This procedure may in principle be applied to 
regenerators of commercial catalytic cracking 
units since the coke burning is first order in the 
oxygen concentration and the chemical kinetics 
are fairly well understood. The main difficulty is 
the accurate measurement of the spread in resi- 
dence times. The problems involved with transient 
mixing tests in large units have been described 
in the papers of Danckwerrts ef al. [11] and 


*The number of reaction units KL/u = 10, where K is 
the appropriate chemical rate constant. Assumptions : 
No net vertical gas flow in the dense phase. No catalyst 
present in the dispersed phase. 














Fic. 8. First-order heterogeneous reaction in a two-phase 


fluidized bed, 6 = 0-3. 


Hanpb os et al. [12]. Some knowledge on the value 
of the main parameters, H, and H,, of the 
model for large-diameter vessels is, however, 
absolutely necessary if one ever wants to design 
large fluidized bed reactors on a more scientific 
basis. 
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APPENDIX 

For the calculation of the standard deviation of the 
residence time distribution use will be made of the method 
indicated by vAN per Laan [10]. 

It will be assumed that the dense-phase gas velocity, v, 
is so small that the term vf (d¢/)z) in equation (3) can 
be neglected. The boundary conditions corresponding to 
a pulse injection of tracer gas are then 
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where 5 (é) is the impulse function, 
When the concentration C, is normalized in such a way 
@ 


that 8 (t) dt 1, then the concentration response at 


. 


the bed exit, a = L, is exactly the residence time distribu- 
tion function, R(t). Introduction of 
P 

S+Fr’ 

and application of the Laplace transformations 
- 
exp (— pr)Cdr, €= | exp(— prjedr 
0 0 

reduces the equations (2), (3) and (A.2) to 


d 
dpU + b my (C 
dg 


(1 — ¢) pé — 


dé 
The solution for (' is obtained by standard methods 


! a, exp (Af) 4 a» exp (AQg) + a exp (Agé) (A.5) 
a + @ + ay wii 


where A,, A, and A, are the roots of the equation 


ne [1 — dp + mas 
¢)p? ng — pry npg = 0 


4 (nm, + dp) #* - 
é( (A.6) 

and 
exp (A3) — exp (Ag) 


nm, + OP + Ay 


as 


(A.7) 
A, exp (A cp (A 
1 XP (Ay) 25 3) (A.8) 
Ag My ? t Ag 


- A exp (A,) - exp (A,) 
Is : \ 


Ay ™ + OP ts 


transformation 


(A.9) 


The Laplace 
distribution is 


of the residence time 


O(¢ =1) = Rip) = 
a exp (Ay) + % exp (Ag) + 43 exP (Ag) 


(A.10) 
% + a + ay 


The moments of the residence-time distribution function 
are found by forming the first, second, etc. differential 
quotients of R (p) with respect to p and letting p go to 
zero. The calculation for the present case requires dif- 
ferentiation of the expressions (A.6) up to (A.10) inclusive. 

In the following differentiation with respect to p will 
be denoted by a prime and limiting values for p = 0 by 
a subscript 0. 

It follows from (A.6) up to (A.10) inclusive 


0, Ag, A 


2 430 = — "Ee Aso + Aso 


Ny X49 = EXP (Aga) — EXP (Ago), X29 = 9, 
go 0, Ry =1 
(A.10) will be written in the form 
RZ a, 2 a, exp (Aj) 
Differentiating (A.12) gives 


t’ Za, 4 


‘ R22! = 


2 a,’ exp (Aj) 4 2 a; Aj’ exp (Ay) 


For p = 0 (A,13) becomes, considering (A.11) 


Oy He’ + ay + ag’ = 


ay’ EXP (Ago) + a3 EXP (Ago) + %4 A,’ (A.14) 
By differentiating (A.6), (A.8) and (A.9) and using (A.11) 


it is found that 


he = —1, 


EXP (Ago) —1 
NE Ny ' 


1 — exp (Ago) 


Log X39 


ne, (A.15) 
and from (A.14) and (A.15) 
R,’ = -1 (A.16) 


(the average dimensionless residence time is equal to 1). 
Differentiating (A.13) yields 
Rv Sa, + 2h Lal + R La,” = La," exp(d) + 
2 Z a,'d,;’ exp (A;) + 2 a;,A;" exp (Aj) + 
Za; (A, exp (A;) 


For p = 0 (A.17) becomes, using (A.11), (A.15) and (A.16) 


(A.17) 


ay Rg” = 2 aq [1 + Ago’ EXP (Ago)] + 
+ 2 a9 [1 + Ago’ EXP (Ago)] + 
b+ gq” [exp (Ago) — 1] + ago” [exp (Ago) — 1} 4. 


+ a : re 4 + 24 ~ 6] » (A.18) 
‘ Ne ny 


Expressions for 29° and a9’ may be derived by differen- 
tiating (A.8) and (A.9) twice. It is found that 
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Mg Ny Aq = 2 tag’ Ang + 2a, fexp[m(q—1)/2] —1} fexp[ —m(q + 1)/2] —1} 


” 
"Ep my a0 - 30 30 


ig [EXP (Agg) — 1] — * Leo 


io [1 — exp (Agg)] + * Ago XP (Aga) exp [nm (q — 1) 2] — exp[— mq + 1) 2] 


\gq Ang t+ 2 %gq' Ago (¢ The variance of the residence time distribution function 


exp ( )] (A.20) is found from 


20") 


2 2 2q 


When the expressions for ayo, % 9°, %o @Nd ag, are 4! "a 4)? 


introduced into (.\.18) one obtains Np 


£773 [EXP (Azo) 1] [exp (Ago) 


NE k 


2 
” gE 


fexp [n,(q ’ I} {exp[ —nm(qe1)2 1} 


a A.26 
Db ' exp [ m, (q ‘ exp | ny (q 1) 2] , 


\ few numerical values are found in the following table. 


Values of ov «2 my) if (F 4 fp 





(30 \gq) (8 \n0 


It follows from (A.6) that 


\ 
1 


Mg + Ay 


OSS WOoSs) O1650 OSISL OLS) 1 


or, by differentiating and letting p go to zero 2 OFO852 O-O822 01608 06-2020 O-5758 O-9a65 


The use of (A.11) and (1.23) finally reduces ((.21) to 


Pp 


, loo \go 4 (10330 OS810 | 01563 2868  O5274 | O8522 

OOo OTOL OF 14904 082646 O 4614 7358 

O“oszl W0OT55 O1STZ2 OVE | O-3689 | O-S6TT 

OOo00 0855 OF 1108 1687 OO 2293 O-Bzs80 
0 0 0 0 0 rT) 





NOTATION 


concentration moles ft' Ny = KL« number of reaction units 

eddy diffusivity ft? sex standard deviation 

volume fraction ‘ mixing slope 

height equivalent to a completely mixed time 

stage ft superticial velocity ft, sec 
height of an eddy diffusion unit ft v. linear velocity ft sec 
height of a transfer unit ft vertical co-ordinate ft 


mass-transfer coeflicient sec™! . 
fraction of tracer injected into the dense 


first-order reaction rate constant seco! 
phase 


height ft 


average residence time 
number of completely mixed stages . 


number of back-mixing units relative standard deviation 


number of eddy diffusion units fraction of gas in the dispersed phase 
number of transfer units ‘ 
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Heat and mass transfer in the presence of an inert gas 
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Abstract—A new graphical method for solving enthalpy transfer problems is described in which 
a modified enthalpy is the transfer potential. The condition line is determined by successive 
approximations. An entire condition line is constructed for each approximation, and the resulting 
approximate line determines an auxiliary line used for the next construction. The method was 
used to calculate enthalpy transfer coeflicients from experimental data on enthalpy transfer in a 
urid-packed tower for the tetrachloroecthylene air system, which are reported. Three types of 
countercurrent operations were studied : 

(1) ‘The vapour content of the gas is increased and the gas is heated by hot iiquid. 

(2) The vapour content of the gas is increased and both the gas and the liquid are cooled. 

(3) The liquid is heated by the bot gas and the vapour content of the gas is either increased 
or decreased. 

The enthalpy transfer coeflicients for operations of the first type were lower than those 
for the other two. The existence of fog during dehumidification had no apparent effect on the 
coetlicients. 

The application of the graphical method to the design of a cooler-condenser is described. 


Résumé Les auteurs donnent une nouvelle méthode graphique de résolution des problémes de 
transfert d’enthalpie dans laquelle le potentiel de transfert est une entropie modifiée. L*équation 
de la droite opératoire est déterminée par des approximations successives. Une droite opératoire 
est construite pour chaque approximation et la droite approximative qui en résulte détermine 
une droite auxiliaire utilisée pour la construction suivante. La méthode est utilisée pour calculer 
les coeflicients de transfert d’enthalpie a partir de données expérimentales sur le transfert d’en- 
thalpie dans une tour a garnissage a grilles pour le systéme tetrachloroethylene-air. Trois types 
dopérations a contre courant ont été étudiées : 

1. La teneur en vapeur du gaz est augmentée et le gaz est chauffé par le liquide chaud. 

2. La teneur en vapeur du gaz est augmentée et le gaz et le liquide sont tous deux refroidis, 

3%. Le liquide est chauffé par le gaz chaud et le teneur en vapeur du gaz est soit augmentée, 
soit diminuce, 

Les coeflicients de transfert d’enthalpie, pour les opérations du premier type sont inférieurs 
& ceux des deux autres, 

L/ existence de brouillard pendant la deshumidification n'a pas d’effet apparent sur les co- 
eflicients. 

Les auteurs décrivent l'application de la méthode graphique au calcul dun condenseur. 


Zusammenfassung—-Es wird eine neue graphische Methode fiir die Lésung von Enthalpie- 
iibergangsproblemen beschrieben, bei der cine moditizierte Enthalpie das Ubergangspotential 
bildet. 

Die Zustandskurve wird durch schrittweise Anniherungen bestimmt. Die gesamte Zustands- 
kurve wird bei jeder Anniherung konstruiert und die hieraus sich ergebende Naherungslinie 
ergibt cine Hilfslinie, die fiir die Konstruktion der nichstfolgenden Anniherung benutzt wird. 
Die Methode wurde zur Berechnung von Enthalpie-Cbergangskoeflizienten aus experimentellen 
Ergebnissen des Enthalpietibergangs in einer Fillkérperkolonne fiir das System Tetrachlorithylen- 
Luft benutzt, die mitgeteilt werden. Drei Typen von Gegenstromverfahren wurden untersucht : 

1, Der Dampfgehalt des Gases wurde erhOht, und das Gas wurde durch die heisse Fliissigkeit 
erwarmt, 

2. Der Dampfgechalt des Gases wurde erhéht, und sowohl das Gas als auch die Fliissigkeit 
wurden gekihlt. 

3. Die Flissigkeit wurde durch das heisse Gas erwiirmt, und der Dampfeghalt des Gases 
wurde entweder erhdlt oder gesenkt. 

Die Enthalpie-UCbergangskoeflizienten waren bei Versuchen vom Typ (1) niedriger als bei den 
beiden anderen. 

Das Auftreten von Nebeln wiihrend der Trocknung hatte keinen augenfilligen Einfluss 
auf die Koeflizienten. 

Die Anwendung der graphischen Methode bei der Auslegung eines Kiihler-Kondensators 
wird beschrieben, 


*Present address: Consolidation Coal Company, Library, Pennsylvania. 





Fr. u 


INTRODUCTION 


Wuen transferred simul- 


taneously 


heat and mass are 
interface it is 
enthalpy is the 


For mixtures 


across a gas-liquid 
convenient to that 
quantity that is transferred [1] 


consider 


of water vapour and air, for which the Lewis 
relationship, hg kg = C,. is valid, the enthalpy 
difference is the potential that establishes the 
rate of transfer. For the general case 
he «aii (=<) 
ky "\Pr 


Pas ta 

(1) 
: I hh 
and the simple enthalpy potential cannot be used 
but a modified enthalpy term that includes the 
Schmidt and Prandt! 


proposed by Mizusmina and Koroo [2] and 


groups is required, as 
Mizusuina and Nakasma [3]. 
é is a modilied enthalpy defined by 


¢=C,¢(t —t) +A H 2) 


Lewis and Wurrr [4] also use a modified enthalpy 
potential for enthalpy transfer, but they arrange 
the groups in another way. 

In this paper a graphical method of solving 
enthalpy transfer problems is described in which 
the modified enthalpy is included. The method 
was used to calculated transfer coetlicients from 
experimental data on the tetrachloroethylene—air 
system, which are also reported here. In the 
graphical method the condition line for the 
enthalpy of the gas at any temperature is assumed 
to start the Method of 
Approximations is used convergence is 
obtained for the assumed and calculated condition 


and the 
until 


calculations 


lines. 
Tueroreticat Discussion 
In the air-water system the slope of the 
condition line for any enthalpy transfer process 
is [1] 
di 
G « 
(3) 
dt. 


For the general case 


dig 


4 
dt (4) 


ba = 
G ta ¥ 


The derivation of equation (4) is similar to that 


15¢ 
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given by Mickiey for the air-water system [1], 
except that equation (1) is used in the derivation 
instead of the Lewis relationship. Equation (1) 
is derived by equating Co_nurn’s factors for 
heat and mass transfer and converting the mass 
transfer coefficient to the use of a potential 
expressed as pounds of vapour per pound of gas. 

To illustrate the graphical method, the coefli- 
cients, h;, kg and he, will be calculated from the 
data on an humidification tower. 
The initial data include the rates of flow of the 
gas and liquid, the temperature of the liquid at 
the ends of the tower, the temperature and humi- 
dity of the gas at the ends of the tower and the 
tower height. The terminal values of ig and ¢¢ 
are calculated from the observed temperatures 
and humidities. 

The graphical calculation is shown in Fig. 1. 
The equilibrium and operating lines are re- 
presented by curves of enthalpy, or ¢, vs. tem- 
perature. Only the 4 equilibrium line is needed. 


adiabatic 


@.OR | BTU PER LB INERT Gas 





TEMPERATURE OF GAS OR LIQUID, “F 
Fie. 1. 


The equation for the enthalpy operating line is 


~%J 


to, — ‘se, Le 
This is an enthalpy balance line which connects 
the points for the enthalpy of the gas and the 
temperature of the liquid at the ends of the 
tower. The slope LC,/G is assumed to be con- 
stant so that the line is straight, since small 
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changes in the liquid rate caused by evaporation 
or condensation and changes in the heat capacity, 
C,, are neglected. 

The ¢ operating line on Fig. 1 connects the 
values of dg and t, at every point in the tower. 
This line is not straight and its construction is 
complicated because its shape is determined 
partly by the condition line, which is unknown 
at first. Only the terminal points on these lines 
can be calculated from the data, but the curves 
are uniquely related at all points. Therefore, 
if a condition line is assumed and plotted on Fig. 2, 
a single ¢ operating line can be determined and 
Rach this 
assumed condition line with a line of constant 


plotted on Fig. 1. intersection of 


= 
| 


GAS ENTHALPY 


¢, on Fig. 2 determines a value for the enthalpy, 
ig, which corresponds to a value of t, on the 
enthalpy operating line. These values of ¢, 
and ¢t, determine a point on the ¢ operating line 
to be constructed on Fig. 1. 

A condition line is now constructed to check 
the one that was assumed. A value for the slope 
of the tie line, — hy ay/kg ay,, is assumed and a 
series of lines with this slope is drawn on Fig. 1. 
Vertical lines are drawn through points E and B 
at which the tie lines intersect the equilibrium 
and the ¢ operating line. Horizontal lines are 
then drawn through points 4, which are located 
on the enthalpy operating line, by the vertical 
lines through points B. The vertical distance 


between points E and B represents the transfer 
potential, 4; — dg, at any point in the tower. 
The distances, (1/%) (¢; — ¢g), are laid off on the 
vertical lines through points £, upward from their 
intersections with the horizontal lines through 
points A. The points D located in this way 
give the so-called 1/y — auxiliary line. 

According to equation (4) the line C, D, 
has the slope of the condition line at the gas 
inlet. The initial slope is assumed to be constant 
for a short segment of the condition line to C, 
in Fig. 1. By trial and error the condition line 
is constructed in increments until it reaches C,,. 
When the correct initial value of hy ay/kg ayy is 
assumed the condition line passes through point 





C, and the line must also coincide with that 


n 


assumed to start the construction. When con- 
vergence has been obtained, kg ay, is found by 


graphical integration of 


Gf le 
kg ay = Z | 6, ae oe (6) 
h, ay is found from the slope of the tie line and 
hg @y is found from equation (1), using the mean 
value of C,. 

The humidification process shown in Fig. 1 
is one of the three types of operations that may 
occur in an enthalpy transfer device. These may 
be classified as follows : 
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1. The gas is heated and the liquid is cooled 


Enthalpy is transferred from the liquid to the 
gas and the vapour content of the gas increases, 
For this case the condition line slopes upward 
and to the right. This is similar to the operation 
of a humidilier or a water cooling tower. Fig. 5 


is an example. 


The 


Enthalpy is transferred from the liquid to the 


gas and the liquid are both cooled 


gas and the vapour content of the gas increases, 
The condition line slopes upward and to the 
left. Fig. 7 is an example. 


3. The gas is cooled and the liquid is heated 


Enthalpy is transferred from the gas to the 
liquid, and the vapour content of the gas may 
The latter case 


the operation of a cooler-condenser. 


increase or decrease. represents 


The condi- 


tion line slopes downward and to the left. Fig. 6 
is an example. 

The slopes mentioned are the over-all changes 
in the condition line; the actual slope may 
change considerably along its path. For example, 
the the 


line may slope upward to the left initially, cross 


in operations of first type condition 
the operating line and then slope to the right as 


the difference between the temperature of the 


ADSORPTION 


a 


IL 


z 
7 = 


TCE TANK = 


ws 


Fie. 3. 


and W, 


B. RerauiicKk 

gas and of the interface changes signs. In opera- 
tions of the third type the condition line crosses 
a line of constant humidity twice when humidi- 


lication of the gas is followed by delhumiditication. 


EXPERIMENTAI 


Individual tilm coetlicients were determined 
for the tetrachloroethylene—air system with the 
The tower was 1 ft. 
ft) 
of No, 6897 drip-point ceramic tiles obtained from 
Each tile 
consisted of sixteen parallel rectangular channels 


The 


non-stagyvered 


apparatus shown in Fig, 3. 


» 75 


square and was packed with six layers ( 
the General Refractories Company. 


separated by corrugated ribs. tiles were 
stacked 
The 


cent 


above each other in a 


insulated with Zin. of 


blocks. 


were placed above and below the packing so that 


way. tower was 


85 per magnesia Glass windows 
the presence of fog in the gas could be observed, 
The liquid was distributed over the packing 
by eighteen nozzles made from jin. pipe plugs 
drilled } in. 


liquid in low velocity streams onto splash bars. 


with holes. These projected the 
The liquid rate was measured with a calibrated 
rotameter. 

\ir at room temperature was blown into the 
bottom of the tower for the runs in which the 


air was heated by hot liquid. For runs in which 


VENTUR! METER 


| STEAM-AIR 
HEATER 


AiR 


ORIFICES 


ADSORPTION 


Gas 
TUBES 
, 
METER 


~ TE TRACHLOROETHYLENE 
~ THEAaMOCOUPLE 


\pparatus used for enthalpy transfer studies. 
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the air was cooled the inlet air was heated with 
When the 
air was to be humidified, liquid was metered 


a steam-air heat exchanger. inlet 
through a calibrated rotameter and sprayed into 
the exchanger where it was evaporated by the 
hot tubes. In 
(No. 34) 


burner was added to the air stream to raise the 


one run at high temperature 


combustion gas from a natural gas 
temperature. The rate of air flow to the exchanger 
was measured with a venturi meter and the air 
and gas to the burner were metered with orifices. 
The venturi meter and orifices were calibrated 
against a calibrated nozzle from the Taylor 
Instrument Company. 

The concentrations of vapour in the inlet and 
exit gas streams were determined by drawing gas 


samples through two U-tubes in series containing 


activated charcoal, which were weighed to an 
OOLe 


5 


The 
10 g. 


gas samples, which was measured with calibrated 


total amount of 


The volume of the 


accuracy of 
vapour adsorbed was 
wet-test meters, was about 1 ft®. The increase in 
weight of a third U-tube in the adsorption train 
was negligible, showing that adsorption of the 
vapour in two U-tubes was complete. The sample 
of the gas at the tower exit was drawn through 
a perforated stainless steel tube which extended 
across the tower and was heated electrically to 
The 


sample at the inlet was drawn through a_per- 


prevent condensation of the vapour, gas 
forated glass tube which extended across the inlet 
gas duct. Temperatures were measured with 
No. 26 


and a Leeds and Northrup portable precision 


gauge chromel-alumel thermocouples 


potentiometer. The thermocouples used in the 
four terminal streams were calibrated in place. 
The liquid leaving the packing was collected 
into one stream which flowed into a_ trough 
below the gas entry point for the temperature 
measurement, in order to avoid the effects of a 
the 


entering gas stream. The temperature measured 


horizontal temperature gradient and of 
in this way was influenced by the tower end- 
effects, but included in the 
balance of the The thermocouple for 
measuring the temperature of the 
liquid was placed in the feed pipe near the entrance. 

The thermocouple for measuring the tempera- 


these were heat 
tower. 


entering 


ture of the gas leaving the packing was placed 
about 8 in, above the distributor; the junction 
was enclosed in a stainless steel tube and the gas 
was drawn past it at high velocity. The tempera- 
ture of the entering gas was measured with a 
The 


duct was insulated and heated electrically with 


bare thermocouple in the gas entry duct. 


resistance wires to a temperature close to the 
measured temperature of the gas in order to 
reduce radiation errors. 

During a run a_ steady temperature was 
necessary in order to obtain a good heat balance 
because of the large heat capacity of the tower 
and the relatively small heat exchange across the 
liquid—gas interface. Temperature and flow rates 
were recorded at least three times during each 
run. Samples of the inlet and exit gas were drawn 
for analysis throughout each run. The duration 
of the runs was 15-20 min. 

The tetrachloroethylene was obtained from the 
DuPont Company. It boiled between 120-5° and 
122-0 C, 
Ref. [5]. 


Other physical properties are listed in 


Resuvrs 
The results of all experiments with satisfac- 


tory heat balances are shown in Fig. 4 and Table 1. 
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Fic. 4. Enthalpy transfer coeflicients for the tetrach- 
loroethylene- air system. 
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Table 1. 





kKe@m hy ay 
Run | Type of for exponent of for exponent of 
no. operation 4 4 4 


| 


Intinite 
2510 


162 


405 


DAS 


2015 Infinite 





In the later runs (No, 31 and higher) the technique 
5 per cent 
in the heat balance were enough to reject the data. 
Coeflicients were 


was improved so that deviations of 


calculated for some runs in 
which the exit gas was saturated and fog was 
present in the gas. Calculation of the coeflicients 
in these cases was possible when the slope of the 
tie-lines could be found, which was not always 
possible when the exit gas was saturated. 

As shown in Fig. 4, all of the coetlicients fall 
on one of two lines. The coefficients for the runs 
in which the gas was heated and the liquid 
was cooled are lower than those for all of the 
other runs. When the liquid was heated and the 
gas was cooled and its vapour content was 
higher at the exit than at the inlet condensation 
could occur in the upper part of the tower. 


Condensation could occur throughout the tower 
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Fic. 5. Construction lines for Run 28, type 1, air 


heated and humidified tetrachloroethylene—air system, 


when the vapour content of the inlet gas was high. 

Figs. 5, 6 and 7 show the graphical construction 
for calculating the gas- and liquid-film coetticients, 
kg @y, and hy ay, for various operating conditions 
with the tetrachloroethylene-air system. Figs. 5 
the construction of an element of 
The former shows the con- 


and 7 show 
the condition line. 
dition when the liquid leaves the tower below 


LIBRIUM LINE 


_— call 
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— 
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Fic. 6. Condition line for Run 40, type 3, showing alternate hum- 
idification and dehumidification tetrachloroethylene—air system. 
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C; Cry = CONDITION LINE 
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@,OR ENTHALPY, BTU PER LB OF INERT CAS 


A, An = ENTHALPY OPERATING LINE Fic. 7. Construction lines for Run 49, type 2, 


air cooled and humidified tetrachloroethylene—air 


-..==<2 CONSTRUCTION LINE system. 


¢-operating line 


1 , 
— — auxillary line 


enthalpy operating line 
condition line 

tie line 
= construction line 








the temperature of the entering air. The air is 
cooled at the bottom of the tower but it is re- 
heated rapidly. The minimum temperature 
occurs where the temperatures of the gas and of 
the interface are equal. 

The runs in which both the gas and the liquid 
were cooled were made to find why the coefficients 
fall on two distinct lines. Under these conditions 
the gas was cooled but its vapour content in- 
creased so much that the enthalpy of the gas 
increased and the liquid was cooled. Heat and 
mass were transferred in opposite directions 
across the interface. Since the data for these 
conditions fall along the upper line, the direction 
of sensible heat transfer across the interface 
affects the value of the coetlicient more than the 
direction of mass transfer. 

For the conditions shown in Fig. 6 the gas was 
first humidified in the tower and then dehumidified 
so that the condition line crosses a line of constant 
humidity twice. It is possible to have humidifica- 
tion of the gas near the bottom and then dehumidi- 
fication even when the inlet gas contains no vapour. 
For Runs 39 and 41 the condition lines show that 
the vapour content of the gas increased through- 
out the tower. There should be no fog formation 
in this case and none was observed. In Run 37 
the condition line for the top of the tower is 
parallel to the constant humidity lines and, as 
expected, no fog was observed. 

In a packed countercurrent tower fog can form 
only when vapour is transferred from the gas to 
the liquid. The condition line indicates that 


dehumidification takes place when its slope is 
greater than the slope of the lines of constant 
humidity. This condition alone is not sufficient 
to cause fog to form but the tendency to form fog 
increases with the slope of the condition line, 
that is, with the value of —dH,/dtg. The 
effect of fog formation on the values of the 
anthalpy coetlicients was not studied thoroughly 
because fog was never formed over the entire 
packing height. However, when fog was formed 
it had no apparent effect on the values of the 
coeflicients. 

For the conditions in Fig. 7 the gas was cooled 
and its vapour content increased with a net gain 
in enthalpy. However, about one-third of the 
total enthalpy transfer occurred after the gas 
had been cooled below the temperature of the 
liquid. In this case the transfer coefficient fell 
nearer those on the lower curve on Fig. 4. The 
extreme curvature of the ¢-operating line causes 
a corresponding curvature in the condition line 
near the end. The condition line is curved because 
the gas was cooled rapidly by the large tempera- 
ture difference near the bottom of the tower. 
The rapid cooling of the gas below the liquid 
temperature indicates that the interface was 
cooler than the liquid. This suggests a measurable 
resistance to heat transfer in the liquid film. 

McKinzey [6] made some runs in which the 
same apparatus was used for cooling and dehumi- 
difying very hot combustion gas. His results are 
shown in Fig. 4 for comparison of the coefficients 
for the tetrachloroethylene- air system with those 
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for the water-air system. The coeflicients are 
approximately proportional to the square root 
of the vapour diffusivity. The exponent is not 
the 
velocity as shown by the divergence of the lines 


constant but increases slightly with mass 


of kg ay, at high gas velocities, 


Discussion 


Errors in the calculated coetlicients are caused 
by errors in the approximations made in the 
method of calculation, as follows : 
the actual 


volume of packing, but there was actually some 


1. The coeflicients are based on 


and mass transfer above and below the 


packing. The 
head was small because of the low velocity of the 


heat 
local transfer at the distributor 
gas and liquid streams. McApams et al. [7] and 
Hensev and Treypat [8] show that the assump- 
tion that the packing has the same area for heat 


and mass transfer is correct only when it is com- 


pletely wet. When the inside surface of a packing 


ring is partially dry and the outside is wet heat 
is transferred through the dry wall by conduction 
as well as across the liquid—gas interface. In the 
present the tile 
bundle of vertical flues and it is not likely that a 


study packing resembled a 
part of the surface was dry as is apt to be the case 


in randomly packed rings. The packing was 
observed visually through the windows and was 
thoroughly wet at the liquid rates used. 

2. Fog formation did not have a noticeable 
effect on the individual enthalpy transfer coetli- 
cients shown in Table 1 or on the over-all coefli- 
cients for the water-air system at high tempera- 
tures. This conclusion was also reached by 
Mizusmina and Nakasima [3] from experiments 
They found 


that the same enthalpy transfer coetlicient applied 


in a double tube cooler-condenser. 
for all three cases: no fog, fog formation only 
in the film and fog formation in the main gas 
stream. 

3. For calculating the moditied enthalpy, the 
exponent on the Schmidt and Prandtl groups was 
assumed to have the conventional value of 2/3. 
This value was verified experimentally in a wet 
and dry bulb apparatus [5]. For 28 runs on the 
tetrachloroethylene—air system the exponent n 
varied between the extremes of 0-51 and 0-74 


and W. 


B. ReracucKk 


with an average value of 0-64. The exponent was 
calculated from 


Se 
(pr 


where the wet and dry bulb temperatures ty, 


166 A Pwr 


(29) (0-24) (ton ” 


two) Pom 


and ty», are the experimental data. The equation 
is derived by equating Colburn’s factors for heat 
and mass transfer. The numbers 166 and 29 are 
the molecular weights of tetrachloroethylene and 
The 


taken as 


air, and 0-24 is the heat capacity of air. 
Schmidt and Prandtl! were 
2-15 and 0-74 in all of this work. Different coe- 
ficients are calculated when the exponent is 
reduced to 4, but they still do not fall along a 


groups 


single line in Fig. 4. 

The negative value for h, a, in Run 27 is 
impossible because it implies that heat is trans- 
ferred against the temperature potential. 

4. At 
enthalpy nor the modified enthalpy represents 


high vapour concentrations neither 
the transfer potential accurately because the 
true potential causing transfer is proportional 
to Ap /pyy. When the vapour pressure and the 
concentration are low py approaches unity and 
AH is proportional to Ap. Criss and Neson 
{9} have evaluated the error at high vapour con- 
centrations and give the correction as a function 
of temperature, 

5. In deriving the basic equations it was 
assumed that the values of C, and of % at the 
interface are the same as those in the main gas 
stream. This is an approximation for which 
the correction factor is not a simple function of 
temperature. It further that 
Pam Pa is unity. This is justified over the range 


was assumed 
where vapour concentration can be used for the 
concentrations, 
must be 


transfer potential. At higher 
the correction of Crips and NELSON 
applied. 

6. The experimental 
sensitive to errors and approximations in the 


values of hy ay are 


graphical calculation because an error in the 
slope of the tie line affects the value of h, ay 
more than kg @y,, particularly at high values of 
the slope. Therefore h, was calculated theore- 
tically to check the experimental values, At any 
point along the falling film of liquid, the point 
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value of the liquid phase heat transfer coeflicient 
is 

B.Th.U, 
hr ft? °F 


h, (dt dx), o (8) 


tt.9 —~ tum 


where 


it 
| i 


the 

gas-liquid interface 

ky, — the thermal 
liquid, 


temperature gradient at the 


re@ 


conductivity of the 


tutk the integrated average temperature 


over the thickness of the film. 


The experimental coetlicient is identified with 
the value obtained by integrating the point coe- 
This 


is also equivalent to integrating over the life 


flicient over the height of the falling film. 


span of the falling film, which is the time for the 
film to descend a single block of packing 4-5 in. 
high. The assumption here is that the liquid is 
completely mixed before forming a film on the 
next lower block of packing. The life span was 
of the 


s 


calculated from the surface velocity V 
falling film {10}. 


(9) 


where ya, and p, are the viscosity and density 
of the liquid and y is the total wetted perimeter 
HO ft 
includes the wetted walls of the tower. 


in this case, and 
Then 


in 1 ft? of tower section. 


for the conditions of Run 51: 


108 F. 
wy — 1-75 1b mass ft hr 
pp = 99 lb ft. 
k, — 0-073 (B.Th.U hr.) ft? (F ft). 
Vo 0.54 ft ‘sec. 
Film thickness 


Life span of film 


ft, (average) 


0.00040 ft. 


0-695 sec. 


Now if the temperature profile across the tilm 


is known at any instant during the life span of 


the film, the instantaneous (point) value of h, 
can be found from equation (8). The profiles 
were calculated by Schmidt's method of finite 
increments, which requires that the quantity 


(Ar)? /22A4@ be set equal to unity, where 


Ax = the increment of film thickness, 4-0 » 10°75 ft 
or one tenth of the total thickness 
ky 
L PL 


the increment of time for a single step, 


8-18 10 * sec, giving a total of 850 steps. 


In using equation (8), (dt/dr), . was approx- 
10” F ft 
constant 


imated as t, 4 —t, t and 


f, 9 Was given an value. 


arbitrary 


Assuming a constant value for} ¢, » will not 
affect the value of the calculated coeflicients [11]. 
The point value of h, was calculated for every 
other step, giving 425 values whose arithmetic 
mean is identified with the experimental value. 


Here is a selection of the calculated coeflicients : 





Coefficient 
Elapsed fraction int grated 


of the life span 


Point coefficient 
(B.Th.U hr) ft? °F over elapsed life 


span 


WoL 70 
Ont 500 
ory Wid 
Os54 His 
14MM) . tho 





The calculations were done on an IBM 650 digital 
computer, 

We ignored the effect of a velocity gradient 
across the film, but the error is small because the 
temperature profile flattens out rather quickly 
when a constant value is assumed fer t, 9. The 
tlow was laminar since the Reynolds number was 
only 120. As a final result the integrated value 
of hy dy is 440 & 40 = 17,600 (B.Th.U /hr) ft® F. 
Similar values would be obtained for the other 
the 
constant between runs and the viscosity has a 


runs because liquid velocity was nearly 
small effect. The lower experimental values for 
h, ay must result from experimental error and 
approximations in the method of calculation. 

7. The gas phase coeflicient was also calculated 


theoretically, from a diffusional model [11]. 


2x 29P Dv 
kg = _ - 
RI aZ 
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velocity of gas through the slots in the 
packing 

absolute temperature 

system pressure 

molecular weight of air 

diffusion coetticient 


for an arbitrary run in which 
G = 300 (Ib hr) ft® 
T = 550° Rankin 
2-74 ft sec 
0.0770 em? sec 
1 atm 
Z +-5 in., for a single block of packing 


The caleulated coeflicient is 


lb vapour 


‘- - 
(Ib vapour Ib air) hr ft® 


- B.Th.U, 
**(B.Th.U Ib air) hr ft? 


so that kg ay, ts 7-2 & 40 <= 288, 


Equation (10) applies only when the gas is in 


non-turbulent plug flow. In all of the runs the 
gas velocity was low enough to prevent turbulence. 
Plug tlow is more likely to occur when the gas is 
being cooled, which corresponds to the middle 
line of Fig. 4. From the middle line at G 300, 
the value of kgay is 270. Such a close check 
may be fortuitous because in the calculation it 
was assumed that the gas is remixed between 
the blocks of packing. Nevertheless, the coefli- 
cient will be lower when the profile has curvature 
because even when the profile is flat, the vapour 
barely penetrates to the centre plane of the slots 
during the passage through a single block of 
packing. This may explain the lower coeflicients 
obtained experimentally when the gas was heated. 
A mathematical model has been designed to 
explore the effect of the velocity profile, and this 
will be presented in a later paper. 

Comparison wire Orwer Mernops 


or CALCULATION 


The graphical calculation described here differs 
from the graphical methods described previously 


in that an assumed. or a corrected. condition 
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line is used to start each trial calculation. 
Mizvusuina and Koroo [2] and Lewis and Wurre 
[4] construct the condition lines directly from an 
assumed slope of the tie line. Consequently, 
they eliminate the preliminary calculations before 
the actual of the condition line 
starts. The present method requires the use of 
Fig. 2 and the 1/™%-— auxiliary line in Fig. 1 to 
start the calculations. However, fewer operations 


construction 


are required for constructing the condition line. 


DesigN or a Cooier-ConpdENSER 


An important application of the modified 
enthalpy potential method ts illustrated by the 
design of a cooler condenser for the recovery of 
vapours from an inert gas. 

Consider the design of a direct-contact cooler 


condenser for which the following are specified : 


1. Temperature, vapour content and flow rate 
of the entering gas. 
Vapour content of the exit gas. 
Temperature of the inlet liquid. 


The exit gas will have a temperature slightly 
above that of the entering liquid ; the difference 
will be indicated by the first trial calculation. 
The exit temperature of the gas, therefore, is 
known approximately as well as its enthalpy 
which is determined by its vapour content and 
temperature. Thus, only the liquid flow rate 
or the temperature of the exit liquid is unknown. 
The flow rate should be sullicient to wet the 
packing completely, while the temperature of the 
exit liquid must be low enough to cause de- 
humidification at the bottom of the tower. Within 
these limits either the flow rate or the temperature 
of the exit liquid is assumed and the other 
quantity is fixed. 

The graphical design procedure is illustrated 
in Fig. 8. With the assumptions made so far it is 
possible to plot the enthalpy operating line, 
A, — A,, the terminal points on the ¢ 
ing line, B, and B,, and the terminal points on 
the condition line, C, and C,. The shape of the 


operat- 


condition line is assumed and the corresponding 
é-operating line is determined by the method 
described under Turorericat Discusstox., The 
coetlicients for the liquid and gas films are known 
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, 8Tu PER (6 INERT Gas 


¢. OR | 





TEMPERATURE OF GAS OR Liquic 


Fic. &. 
cooler condenser using the moditied enthalpy potential. 


Graphical construction for the design of a 


for the chosen operating conditions and likewise 
the slope of the tie lines, The 
1/%- auxiliary line is laid out on Fig. 8, and a 


hy an ke Uy: 


trial condition line is constructed to some point 
C,, with the same enthalpy as point C,. The 
vapour content of the gas at point C,,’ is higher 


than the assumed value if C,’ lies to the left of 


C,,, and vice versa. Fig. 9 illustrates these possibili- 
ties and the situations causing them. A new exit 
liquid temperature is assumed for the construction 
of another condition line, the shape of which ts 


estimated from the one just constructed. A new 


¢, OR \, BTU PER LB INERT Gas 


TEMPERATURE OF GAS OR LIQUID, “F 
Fic. 9. Graphical design of a cooler condenser, 
showing the effect of varying the effluent liquid 
temperature. 


¢-operating line and a new 1/¥”% ~ auxiliary line 
are used, 

It is not necessary to retain the temperature 
of the exit gas assumed originally since only 
the trial 
condition line terminates at a point representing 


vapour content is important. If a 
a higher gas temperature and a lower vapour 
content than those assumed originally, this is 
the final this 
condition line is used to determine the ¢-operating 
the 1/e 


new temperature for the exit liquid 


accepted for calculations and 


line and auxiliary line for the next 
trial. <A 
need not be assumed, but the liquid flow rate 
must be adjusted to close the heat balance. 
In the final calculation it is only necessary for the 
constructed condition lines to 


assumed and 


coincide reasonably well and terminate at a 
point representing a satisfactory vapour content. 
The tower height is found by graphical integ- 
ration of 
G [ di 
Z (11) 
keay J dq — %% 


The term, by 
the tie line. 


¢;, is the vertical component of 


The design should be checked to ensure that no 
fog will form in the condenser. It is possible for 
fog to form at the gas-liquid interface or within 
the gas film even when supersaturation does not 
occur in the gas. This is possible because both 
temperature and concentration gradients exist 
across the gas film. Fog can form within the film 
when the profiles of actual temperature and 
saturation temperature intersect or become tan- 


gent within the film. Jounsvone et al. [12] give 


equations for the fog limitation lines such as the 
one shown in Fig. 9. These are the loci of all bulk 


gas conditions which cause incipient fog forma- 
tion at the interface. If the conditions of the gas 
lie above this line fog can form ; if they lie below 
the line fog cannot form. For constructing a 
condition line on Fig. 9 conditions at the interface 
and in the main gas stream are defined at all 
levels in the tower. If any point on the condition 
line lies above the fog limitation line constructed 
from the corresponding point at the interface fog 
“an form at the interface. 

Fog formation will always occur first at the top 
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of the tower. If there is no fog here there can 
be none elsewhere, either in the gas film or in the 
main gas stream. The fog limitation line shown 
in Fig. 9 is the critical line for the interface 
temperature at the top of the tower. If the 
condition line terminates below the line no fog 
the 


terminates above the line the distance between 


will form in tower. If the condition line 
the line and the point representing the exit gas 
conditions is a measure of the potential to form 
fog. When the potential is small fog can be 
prevented by raising the temperature of the inlet 
gas, as shown in Fig. 9. If the inlet gas is heated 
at constant humidity its new conditions may be 


represented by some point, such as C",. The 


new exit gas conditions will be represented by 
some point, such as (",. 
point C* 


The displacement of 
is smaller than that of C",, and con- 


n 


siderable heating may be required to avoid fog. 


and W. B. 


ReTALLICK 


The fog limitation lines depend very much on the 
liquid temperature so that sometimes fog may be 
prevented by a slight increase in the inlet liquid 
temperature, This increase must be small because 
a low inlet liquid temperature is necessary for 
thorough dehumidification. 

Fog will not necessarily form in the gas film 
when supersaturation exists. For a small degree 
of supersaturation the presence of condensation 
nuclei is necessary for a finite rate of fog forma- 
tion [12]. Bras [13] points out that in the 
absence of nuclei a large supersaturation is 
necessary for a finite rate of fog formation. For 
water vapour the rate of nucleation is negligible 
when the dimensionless ratio p,/p_ is less than 3. 
The rate increases by a factor of 10° when p,/p, 
is increased from 4 to 5. Furthermore, small film 
high turbulence the 
tendency to form fog in the film. 


thicknesses or decrease 
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The unrestricted engineering Bernoulli equation 
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Abstract—Based on a force analysis, Birp recently derived restricted * engineering Bernoulli 
equations." The procedure consisted of forming the scalar product of the velocity with the 
Navier-Stokes equations of motion and then integrating over the volume of the open system to 
obtain first a general mechanical-energy balance. For the restricted cases of isothermal and 
isentropic flow at constant composition (or else at chemical equilibrium), and with each species 
experiencing the same conservative body forces, Birp converted the mechanical-energy balance 
into “* engineering Bernoulli equations.” In the present paper a Bernoulli equation will be derived 
by a generalized method ; the only condition will be that of identical body forces upon the con- 
stituents. * 


Résumé—Birp a récemment obtenu une équation réduite qu'il appelle “* Engineering Benoulli 
équation ” basée sur une analyse des forces. Le procédé consiste a faire le produit scolaire de 
la vitesse avec les équations de mouvement de Navier-Stokes, et d'intégrer ensuite au volume du 
systéme ouvert pour obtenir un bilan général de énergie mécanique. Dans le cas particulier 
d'écoulement isentrope et isotherme a composition constante (ou a léquilibre chimique) Birp 
expérimente pour chaque espéce les mémes liaisons internes et convertit le bilan énergie mécanique 
en “ équation mécanique de Bernoulli.” Dans la présente étude on obtient une équation de 
Bernoulli par une méthode généralisée, la seule condition étant que des forces identiques s’exercent 
sur les constituants. 


Zusammenfassung— Auf der Grundlage einer Vektoranalyse hat Bird kiirzlich spezielle 
“engineering Bernoulli equations” abgeleitet. Diese Ableitung bestand aus der Bildung des 
skalaren Produktes aus der Geschwindigkeit mit der Navier-Stokes’schen Bewegunygsgleichung 
und anschliessender Integration iiber das Volumen des offenen Systems, wodurch zunichst 
eine allgemeine mechanische Energiebilanz erhalten wurde. Fiir die auf isotherme und isentrope 
Stroémung beschriinkten Fille bei konstanter Zusammensetzung (oder auch bei chemischem 
Gleichgewicht) und in jedem Falle, in dem dieselben Massenkriifte ausgeiibt werden, formte 
Bird die mechanische Energiebilanz in * engineering Bernoulli equations * um. 

In der vorliegenden Arbeit wird cine Bernoulli-Gleichung mit Hilfe einer verallgemeinerten 
Methode abgeleitet ; die cinzige Bedingung ist die der gleichen Massenkriifte auf die einzelnen 
Komponenten. 


*The treatment is readily extended to the case where the body forces are not identical: In Birp’s equations (20) 
and (21) replace pF, by 2p; Fj = 2p; (Fj, + Fj) where Fj, are conservative body forces [equal to — V4; (x, y, 2)), 


and F,, are time-dependent body forces. Then, the integration | [v - (2p; Fj,)] dv may be carried out in a manner 
entirely analogous to that of Birp (his equation 26). This results in 

, iE 

[ v-zmRoay =» 


+. 
Y Y 


peg [(V— Vy) nn] dS — Z oj (jj-nydy 
y, 


z+ | r¢,4v + 2 i Vepdv + | v-(2p, Fav (£1) 


. . . 


y y a 


where ¢; is the potential energy of the i-th constituent, the energy resulting from the location of ¢ in the potential field 


¢;(#. y, 2). Also Ey = 2 po, dy and e, = (1/p) 2p dj: 1; is the rate of production of constituent i, per unit 
j 

volume. The last four terms of this equation are zero when all the body forces are identical and conservative; inclusion 

of these terms renders the equations to follow unrestricted. By rewriting equation (f.1) as an expression for dE, /dt, 

each of the other terms may be rationalized into a specific contribution to the total change in E,. 
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MECHANICAL-ENERGY 
AND BERNOULLI EQUATION 


BALANCE, 


Tue general mechanical-energy balance developed 
by Brrp [1] (the integral of his equation 21 with 


equivalent to 


dE, : dE, _ 
dt “dt 
| p (¢, i eg)[(v Vy) , n| dy 


of . 
liv. v par | [(> v) n|dy 
, a 


+t (Tt: Vv) dy (l.la) 


. 


’ 


This balance may be transformed into one that 


facilitates solution. Rewrite, using the Gauss 
divergence theorem, 


1 (Vv - Vp) dy 


y 


| «v pyv)dy¥ 4 [ piv. war 


! [ piv way 


y 


| (pv -n) dy 
+ 


In the equation 


| [piv 


[ (ow may -Vy)+n}dy 


>. . 
¥ / 


ah | (pVy n)df (1.2a) 
y* 


the first term on the right-hand side is the pressure 
flow-energy* transferred across the boundary of the 
system ; the second term is the work rate resulting 
from pressure forces on the system boundary 
(and is designated dW), /dt). Similarly, the right- 
hand side of 


*Birxp (and others) often call equations (1.2 a, b) the 
pressure work and the shear work, respectively ; in this 
paper all forms of energy transferred across the boundary 
because of material flow are defined flow energy. This 
procedure and most of the symbolism follows the termin- 


ology of Osrrr [2). 


[lc -v)- nj] df = 


| {[r-(v —vy)]-n}dy 
gy 
t [(r Vy): n| df (1.2b) 
¥ 
consists of the shear flow energy and the work 
rate of shear (dW,/dt). Thus, equation (1.1a) 
vields another mechanical-energy balance : 
dk, | dE, dW, dw, 
dt dt dt dt 


(t: Vv) dr 4 [ piv -v)d¥ 


f » ' 
[ tote, T” €.4 


7 
Vy): n] 
g ' |r (v 


~p p)[(v 


vy)|-n}dy (1.1b) 


Yet another variation of equation (1.1) is obtained 
if, using the equation of continuity dp/M + V- pv 
0 and the definition of the Stokes derivative, 
one rewrites p(V ~ (p/p) (Dpe/ Dt). 
The extension of Birp’s procedure to obtain 


V) as 


an unrestricted Bernoulli equation depends solely 


upon the evaluation of the integral fo Vp) dr 


’ 
in equation (I.la). A generalized procedure will 


be treated here. With the thermodynamic 


relationship 
v dp dh T ds 


and the definition of the Stokes derivative, 


[ ww - Vp) dy 


2 My dx;, 
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By the three-dimensional Leibniz 


theorem, 


invoking 


| (v- Vp) dy 


| ph(vy,-n)ds 

~ 

[ (v. Vaya + 4 [ ly 
) vif ; - d¢ 
, dt. I 


’ ’ 


| piv, - nd (a) 
$ 


The continuity equation, along with thed iver- 


gence theorem, allows 


j ) 

{| a p(w Wh)| ds 
. ot 
~ 


= |[—4(V- pv) — p(v- Thy] dy 


| (Vv. pvh) dy 


5 


| ph iv -n)dy 


. 
¥ 


° 


| pu dy 


y 


1 
h) dy : 
pn) at 
dW, 


Vy: ly = 
P\ a dt 


gy 
With these 


: [ wv. vp) ar [o|r 


substitutions in equation (a), 
Ds > Da; ' 
aise Su; —| dv 
pi * a] 

| ph|(v —vy)-n]dy 
g 

dU - dW, 

dts dt 


(1.3) 


Upon introducing equations (1.3) and (1.2b) 


into equation (1.la), the general equation is ob- 
tained 

dE 

dt 


+, + e4)[(V —Vy)-n| 


vy)|-n}dy (r+: Vv)d¥ 


. 


y 


aw,  aW, 


peste 1.4 
dt dt ( 


= | PR; 


¥ 


And this represents the unrestricted “ engineering 
Bernoulli equation ’’—-a generalization of those 
presented by Birp 
science strives to formulate descriptions of the 
most general situations (this equation does have 
the condition that the body forces on each species 
are identical and conservative, but see footnote 


which is desirable since 


on p.1). 

Equation (1.4) is readily restricted to the iso- 
thermal and isentropic constant-composition (or 
chemical equilibrium) flow treated by Birp. To 
specialize equation (1.4) for the case when temp- 
erature is independent of position, consider first 


the term 





— s(V-pv) | dy 
With the divergence theorem, the continuity 


equation dp/dt v. and the Leibniz 


theorem, 


pV; 


r | ps(v-n)dy 
gf 

T : [ os av + T ps (Vy- nj df 

dt , Z ‘ 

, ¥ 

of dS 
—7 | ps[(v —Vy)-n]d¥% —7 + 
¥ 
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Upon substituting in equation (1.4), 


| {ech Ts + & Heg[(v Vy) n| 
¢ 
t [r (Vv vy)} nidy¥ 4 | (x: Vv)ady 
dW, F, 
dt r 


Vé 


Tia) (1.5) 


Wis the mechanical work of 
T's) equals the Gibbs free 


where W,, — W,, 
the system. Also, (/ 
energy g. Since (Vv Vy) Is zero except at a flow 
boundary of the system, the surface integral 
represents the transport associated with material 
For the 


case with one flow stream of constant composition 


flow across the boundary. isothermal 
(or chemical equilibrium) and negligible change 
in shear flow energy, equation (1.5) reduces to 


dd dk, . dk, 
dt dt dt 
- Ud 
constant 


A [ rity (g + & es)| r 
' iw _ dD 
| (tr: Vw) dy all PITS Ju 
. dt nt 

’ Anmpe, = 0 


(1.6) 





TS is the Helmholtz free 


energy (and this form is equivalent to Birp’s 


where A Uw 


equation 28). For isentropic flow equation (1.4) 
becomes 


iE ; 
5 | 1p (h 


¥ 


+e +e) [(V —vy)-n] 


-[t-(v—vy)]-n}dy 4 | (+: Vw)dy 


y 


dW, 
dt~ 


Dr, dy F; v¢ 


s constant 


(1.7) 


+2 | PH, 


’ 


And, with one flow stream of constant composition 
(or chemical equilibrium) and negligible shear 
flow energy change 


. Gaceru 


Vo 
constant 


A[ rity (h ey + ¢,)] 


dW,, 


di Am; e, 0 


[ (+: Vv) dv — 





And this form is equivalent to Brrp's equation 


(29). 


2. FRIcrion EQUATION 


As a matter of interest, equation (1.4) can be 
transformed into what might aptly be called the 
Friction Equation. Notice that 


De, 
Dt 


Do, 


BK; at BK; " 


Dp 


PR, ' Wt 


D 
rT py v- Vp, By; in 


With the continuity equation for constituent ij, 


Ve py V - jj. 


where (j;/p;) — (u; — V) is the difference between 
the velocity of constituent 7 and the mass-average 
velocity (Ref. [3], pp. 295-6, 344-5), 


Da, 
Dt 


Dp 


V - x 
Ki i, iS Dy 


PH, ri wy Hi pi (VV) —~ pw 


Zon, Da, "7 


-Em+B 
Dt mie VAs 


= "% My 


“¥ 
ot 


c 


pg (V -v) gv Vp 


= "By = wyV . i, 
g(pV -v —V-pv +v- Vp) 
LV wi + Li Vay (a) 


& ify 


Ts;. 


partial specific quantities 


Also, with p; = Ah where fh; and s; are 


! 1 
7p hide VT pr hide VT + 2s,j,-VT =0 


Subtracting this from equation (a), and noting 
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Tj, - V (u,/T) 
Da 

Dy 
~ PED 
+ LY Tj, - V (n;/T) 


i, Va; — O/T) wi, VT. 


2 3,j,-VT 4 . . (b) 


Finally, substituting (b) in equation (1.4), and 


rearranging, 


dE dW, 
dt dt 


EY (Twi) a 


2 [(h, T)i,: VT 


Tj, - V (ue, T)| + fr; y;} dy (2.1) 
The integrand of the right-hand side of equation 
(2.1) here called the local friction rate, f. is equal 
to T times the rate of entropy production per 
unit volume (Hirscurenper, et al. [4], p. 702).* 
The integral can be interpreted as being the total 
friction of the irreversible flew process. The first 
term represents the degradation of mechanical 
energy into thermal energy — the mechanical frict- 
ion accompanying fluid motion, The second term 
is the friction from heat transfer through tem- 
perature gradients —an evaluation of the mechan- 
ical energy which could have been obtained had 
the heat passed through reversible cycles rather 
than through irreversible temperature gradients. 
The third term is the friction of diffusion and the 
last is the degradation of chemical energy because 
of irreversible reaction. 

*The term n;V,;-X; = 2j,- F, of Ref. 4 is zero here, 
since all body forces were assumed equal, and 2j; — 0. 
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NOTATION 
Helmholtz free energy of mass contained 
in the system, at time ¢ 
energy (internal + kinetic + potential) of 
mass contained in the system at time 
flow energy per unit mass at position 
(az, y, =) and at time ft 
kinetic energy per unit mass 
kinetic energy of mass contained in the 
system, at time ¢ 
potential energy per unit mass 
potential energy of mass enclosed, at time ¢ 
shear component of flow energy 
local friction per unit volume and unit time 
Gibbs free energy = (1/p) 2 pig; 

(1/p) 2 pip; 
enthalpy at (7, y, =) and at time ¢ 
contribution of the i-th constituent to 
h(a, y, 2,0): the partial specific enthalpy 
flux of relative to 
mass-average flux 


y, =, t) 


j; (*. y. =, 0 = mass constituent i 
m (t) 


my (t) 


mass contained in the system at time ¢ 
mass flow rate 
unit vector normal to a surface pointing 
outward from the enclosed volume 
pressure 
= entropy per unit mass 
contribution of the i-th species to s(2,y,2,¢): 
the partial specific entropy 
entropy of mass in the system at ¢ 
surface of system at ¢ 
time 
T (a, y, temperature 
internal energy per unit mass 
internal energy contained at ¢ 
U; (a, Y, 2. absolute velocity of the i-th constituent 
Vv (a, y. 2, velocity relative to (7, y, =) att 
volume of system at time ¢ 
the total work of the system from time zero 
co-ordinate fixed in space 
concentration of the i-th constituent at 
(a, y, z) and ft; 2; (p;/p) 
Greek Symbols 
AN = the difference between the quantity X at 
the entrance and exit of an open system 
with one flow stream 
heat flux at (7, y, z) and time ¢ 
density 
chemical potential of the i-th constituent 
at (av, y,2) and ¢: partial specific Gibbs 
free energy 
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tT (2, y, =, ) = shear stress tensor (see Ref. [4], p. 699) vt shear 
¢; (2, y, =. 1) = potential for i-th constituent ¢é potential 


Marks above symbols 


Subscripts 


i-th constituent pert meme 
. . . z t » ~ “ . 
kinetic derivative with respect to time 


Stokes (substantial) derivative 
mechanical 


pressure Df (a, 4. = 0) yf 
surface Dt 
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Thermal cracking of acetone—lI. 


Influence of the process variables on the product distribution 


G. Fromenr, H. Pucke and G. Gorrnats 
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(Received September 15, 1960) 
Abstract—The thermal cracking of acetone into ketene was studied by flow-experiments. The 
effect of the feed-rate on the product distribution was investigated in the temperature range 
710-750 © and for residence times below 1 sec. The conversions and selectivities are correlated 
vs. the total acetone conversion in a diagram valid for the whole temperature range investigated. 


Résumé 


régime dynamique. On a déterminé principalement Vinfluence du débit d’alimentation sur la 


La décomposition thermique de lacétone en cétéene a été étudiée par des essais en 


distribution des produits a des températures allant de 710 a 750 °C et pour des temps de séjour 
en-dessous d'une seconde. Les différentes conversions et sélectivités ont été portées par rapport 
& la conversion totale de lacétone dans un diagramme valable dans le domaine entier des temp- 
peratures considérées, 


Der thermische Krackprozess von Azeton ins Keten wurde durch Fliess- 
Der Eintluss der Dosierung auf die Produkt-verteilung zwischen den 
und fiir Reaktionszeiten unter eine Sekunde, wurde unter- 


Zusammenfassung 
experimente studiert, 
Temperaturgrenzen von 710-750 C, 
sucht. Der Zusammenhang zwischen Umwandlungen und Selektivitéten im Verhdltnis zu der 
totalen Umwandlung von Azeton ist in einem, fiir den ganzen untersuchten Temperaturbereich 
giltigen Diagramm, dargestellt. 

1. INrropuction presence of ethylene and of part of the carbon- 
Acetic anhydride is used in large quantities in monoxide has been explained by the apparent 
the cellulose acetate plants and is produced by over-all reaction [1] : 
reaction of ketene on acetic acid. Until recently 
the 


In recent years 


CO + CH, 


2CH, CH, + 2CO (2) 
ketene has been obtained exclusively by 
or the intermediate 


Whatever the mechanism 


steps of this reaction may be, the decomposition 
of ketene is a fact, that is reflected in the very 
short residence times and the consequent low 


thermal cracking of acetic acid. 
however increased attention has been paid to 
the use of acetone as a starting material as a 
result of its present surplus. 


When submitted to pyrolytic conditions acetone conversions per pass (< 25 per cent) that have 


decomposes according to an endothermic apparent been adopted in industrial practice [1]. 


over-all reaction : 


CH, CO CH, + CH, — CO + CH, (1) 
that may be considered irreversible in the range 
of practical interest (700-750 C). Ketene and 
methane are not the only products, however: in 
the rang considered ethylene, carbon monoxide 
dioxide, hydrogen and carbon are also 
The mechanism of formation of these 


and 
obtained. 
byproducts is far from being elucidated. 


The formation of hydrogen and of cokes has 
been attributed [1] to the reaction 


CH,COCH, - 3H, + CO + 2C€ (3) 


which consumes acetone in parallel with (1). The 
prevention of this coke formation, which is 
strongly catalysed by the presence of nickel, has 
received considerable attention in the patent 
literature, 

If, for reasons of strength, the usual types of 
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stainless steel are preferred to the nickel free 
chromium steels, a special pretreatment that forms 
a carbide film [2] or the addition of small quantities 
of a sulphur containing compound, e.g. carbon 
disulphide [3], is required to reduce substantially 
the carbon formation. 

Practically no quantitative data have been 
published on the influence of the process variables 
on the rate of the reactions and the product 
The data available were obtained 


distribution. 
in static systems and under conditions that are 
not useful for industrial practice. It was not the 
purpose of the work reported here to study the 


reaction mechanism in detail: a chemical engin- 
eering approach turned toward design was adopted. 
A flow apparatus was used and the conditions 
of temperature and residence times were chosen 
to be as close as possible to industrial conditions. 
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2. Apparatus AND EXPERIMENTAL 
PROCEDURE 

A scheme of the apparatus is given in Fig. 1. 
The acetone is fed by gravity and under a constant 
head through a flowmeter to the pyrolysis-tube. 
This is a chromium steel (16 per cent Cr, no Ni) 
spiral-wound tube of 10 mm outside diameter, 
6mm inside diameter and of length 1-20 m, a 
detailed scheme of which is given in Fig. 2. This 
tube is placed in an electric oven, filled for the 
sake of equalizing the temperature with copper 
shots. Besides of the main heating wire, proper 
to the oven, secondary resistances are provided 
around the inlet and outlet sections of the tube, 
in order to compensate heat losses to the exterior 
and to reduce the preheating section to a mini- 
mum, 

The inlet temperature of the acetone vapours is 


|-Calibrated Mariotte fiasks 
2 - Stopcock 
3- Acetone fiow- meter 
4- Oven 
5 - Reactor - tube 
6-Cooler 
7-Pump 
8 - Thermostat (70°C) 
9- Transfo 
10-10 - Absorption - towers 
11-1 Traps 
12- Gas- flow - meter 
13-Pressure regulating device 














Flow sheet of the apparatus. 
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Fic. 2. Oven with reactor tube. 


the wall 
temperature of the tube is measured at six 
different spots by Ni-NiCr thermocouples. 


At the outlet the gases are drastically cooled 


measured by thermocouple 1, while 


in order to limit as much as possible the conversion 
to the reactor itself. Because ketene polymerizes 
rapidly in liquid acetone the temperature is not 
lowered beneath the dew-point of acetone in the 
gases. To be safe the lines were kept at 70° by a 


thermostat or by heating wires. The gases than en- 
ter one of the two absorption towers where they are 


scrubbed by a known volume of standard sodium 


hydroxide. The acetone is absorbed while the 


ketene reacts under the formation of sodium 
acetate. 

A trap filled with 0-01 N sodium hydroxide 
serves as an indicator that no ketene escapes. 
The volume of the remaining gases is measured 
in a wet gas meter, The pressure in the apparatus 
is controlled by a pump and the bottle 13. 

The unreacted acetone present in the sodium 
hydroxide liquor was titrated by Messinger’s 
the The 


ketene present in the gases was determined by 


method, based on iodoform reaction. 
titration of the excess sodium hydroxide in the 
absorption liquor. The gases were analysed with 
an Orsat apparatus, while additional information 
was obtained from gas chromatography and mass 
spectrometry. 

absence of carbon 


the nickel, 
deposition did occur and it was necessary to burn 


In spite of 


it off from time to time. After this operation the 
apparatus had to run some time before reprodu- 
cible data could be obtained. After this start-up 
period neither the conversion, nor the product 
distribution were dependent on time. 


3. DEFINITIONS 
The conversions and the product distribution 
of this complex system have been characterized 
as follows: 
The total conversion of acetone has been defined 
as: 
moles of acetone decomposed (4) 
moles of acetone fed 
The conversion of acetone in side-reactions (e.g. 
equation 3) as 
Vas 
moles of acetone decomposed by side-reactions (5) 
5 
moles of acetone fed 
The conversion of acetene to ketene, according 
to (1) as 
moles of ketene formed 


Y. . 6) 
- moles of acetone fed \ 


The conversion of ketene itself, according to (2) as 


moles of ketene decomposed 


. 
kal 
. moles of acetone fed 
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ketene concentration, measured 


The 


effectively, is also obtained through a conversion, 


net as 


defined as 


moles of ketene measured 
Dice : ; (8) 
moles of acetone fed 


(4). (5) and (6) are related by 


(6), (7) and (8) by 


Vak lke "Ka 


Furthermore, we define the yield by 


the acetone selectivity by 


moles of ketene formed 1 (12) 


moles of acetone decomposed 


Table 1. 


Experi 


the ketene selectivity by 


moles of ketene measured — a, 


(13) 


moles of ketene formed Tek 


the measured amount of 
acetone fed and from titration of the acetone in 
the absorption liquor; 2, from the titration of 


NaOH in this liquor; 2,4 from the ethylene 


r,, is obtained from 


content of the gases and their volume, according 
to reaction (2); 7, from (10); 2, from (9). 
Resu.ts 


4. EXPERIMENTAL 


The results of the series of experiments at 710, 
3, 
t and 5 respectively. The conversions are given 
in function of the time factor Vp, F, rather than 


730 and 750 C are given in Table 1 and Figs. 


of the residence time itself, which is a complicated 
function because of expansion. The Vy, F values 
are such that the residence times varied between 
0-1 and 0-6 sec. 


mental results 
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Fic. 5. Experiments at 750 °C. 


4 


5. Discussion 


Qualitatively the results agree with the picture 
given by equations (1), (2) and (3): 2, 
increase with Vy, F; a, tends to a maximum, 
that is best seen from the data at 750 C, where 
the conversions are highest and will probably 
decrease at higher Vy, F because of the increasing 
importance of the decomposition reaction (2). 

The attention should be drawn on the incom- 
plete picture of the consumption of acetone in 
side-reactions given by equation (3). 

Indeed, material balances on the acetone con- 
sumed in side-reactions may be obtained in two 
the 
sodium hydroxide content of 
liquor or from the CO content of the gases after 
deduction of the CO generated by the ketene 


ways: either from acetone and residual 


the absorption 


decomposition. 

These balances do not check, although the 
apparatus was controlled for tightness before 
every experiment. Clearly then, equation (3) 
does not account for all the acetone decomposed 
in the side-reaction, either because other reactions 
also occur or because it is only the last of a series 
of steps which have not come completely to their 
achievement. 

For the experiments at 750 C for example, the 


percentage of acetone consumed in side-reactions, 
based on the content of the gas, showed a very 
detinite trend and increased from 28 per cent of 
17 


This seems to 


the value obtained by titration at V,_ F 
to 70-4 per cent at Vp, F = 71-50. 
support the latter hypothesis. 

It is seen also how the curves of the Figs. 3, 
t and 5 do not extrapolate through the origin. 
This is characteristic for homogeneous reactions 
carried out in a flow reactor and results from the 
fact that not all of the volume accounted for is 
at the reference temperature considered, This 
will be discussed in detail in Part IL, where the 
curves will be corrected for this effect. 

In a first attempt to eliminate this effect and 
also to relate the different conversions to one 
another, these have been plotted in Fig. 6, together 
with the vield and the selectivities, vs. the total 
acetone conversion. 

It is seen how «,,. 1, and consequently » may 
be correlated as unique functions of 2,,. The dia- 
gram indicates clearly how the vield » decreases 
as the acetone conversion increases. The diagram 
is limited to the ascending part of the curve 
Ve Vy F. 


considered a very sharp decrease in the 2, is 


If values past the maximum were 


to be expected. As is to be expected also from a 
of the ketene 


selectivity = decreases as x,, increases and = vs. ay) 


system consecutive reactions, 


~ 


is very nearly independent of the temperature. 
A closer look however shows how the values at 
7 C 750 CC. 
This seems to indicate that the activation energy 


are somewhat below those at 


of the ketene decomposition (2), although not too 
different from that of the formation (1), is lower. 

The y. 
710 C are also slightly higher than at 750 °C, 
Nevertheless, these differences are very small and 
the diagram, which covers the field of interest 


values of the acetone selectivity, at 


for industrial practice, is therefore believed to be 
general, i.e. not limited to the apparatus used 
with its characteristic temperature and conver- 
sion profiles. 

Presented as they are the results permit the 
prediction of the product distribution when the 
total of is known. The 


conversion acetone 


prediction of this conversion is a problem of 


kineties which makes use of w vs. Vp F data. 


8 
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Fic. 6. Conversions, yields and selectivities vs. the total acetone conversion. 


It is not possible however to derive a rate gus constant 
error-variance 


equation from the data as they are presented in 
- . at ete . temperature of the gases at a distance L 
the Figs. 3, 4 and 5, since they are too specitic 
a Be ; 4 entrance temperature of the gases 
for the apparatus used. The derivation of rate t, = wall temperature 
equations from experiments of the type reported T — temperature 
here will be tackled in Part II of this paper. reactor volume 
= ‘py = equivalent reactor volume 
NOTATION conversion of acetone to ketene 
= specific heat conversion of acetone in side reactions 
= activation energy total acetone conversion 
feed rate ketene conversion 
= thermal conductivity of the gas 2 net acetone conversion to ketene 
rate-constant 7 
length ' acetone selectivity —** 
= mass feed rate “at 
. ti ae 
= order of a reaction ketene selectivity "ke 
partial pressure Tak 
reaction rate = yield 
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Deriving rate-equations from non-isothermal flow experiments* 


G. Froment, H. Pucker and G. Gorrnats 
Laboratorium voor Organische Technische Chemie, Rijksuniversiteit Gent, Belgium. 
(Received September 15, 1960) 

Abstract—When homogeneous reactions are carried out in a flow apparatus a clear-cut separation 
between preheat and reaction sections is impossible, so that the latter shows in fact a longitudinal 
temperature gradient. This situation greatly complicates the derivation of kinetic data. With 
the help of the equivalent reactor volume concept and a proposed method for estimating the 
activation energy prior to the knowledge of the rate constants, the conversion vs. V,/F data 
of Part 1, which are too specific for the apparatus used, have been corrected. This enabled rate 
equations to be derived for the ketene formation and decomposition and from these rate equations 
for the net ketene production. Used in conjunction with the selectivity diagram of Part I this 
result permits the prediction of a fairly complete product distribution. 


Résumé 
définir nettement les sections de prechauffage et de réaction. Tl existe done en cette derniére un grad- 


Lorsque des réactions homogénes sont exécutées en régime dynamique il est impossible de 


ient longitudinal de température qui complique séricusement étude cinétique. A Taide du concept 
du volume équivalant et dune methode originale permettant d’évaluer lénergic d' activation, avant 
- Vy F de la 1° partie, 
trop spécifiques, ont été corrigées. On peut alors deériver des équations cinétiques pour la formation 
utilisé en conjonction avec 


méme de connaitre les constantes cde vitesse, les données conversion 


et la décomposition du céténe et done pour la production effective : 
le diagramme des sélectivités de la 1°°° partie ce résultat permet de prédire la distribution quasi 
compléte des produits du cracking. 


Zusammenfassung —Wenn in cinem Fliess-apparat homogene Reaktionen durchgefihrt wird, 
eine volistindige Trennung der Vorwirme- und Reaktions-scktionen ist unmédglich, so dass dic 
Reaktions-sektion in der Tat cinem longitudinalen Temperatur-gradient aufweist. Dieser Umstand 
kompliziert sehr die Ableitung von kinetischen Daten. Mit der Hilfe des Begriffes des aquivalenten 
Reaktions-volumens und ciner vorgeschlagenen Methode zur Abschiitzung der Aktivierungs- 
energie, bevor man die Geschwindigkeits-konstanten der Umwandlung kennt, 1, F Daten im 
ersten Teil, welche viel zu spezifisch fir den verwendeten Apparat sind, wurden korrigiert. In 
dieser Weise es ist méglich dic kinetischen Gleichungen fiir die Formung und Zerlegung der effek- 
tiven Keten-produktion abzuleiten. Verwendet zusammen mit dem Selektivitét-diagramm vom 
ersten Teil, dieses Ergebnis ermoglicht dic ziemlich gute vorherige Abschitzung der Krackprozess- 
verteilung. 


1. INrropUucTION the determination of the form of the rate equation 


Tue design of a chemical reactor requires in the 


most general case the simultaneous integration of 


the continuity, energy and momentum equations. 
Even if the physical properties are taken to be 
constant analytical integration of this system is 
excluded, because of the non-linearity of the equa- 
tions. Numerical integration however is straight- 
forward, 


The reverse problem, namely, the derivation of 


kinetic data is far more involved, This comprises 


*Part I: This issue, p. 173. 


itself and of a number of constants-for homo- 
the order. the rate constant 
which 


geneous reactions ; 
the 
temperature dependence of the rate constant. 


and activation energy, gives the 
This can be done by successive trials, postulating 
rate equations and trying values of the constants 
until the calculated conversion, temperature and 
pressure profiles reproduce the experimental. 
Even with the help of statistical methods and 


computers this is a formidable task however, 
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especially when no first approximations are at 
hand. The procedure most often used is than to 
simplify the physical picture. If, by reducing the 
size to a convenient scale, the pressure drop 
becomes negligible, then the momentum equation 
need not be considered and if it is possible to 
eliminate all temperature gradients, the energy 
The 


mathematical model has then become sufliciently 


equation need not be taken into account. 


simple to permit, in most cases, the derivation of 


a satisfactory rate equation. 

The condition of isothermicity is a very drastic 
one, however: the heat effect of the reaction is 
often considerable and cannot always be com- 
pensated entirely. 

Furthermore, in some cases it is inevitable that 
part of the reactor is used to bring the feed to 
the desired temperature. This is the case with 
flow reactors for homogeneous reactions, where, 
in contrast to catalytic reactors, no clear-cut 
separation between preheat and reaction section 
is possible. If the rate is to be determined at a 
reference temperature, say 7',, and if the reaction 
volume is counted from the point where 7, is 
the 
section, that cannot be excluded, is not accounted 
for. 
where, while being cooled, the Conversion COon- 


reached, then conversion in the preheat 


Similarly at the outlet there is a section 


tinues to some extent. 

This situation, from which there is no way out, 
complicates seriously the task of deriving rate 
equations for homogeneous reactions. The idea 
of this paper is to take the situation as it is and 
nevertheless to find a method for deriving at 
least approximate values for the rate constant, 
the order and the activation energy, the refine- 
ment of which, if necessary, would be left to the 
computer, 

The essential features of the method are the 
use of an equivalent reactor volume, as introduced 
by Hovcex and Warsow [2] and of a proposed 
short-cut procedure fora preliminary determination 
of the activation energy. This will be illustrated 
on the data of Part I. 

2. The 


EQUIVALENT KEACTOR-VOLLU ME 


The equivalent reactor volume has been detined 


as that volume, Vy. which, at the reference 


Part Il 


temperature T's: would vive the same conversion 


as the actual reactor, with its temperature profile.* 


It follows that 


from which 


| R | exp 


I 
, ue (1) 
1 


| Oa | 
R T 


Once Vy, has been derived the calculation of the 


rate constant is straight-forward, as for isothermal 
experiments and solely based on the continuity 
equation 


F da 


rd Vp 2) 


The calculation of V;, requires the knowledge of 
the temperature profile along the tube and of the 
activation energy E. 

Note that 
occuring simultaneously the dependence of Vie 


also where several reactions are 
on E leads to different values of Vie for each of 
the reactions considered. 

ACTIVATION ENERGY 


3. ESTIMATION OF THE 


However, as the activation energy ts not known 
a prior, the use of the equivalent reactor volume 
concept leads to a trial and error procedure: a 
value of E is guessed and with this value and the 
measured temperature profile, Vi, is derived from 
(1) by graphical integration. Then, from (2), the 
rate constant is derived. 

This procedure is carried out at several temp- 
erature levels and from the temperature depen- 
the 
Arrhenius’ formula, a value of EF is obtained. If 


dence of rate constant, expressed by 
this value is not in accordance with that used in 
the calculation of V;, the whole procedure has to 
be repeated with a better approximation for EF 
and thus for Vy. 

In this paragraph a short-cut method for the 
first estimate of E. not based on Arrhenius’ 
formula and therefore not requiring knowledge of 
the rate constant, is proposed. It is based on 
the following considerations. 


Consider two experiments carried out in a flow 


*The notation used is given at the end of part L. 
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reactor, one at a reference temperature 7',, the 
other at the reference temperature 7, and let the 
conditions be such that the temperature difference 
AT =T,— T, is the over the whole 
length of the reactor. The reaction taking place 
> B. If the 
feed-rates are adjusted so that equal conversions 
are obtained, then the x, i.e. the p, vs. Vp profiles 
Then 


same 


is homogeneous and of the type A 


are identical in both cases. 


= 


dx, 
dv,” ‘Pah 


(Pare 


at all points, while AT is independent of Va. 


From the continuity equation, 


E — 
F da x exp | =) (p,)"dVy 


applied to both experiments, it follows that 


fA 
1 


from which E = 23 R (3) 


This means that the activation energy may be 
derived different 
temperatures the rate 
constants, provided the conversions are the same 


from two experiments at 


without even knowing 
and the temperature profiles, plotted vs. V,, are 
parallel. 

The latter condition is not always fulfilled in 
practice. It requires that the heat effect of the 
reaction is negligible or entirely compensated for 
by the heat flux from or to the surroundings or 
(and) that the specific heat of the gases is very 
large. 

If the reactant A is consumed by more than one 
reaction than, at equal conversion to the product 
of interest, B, p, is only equal in both experiments 
when the activation energies of the parallel reac- 
tions are equal. If not, the approximation is 
better the more the principal reaction prevails 
over the side reaction(s). 


4. KINETIC ANALYSIS OF CRACKING OF 


ACETONE 


THE 


It was mentioned in part I how the data, 
presented as they are, permit the prediction of 
the product distribution provided that one of the 
main conversions is known. Extrapolation of the 
conversions from the Figs. 3, 4 and 5, Part I, is 
not allowed, however, as they were obtained with 
temperature profiles too typical for the apparatus 
used, The correct procedure is to extract from 
these data what is general: i.e. the rate equations. 

The data available permit a kinetic analysis of 
reactions (1) and (2) of Part I. The prediction of 
the total acetone conversions would also require, 
strictly speaking an analysis of the acetone 
consuming side reaction for which reaction (3) 
of Part I was shown to be a far too simplified 
scheme. In the range investigated however, the 
total acetone conversion may be related uniquely 
through the ketene selectivity parameter to the 
conversion to ketene (Fig. 6, Part I) so that a 
fairly complete picture is obtained with rate 
equations for reactions (1) and (2) of Part I, only. 

Prior to the derivation of the rate-equations, 
however, the temperature profile of the reaction 
mixture itself had to be calculated. Measuring 
the gas temperature itself would have required an 
inner tube to hold the thermocouples, leaving an 
annular reaction space, which is not likely to be 
adopted for pilot-plant work, a stage that was 
continuously kept in mind in this investigation. 

With the wall temperature and the flow regime 
known and the physical properties of acetone 
and reaction mixtures of different compositions 
taken from the literature or estimated, it is a 
classical problem of heat transfer to calculate the 
temperature of the gases at any point along the 
tube. As the flow regime was laminar, use was 
made of the experimental Sieder—Tate equation, 
represented in a convenient form as 


[3] 


In the reaction section the heat effect of the reac- 
tion was not taken into account. There is ex- 
perimental evidence however that the heat flux 
from outside was so high that the effect of the 
endothermicity of the reaction (18600 cal ‘mol at 
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700-750 C) was easily compensated, so that the 
computed profile is not too far off. In the exit 
section, where the wall temperature decreases 
sharply as a result of the sudden cooling, small 
increments were In Fig. 2 a typical 
temperature profile (experiment 6) measured in 
the wall and the computed profile of the flowing 
reaction gases are shown. Now rate equations for 
reactions (1) and (2) of Part I may be derived. 


che ycnh. 


4 (a). The formation of ketene 

The first step in the kinetic analysis of the 
ketene formation consists of calculating the equiv- 
alent reactor volume. The method proposed in 
section 3 was used to estimate the activation 
energy. In the reaction and exit sections the 
temperature profiles at different reference tem- 
peratures and feed rates were very nearly parallel. 


The largest deviations were observed in the 
preheat section, but their importance is relatively 
minor for the present calculaticn as little reaction 
is taking place in this section. 

In Fig. 1, 7, and p,, are plotted vs. log F for 
two series of experiments, one at 750°, the other 
at 710 °C. The two «,, lines are parallel, while 
although side-reactions do occur, equal p,. corres- 
pond very nearly to equal x,,, as is shown by the 
dotted lines. The conditions for a satisfactory 
estimate of E are fulfilled. 

The the 
parallel x,, lines leads to a unique value of the 
51800 cal mole. This value 
compares quite favourably with the activation 
reactions 


horizontal distance between two 


activation energy : 


energies of a number of cracking 
mentioned by Hixsnetwoop [1] and Murpocn 


and Ho..tanp [4]. 
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Estimation of the activation energy. 


Chem. Engng. Sci. Vol. 13, No. 3. February, 1961. 


183 





G. Froment, H. Pucker and G. Gorruars 








Distance oaiong the tube, cm 


Fic. 2. Wall temperatures, computed temperature profile 
of the gases and equivalent reactor length for experiment 6. 
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Cracking of : 

: 60,500 cal mole 

: 58,500 cal 
47,000 cal 
$5,500 cal 


dipropylether 
diethylether 
methylethylether 


mole (1) 


mole 


acetaldehyde mole 


yropane : 51.300 cal mole 4 
pro} 


With this value of F and the temperature profiles 
of the experiments approximated as described 
above, the equivalent reactor volume may be 
obtained by graphical integration as indicated in 
Fig. 2. With the Vy, F reported in Part I cor- 
rected for each experiment it is possible to draw 
corrected curves for the 7, vs. VV, F data of the 
Figs. 3, 4 and 5 of Part I. The result of this is 
shown in Fig. 3. Compared with Fig. 5 (Part I) 
it is seen how there is no longer any break 
in the curve corresponding to the experiments at 
710 C and how the curves may now be extra- 
polated through the origin. The other conversions 
shown in Figs. 3-5 (Part I) are also better cor- 
related vs. V, F, but they are not shown in 
Fig. 3 as they refer to reactions with activation 
energies which may differ from that used in 
With this corrected 


computing Vi, F. Tak VS. 


V,, F-diagram established, the derivation of a 
rate equation may be tackled. Two ways were 
followed here: the integral and the differential 
method. 
When a rate equation of the form 
k pm (4) 


ac 


is postulated, the continuity equation becomes 


Vr l j dX. 
F kd py” 
0 


(5) 


The differential method is based on (4), the integral 
on (5). 

Integral method, Before the integral in equation 
(5) may be worked out it is necessary to express 
Poo &S a function of x7,,. A rigorous expression 
would only be possible if all reactions taking place 
were exactly known. Therefore an empirical fit 
of the function p,. vs. “7, was undertaken. In 
Fig. 4 it is shown how, at the three temperatures 
considered, py. may be expressed as a unique 
function of x, : 
21ay + 11a (1 


Poo = 1 - 1-05 x)? (6) 





00 











Empirical equation for the p,.. vs. 2, data. 
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Equation (5) then becomes 


1 
for n # > :k= 


moles ('sec' otm” 


Kos 


motes ('sec' atm” 


° 
° 


ae 
°o 


l 
1-05 (VF) (2n — 1) 
l 
(1 — 105 2,,)"""' 


1 
for -:k = — In ——_——-_ (8 
ones F 1 105 2, (8) 
It follows that, when the values of x7,, and V;,/F 
corresponding to the different experiments are 
substituted in equations (7) or (8), k becomes, for 
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5. Determination of the rate-constant, Ky. and of the order, n, of the reaction. 
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Fic. 6. 


each experiment, a function of n only. The point 
of intersection of the k vs. n curves should give 
the value of k at the temperature considered and 
the unique value of n. The k vs. n curves are 
shown in Fig. 5. It is seen how the experiments 
at 750° and 730 °C lead to an order of 1-50 and 
values for k of 0-040 and 0-024 moles 1. sec atm.?® 
The largest deviations amount to 0-005 at 750 
(12 per cent) and 0-0025 moles 1. sec atm'® at 
730 °C (10 per cent). The curves relating to the 
experiments at 710°, which are not so steep, do 
not intersect in one point and do not permit the 
accurate determination of the order at that temp- 
erature. It can only be stated that the relative 
spread on k is lowest between n = 1-0 and 1-5. 
However, if n is also taken to be 1-5 the value 
of 0-013 moles,1. sec. atm'* + 0-001 (8 per cent) is 
obtained for k. When these k values are plotted 
vs. the reciprocal of the absolute temperature, 
1 T, as in Fig. 6, they fall on a straight-line from 
which an activation energy of 52,900 cal /mole is 
derived. It may be concluded that the mechanism 
of formation of ketone from acetone is not altered 
between 710° and 750°C and that the rate of 
this process may be described by 


26,600 
r exp (22-780 - 7 ) pt" (9) 


The value of 52,900 cal, mole for E is very close 


Arrhenius plot for the final determination of the activation energy. 


to that obtained by the short-cut method of 
section 3, Le. that no 
iteration is required. This proves the value of the 
proposed method. 

Differential method. This method is based upon 
equation (4). The rate r, has to be derived from 
ap, VS. Vz/F plot, e.g. by graphical differentia- 
ion. 

The values of k and n are then obtained from 
a logarithmic plot of equation (4) by means of 
a least-square fit of the points to a straight line. 
The results are given in Table 1 together with the 
error variance on 7. 


51,800 cal mole. so 


Table 1 





*c | k 
| 


70 0-041 | 
70 0-023 
710 0-014 





This method, in contrast to the integral method, 
leads to a clearly defined value for n at 710°. 
The value at 730° is less reliable however, because 
due to the deviation of the point related to 
experiment 10 (V, F = 16-5; 2, = 0-233) no- 
ticed in the integral method (Fig. 5), the a vs.Vp/F 
curve cannot be drawn accurately for the higher 
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Vy, F region. Nevertheless, these results confirm 
entirely those obtained by the integral method. 


Kh). The decomposition of keten 


The prediction of the effective ketene content 
at the outlet of the reactor requires a rate equation 
for the decomposition of ketene (reaction 2, 
Part 1) in addition to the equation for the rate 
of formation (equation 9, Part II). 

The procedure described above may be repeated 
here. Again a value of E has to be estimated. 
In an investigation of the ketene-decomposition 
by static experiments (p — 100mm Hg; 510 
590 C) Youne found that at the higher tempera- 
tures the activation energy was about 50,000 
cal mole [5]. From the correlation of the results 
as given in Fig. 6, Part I, it was deduced that the 
different 
from that ef the ketene formation. Therefore E 


51.800 cal mole. 


activation energy was not too much 


was estimated to be which is 
the value used for the determination of the equiv- 
the ketene 
In this way the V, F reported in section 4(a) 


alent reactor volume of formation. 
could be used again as a first estimate, saving 
thus a substantial amount of calculations. 

In Fig. 7 Mya IS plotted Vs. Vp, F at the three 
The 


may be represented by straight lines through the 


reference temperatures investigated. data 
This means that the decomposition is 
respect to ketene. The rate 
then O-O0044, 0-00262 and 0-00153 


origin, 


zero order with 


constants are 


moles |. sec at 750, 730 and 710 °C respectively, 
These values lead to an activation energy of 
51,000 cal mole, which is sufliciently close to the 
estimated value to render a second approximation 
supertluous. This activation energy is also very 
close to that for the ketene formation, as predicted, 
From the standpoint of kinetics, however, no 
detinitive conclusions should be drawn from these 
results. The linear change of x,, with V,/F 
which leads to an order zero, is probably due only 
to the rather narrow range over which the partial 
pressure of ketene was varied. It is a fact however 
that, in the range of practical interest investigated 
in this work, the decomposition of ketene may be 


described with excellent ACCUTACY by 


i moles 


exp (19-617 . 


) (10) 


l. sec 


Substitution of the rate equations (9) and (10) 


into the continuity equation leads to 


26,600 


F da, p,.® 


exp (22-780 


25.660 


T 


exp (19-617 


yan (11) 


where p,. is given by (6). 


In this equation the same value of Vp was used 
for both the formation and decomposition of 
ketene, although this is only an approximation as 


already explained 
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—— raperimental 
---- Caiculoted 








Comparison of e¢ 


When integrated simultaneously with the energy- 
and eventually the momentum equation, the net 
ketene conversion is obtained. From Fig. 6, part 


I, the corresponding total acetone conversion, the 


acetone lost in side-reactions and a fairly good 
estimate of the amount and composition of the 
cracking gas are obtained. 

For non-isothermal 


isothermal conditions or 


conditions which have been corrected by means of 


the equivalent reactor-volume concept, equation 


40 60 70 
Va/F 


xperimental and calculated profiles. 


(11) 


give 


Vie 0-995 1 (2-1 ste 1) ’ 


may be integrated, remembering (7). to 


k, Ve 
27 


(12) 

Fig. 8, finally, is a graphical representation of 
equation (12) which illustrates how the kinetic 
constants derived lead to an excellent reproduction 
of the corrected experimental results. 
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Skewed holding time distributions 


(Reevived 18 June 1900) 


Tur use of the holding or residence time distribution to 
characterize continuous-flow equipment, such as chemical 
reactors, percolators, chromatographic columns, has been 
discussed in the studies of {1} «ned 
KLUINKENBERG and Ssenrrzen [2), 


DANCKWERTS 


DANCKWrETs paper introduces the problem in its 
general form to the chemical engineer, KuNkKeNseRG and 
SIENITZER Investigate under what conditions a normal or 
Gaussian distribution is obtained. It appears that this 
type of distribution is found in many instances where the 
spread in holding times becomes relatively small, such as 
in pipe flow, flow in packed beds, flow through a series of 
identical units. 

This note deals with a few typical cases, where the 
holding 


skewness. 


time distribution may possess a considerable 


Furthermore it will be investigated whether 
the skewness is characteristic of the physical properties 
of the system. 


The skewness, g, is defined by 


» 
* 


2 (t —~ OP RW at, 


where R(t) is the holding time distribution function and 
@ is the average holding time. Standard deviation and 
skewness may be calculated from the Laplace-transformed 
holding time distribution function by a differentiation 
procedure [4]. 

The accompanying table and graph summarize the 
results of the calculations. 

Only for the first case is the skewness, y, proportional 
to the standard deviation, ¢ For all other cases y is 
proportional to o 4/3, at least for low values of «. The 
first distribution is skewed for all values of o, the others 


become symmetrical as o approaches zero. It has been 


Table | 





Relative 


. One completely mixed stage in series 


with a piston flow stage On 


. Equal completely mixed stages in 


series 


2{1 — exp ( 


%. bedely 


flow* 


diffusion in the direction of 


"E 


Sk 


. Mass transfer between a moving and 
a stationary phaset 


standard deviation @ 


= for large np 


JK +F Jz (so) 


Relative skewness y 
0 


Om 


0, 


ny))}'* 


1/3 
+ 3( inp + lpexpt vs) 


31/3 94 3 for large np 


1/3 
of 





*The expression for the standard deviation has been given by van per LAAN [4). 
tSolutions to the underlying differential equations are available from the literature [3). 
The Laplace transformed holding time distribution function for this case reads 


Lp SkKLp 


exp 


u FPutk + {Kp) 


| (p corresponds to f). 
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When the higher values of F /K are left out of considera- 
tion, it appears that between o = 0-5 and o = 1 (1-4 
mixed stages) the differences in relative skewness at equal 
o between the various cases become small. This means 
, that under many practical conditions the skewness is not 
—— ow i very characteristic of the physical process in the system. 

It also follows that in many cases the standard deviation 
(or, alternatively, the number of mixed stages) as a single 
parameter is sufficient to characterize a skewed holding 
time distribution up to the third moment. 


NOTATION 
eddy diffusivity m?/sec 
volume fraction of stationary phase 


ONE STAGE / 7 volume fraction of moving phase 


in SERIES — 4 / 
WiTM PISTON / skewness sec 


F . 
=e mass transfer coeflicient sec™! 
at 

. —— equili- 
concentration in moving phase | prigm 


concentration in stationary phase 


s/f 


a ; 
fe. S/@ “STAGES wn SEMES - height of vessel 
/ number of completely mixed stages 


4s 
fr EDDY OFF USION . : ; 
number of diffusion units 





-. 1 
02 os 





number of transfer units 


wNESS F T : : — is ‘ . 
Sue SS AS A FUNCTION OF THE STANDARD DEVIATION holding time distribution function 


Fic. 1. Skewness as a function of the standard deviation. standard deviation 
time 
linear velocity 
pointed out earlier [2] that the latter distributions become 
Gaussian under this condition. = relative skewness 
For the case of mass transfer, values of ¢ > 1 are 
possible when n, becomes very small. Physically this 
means that the stationary phase acts as a dead zone. This 
leads to tailing of the holding time distribution and 
consequently to large values for ¢ and y. Furthermore the 
skewness increases for increasing capacity F/f/K of the 
stationary phase. 
Shell Oil Company, J.J. vaN Deemrer* 
Houston Research Laboratory 
Houston, Texas, U.S.A. 


average holding time 

average holding time of completely mixed 
stage 

residence time of piston flow stage 


relative standard deviation 


*Present address: Koninklijke Shell-Laboratorium, Amsterdam, Badhuisweg 3, Amsterdam-Noord. 
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Book Reviews 


Principles 


Lesaren A. Hincer and Rourre A. Herser 
of Chemistry. MeGraw-Hill, New York 1960, 


pp. $7.75. 


\V 


Tuts textbook has been written for use in courses empha- 
It 


interest 


sizing the quantitative aspects of modern chemistry. 
is intended primarily for students whose major 
As stated in 


the preface, it is planned for that particular type of student 


is in the physical sciences and engineering 


who will have had an above-average course in high school 


chemistry, will understand and be prepared to use 


first-vear college-level mathematics and will be interested 
in both the “ how ~ and “ why ~ of chemistry. The organ- 
ization and presentation of the subject seems, in general, 
well-suited for the type of student described. Students 
without superior background and enthusiasm for the sub- 
fined the difficult. This 


reviewer believes that, at the present time, the courses 


ject would probably work too 
in mathematics and chemistry taught in many secondary 
schools fail to provide adequate background for the rig- 
orous treatment presented in this book. 

The material covered is divided into two major parts: 
* The 


properties 


“The Properties of Chemical Substances and 
Study of the 


of matter includes emphasis on understanding underlying 


Dynamics of Chemical Change.” 


principles. These principles are developed primarily as 
logical conclusions based on theoretical considerations 
text emphasizes mechanisms by 
which Related 


subjects are grouped together within chapters which are 


The second part of the 


changes in matter are brought about. 
relatively self-contained. 

Following a brief review of simpler fundamental concepts 
the periodic classification of the clements is derived from 
the orderly sequence of atomic numbers. The distribution 
of electrons surrounding atomic nuclei into shells and 
subshells, or orbitals, is deduced by reasoning based on 
the periodic repetition of properties among the clements 
The probabilities of location and the properties of clectrons 
in atomic orbitals are described and explained by four 
quantum numbers defined by « solution of the Schrodinger 
wave equation. 

Electron pairing and the nature of different types ol 
covalent bonds are discussed with emphasis on orbitaf 
theory. This is followed by a chapter dealing with the 
structures of organic compounds and types of reactions 
they exhibit. This section is descriptive in nature, with 
little emphasis on underlying principles. Electron transfer 
and the electrovalent bond, together with the properties 
and reactions of ionic substances are next considered. 

Chapter 5, which follows, is devoted to the discussion of 
a limited number of clementary substances and inorganic 


compounds. The reason for the selection of these particu- 


Pergamon Press Lid., London 
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lar compounds in preference to others is not evident. 
Theoretical considerations are not stressed. 


Chapters 6. 7 
changes in state. Crystal structure and lattice energy are 


and 8 deal with the states of matter and 


emphasized in the discussion of solids. The properties of 
gases are related to randomness of distribution in space 
and the knietic energy of the particles of which gases are 
composed, The liquid state is pictured as intermediate in 
structure and properties. The treatment of the subject 
matter in these chapters is excellent. 

Part 2 — ** The 


in successive chapters with: 


Dynamics of Chemical Change ~ — deals 
the of 


equilibrium ; heterogencous equilibrium and the properties 


nature chemical 


of solutions ; ionic equilibria in aqueous solution ; chemical 


thermodynamics ; oxidation reduction reactions ; chemical 
kinetics and catalysis ; and radiochemistry and the nuclear 
properties of matter, 

The discussion in each of these topics is clear and com- 
prehensive, Underlying principles are utilized in the logical 
development of the subject matter. The student will find 
that sustained concentration is required to secure sound 
understanding of the material presented ; and in a few 
instances the mathematical demands may prove excessive 
for first-year college students, it 


probable that some of the advanced topics, particularly 


This reviewer thinks 


in the chapters dealing with chemical thermodynamics 
and chemical kinetics, will prove to be beyond the grasp 
of many students at this stage of their work. 

Questions and problems at the end of each chapter will 
be helpful to students in the preparation of their work. 
of mathematical and a selected 


\ review operations 


bibliography, grouped according to chapters in the text, 
will be found in the appendix. 

This book 
textbooks available 


In addition to its suitability for use as a text- 


the 
general 


is an important contribution among 


for a rigorous course in 
chemistry. 
book in courses for students with superior background, it 
will be of great value to many teachers in the organization 


of their own courses. 


Sruar’ R. Brinkiey 


Non-Newtonian Fluids. London: 
138 pp. 37s. Gd. 


W. L. WiLkrNson : 
Pergamon Press, Oxford, 1960. 


AN ELEMENTARY textbook on non-Newtonian fluid flow 
iso erdue, and it is pleasing to review a book that does not 
cover ground already adequately covered. 

Because of the immense range of non-Newtonian effects 
and the very real distinctions between various types of 
behaviour, any author must restrict himself to a small 
number of topics: this book is built round the mechanics 


of pipe and channel flow, on the assumption that this is 





Book Reviews 


the facet of greatest importance to chemical engineers. 
On this topic the book is sound, readable and up-to-date 
for both laminar and turbulent flow. 

On other topics, the author is on less sure ground. 
Although his classification of non-Newtonian fluids follows 
conventional lines, it fails to make any reference to the 
necessarily tensorial character of any generally-applicable 
equation of state; for teaching purposes, however, his 
restriction to plane shear will probably be welcomed by 
students in that it avoids undue mathematical complexity. 
The sections devoted to linear models for visco-elastic 
materials (more solid than fluid in their behaviour) prove 
to be of little practical use, and might well have been 
left out 
lack of utility is the fault of engineers who are slow to use 
the mechanical properties of high polymers in design. 


though physicists might argue that this apparent 


The chapters on experimental characterization and on 
viscometric apparatus will probably be useful to those 
who are not familiar with the standard physical measure- 
ments that are commonly made. 

Although well printed and bound, the book is rather 
expensive for the student, who might well expect a more 


comprehensive volume for 37s. 6d. 
J. R.A. Pearson 


C. W. Beckerr, W.S. Benepicr, L. Fano, J. Hinsenraru, 
H. J. Hoge, J. F. Masi, R. L. Nuvrats and Y. 5, 
TovuLovukian : Tables of Thermodynamic and Trans- 
port properties of Air, argon, carbon-dioxide, carbon 
monoxide, hydrogen, nitrogen and Steam. Pergamon 
470. £7. 


Press, London 1960. pp. xiii 


THe present volume is, according to the tly-leaf, a revised 
edition of a series of compilations by the American National 
Bureau of Standards and was originally published by the 
U.S. Department of Commerce as N.B.S. Circular 564. 

The material presented here, which appears to be 
identical in practically all respects with that of Circular 
564, is the compressibility factor, density, specific heat, 
enthalpy, entropy and velocity of sound at pressures of 
0-01, 0-1, 1, 10, 40, 70 and 100 atm from low temperatures 
up to temperatures of 3000-5000 Kk. 


viscosity, thermal conductivity and values of the vapour 


Also given are the 


pressure for all the materials except air, over a wide range 
of temperatures and at close temperature intervals. The 


thermodynamic properties have been obtained using 
spectroscopic data to give the ideal gas state functions 
and the virial equation of state to calculate the real gas 
Tables of ideal gas state thermodynamic 


for further 


functions. 


functions and virial coefficients suitable 
calculations are included. 

The figures presented refer to correlated experimental 
data and have been arrived at by careful weighting of 
experimental measurements. In a number of cases use 
is made of the Lennard-Jones intermolecular potential as 
an interpolation and extrapolation formula, and the best 
The tabulated 


properties of steam have 


experimental results are fitted to this. 
values of the thermodynamic 
been obtained by correlating experimental data with the 
aid of the Keyes equation. The figures given in the Tables 
are expressed as ratios of the properties under standard 
conditions (e.g. 0 °C and 1 atm) and the absolute values 
of the standard quantities are then clearly given in both 
«.g.8. and engineering units, so that a single multiplication 
gives any property in the desired units. 

The Tables, which first appeared in 1955, have been 
found of great value by the reviewer for physical and 
chemical laboratory work and they should be of consider- 
able use in design calculations of all sorts for engineers. 
The production and lay-out of the material is good. 

The original text was sold at 34s. and though this 
presumably represented a subsidised price it is difficult 
to see how the publishers justify an increase to £7 for 
book identical in practically all tabular 
matter with the original. Owners of the earlier volume 
appear to have made a substantial capital gain. 

Although the publishers mention revisions, these are 
difficult to find. 
arises in a table of temperature conversions at the end of 


a re-printed 


However, one revision in particular 


the book which has been altered to bring it into accord 
with the 273-15 K, the previous 
convention in the earlier text being 0 ¢ 273-16 “K. It 
is unfortunate however, that all the enthalpy tables of 
the earlier text have been reproduced in the present work 
divided by RT, where the value of T, = 0 C = 273-16 °K, 

One can only regret that the U.S. Department of Com- 
merce which has given us a number of worthwhile Tables 


convention 0 ¢ 


of thermodynamic data at reasonable cost should have 
found it necessary to discontinue their financial support 


for this valuable set of Tables. 
J. C. McCousrey. 
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Letters to the Editors 


I xerer to the paper “ Transfer from a Sphere into a I have also been told [3] that in his book, Levicu [4] 
Fluid in Laminar Flow ~ by Linton and SurHer.anp [1]. has obtained an “ exact " expression for Npg > © using 
In their paper, the authors refer to an “ inferior approx- a thin boundary layer approximation. His result is 
imation " in my paper [2] to a method used by AKseL’rup. Ngy = 0-097 Npy?! . 

More careful consideration will show that my method, S. K. FriepLanDEeR 
developed independently of AsKEL’RUD’s approximation, 


was actually used to solve a different problem and should ; set ; 
Department of Chemical Engineering 


be judged in that context. My approach permitted the 
The Johns Hopkins University 


determination of the complete Ngyy Npp curve rather than 


the limiting behaviour for Npp —> @, alone. Hence it Baltimore 18 
allowed an estimate of the region where the limiting law Maryland 


begins to break down. U.S.A, 
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